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ABSTRACT 

Continuum strong QCD is the application of models and continuum quantum field theory to the study 
of phenomena in hadronic physics, which includes; e.g., the spectrum of QCD bound states and their in- 
teractions; and the transition to, and properties of, a quark gluon plasma. We provide a contemporary 
perspective, couched primarily in terms of the Dyson- Schwinger equations but also making compar- 
isons with other approaches and models. Our discourse provides a practitioners' guide to features of 
the Dyson-Schwinger equations [such as confinement and dynamical chiral symmetry breaking] and 
canvasses phenomenological applications to hght meson and baryon properties in cold, sparse QCD. 
These provide the foundation for an extension to hot, dense QCD, which is probed via the introduc- 
tion of the intensive thermodynamic variables: chemical potential and temperature. We describe order 
parameters whose evolution signals deconfinement and chiral symmetry restoration, and chronicle their 
use in demarcating the quark gluon plasma phase boundary and characterising the plasma's properties. 
O ■ Hadron traits change in an equilibrated plasma. We exemplify this and discuss putative signals of the 
'effects. Finally, since plasma formation is not an equilibrium process, we discuss recent developments 
in kinetic theory and its application to describing the evolution from a relativistic heavy ion collision 
to an equilibrated quark gluon plasma. 
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1 Prologue 



Continuum strong QCD. The phrase comes to us via Ref. ||l||, although that is likely not its debut. In our 
usage it embraces all continuum nonperturbative methods and models that seek to provide an intuitive 
understanding of strong interaction phenomena, and particularly those where a direct connection with 
QCD can be established, in one true limit or another. The community of continuum strong QCD 
practitioners is a large one and that is helped by the appealing simplicity of models, which facilitates 
insightful contributions that are not labour intensive. 

Our goal is to present an image of the gamut of strong QCD phenomena currently at the forefront of the 
interaction between experiment and theory, and of contemporary methods and their application. Our 
portrayal will not be complete but ought to serve as a conduit into this field. We judge that the best 
means to pursue our aim is to focus on a particular approach, making connections and comparisons with 
others when possible and helpful. All methods have strengths and weaknesses, and, when modelling 
is involved, they are not always obviously systematic. That is the cost of leaving perturbation theory 
behind. They can also yield internally consistent results that nevertheless are without fidelity to QCD. 
Vigilance is therefore necessary but herein lies another benefit of a leagued community. 

Everyone has a favourite tool and the Dyson-Schwinger equations [DSEs] are ours: they provide the 
primary medium for this discourse. The framework is appropriate here because the last decade has 
seen something of a renaissance in its phenomenological application and we chronicle that herein. 
Additionally, the DSEs have been applied simultaneously to phenomena as apparently unconnected 
as low-energy vrvr scattering, B ^ D* decays and the equation of state for a quark gluon plasma, and 
hence they provide a single framework that serves to conduct us through a wide range of hadron physics. 
Focusing on one approach is no impediment to a broad perspective: the continuum, nonperturbative 
studies complement each other, with agreement between them being a strong signal of a real effect. 
In this and their flexibility they provide a foil to the phlegmatic progress of numerical simulations of 
lattice-QCD, and a fleet guide to unanticipated phenomena. 

To close this short introduction we present a list of abbreviations. They are defined in the text when 
first introduced, however, that point can be difficult to locate. 



BSA . 


. Bethe-Salpeter amplitude 


BSE . 


. Bethe-Salpeter equation 


DCSB . 


. dynamical chiral symmetry breaking 


DSE . 


. Dyson-Schwinger equation 


EOS . 


. equation of state 


LHC . 


. large hadron collider 


QCD . 


. quantum chromodynamics 


QC2D . 


. two-colour QCD 


QED . 


. quantum electrodynamics 


QED3 . 


. three-dimensional QED 


QGP . 


. quark gluon plasma 


RHIC . 


. relativistic heavy ion collider 



2 Dyson-Schwinger Equations 

In introducing the DSEs we find useful the brief text book discussions in Chap. 10 of Ref. and Chap. 2 
of Ref. [0], and particularly their use in proving the renormalisability of quantum electrodynamics [QED] 
in Ref. [Q], Chap. 19. However, these sources do not cover the use of DSEs in contemporary nuclear 
and high-energy physics. Reference 0] ameliorates that with a wide ranging review of the theoretical 
and phenomenological applications extant when written, and it provides a good foundation for us. 

The DSEs are a nonperturbative means of analysing a quantum field theory. Derived from a theory's Eu- 
clidean space generating functional, they are an enumerable infinity of coupled integral equations whose 
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solutions are the n-point Schwinger functions [Euclidean Green functions], which are the same matrix 
elements estimated in numerical simulations of lattice-QCD. In theories with elementary fermions, the 
simplest of the DSEs is the gap equation, which is basic to studying dynamical symmetry breaking in 
systems as disparate as ferromagnets, superconductors and QCD. The gap equation is a good example 
because it is familiar and has all the properties that characterise each DSE: its solution is a 2-point 
function [the fermion propagator] while its kernel involves higher n-point functions; e.g., in a gauge the- 
ory, the kernel is constructed from the gauge-boson 2-point function and fermion-gauge-boson vertex, a 
3-point function; a weak-coupling expansion yields all the diagrams of perturbation theory; and solved 
self-consistently, the solution of the gap equation exhibits nonperturbative effects unobtainable at any 
finite order in perturbation theory; e.g, dynamical symmetry breaking. 

The coupling between equations; i.e., the fact that the equation for a given m-point function always 
involves at least one n > m-point function, necessitates a truncation of the tower of DSEs in order 
to define a tractable problem. One systematic and familiar truncation is a weak coupling expansion 
to reproduce perturbation theory. However, that precludes the study of nonperturbative phenomena 
and hence something else is needed for the investigation of strongly interacting systems, bound state 
phenomena and phase transitions. 

In analysing the ferromagnetic transition, the Hartree-Fock approximation yields qualitatively reliable 
information and in QED and QCD its analogue: the rainbow truncation, has proven efficacious. How- 
ever, a priori it can be difficult to judge whether a given truncation will yield reliable results and a 
systematic improvement is not always obvious. It is here that some model-dependence enters but that 
is not new, being typical in the study of strongly-interacting few- and many-body systems. To proceed 
with the DSEs one just employs a truncation and explores its consequences, applying it to different 
systems and constraining it, where possible, by comparisons with experimental data and other theo- 
retical approaches on their common domain. In this way a reliable truncation can be identified and 
then attention paid to understanding the keystone of its success and improving its foundation. This 
pragmatic approach has proven rewarding in strong QCD, as we shall describe. 

2.1 DSE Primer 

Lattice-QCD is defined in Euclidean space because the zero chemical potential [fi = 0] Euclidean QCD 
action defines a probability measure, for which many numerical simulation algorithms are available. 
The Gaussian distribution: 

/C,(g, q') = {2nt)-^/' exp[- (g - qf /(2t)] , (2.1.1) 

defines the simplest probability measure: dq' }Ct{q,q'). K,t is positive and normalisable, which allows 
its interpretation as a probability density. An heuristic exposition of probability measures in quantum 
field theory can be found in Ref. [Q, Chap. 6, while Ref. [0|, Chaps. 3 and 6, provides a more rigorous 
discussion in the context of quantum mechanics and quantum field theory. 

Working in Euclidean space, however, is more than simply pragmatic: Euclidean lattice field theory is 
currently a primary candidate for the rigorous definition of an interacting quantum field theory [^^ 
and that relies on it being possible to define the generating functional via a proper limiting procedure. 
The moments of the measure; i.e., vacuum expectation values of the fields, are the n-point Schwinger 
functions and the quantum field theory is completely determined once all its Schwinger functions are 
known. The time-ordered Green functions of the associated Minkowski space theory can be obtained 
in a well-defined fashion from the Schwinger functions. This is one reason why we employ a Euclidean 
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formulation. Another is a desire to maintain contact with perturbation theory where the renormahsation 
group equations for QCD and their solutions are best understood 0. 

To make clear our conventions: for 4-vectors a, b: 

4 

a-h := a^K5^y ■.= ^ aihi, (2.1.2) 

i=l 

SO that a spacelike vector, Q/i, has > 0; our Dirac matrices are Hermitian and defined by the algebra 

{7^,7^ = 2 5^,; (2.1.3) 

and we use 

75 := -71727374 (2-1-4) 

so that 

tr Y^bl^l'lu1pl<T] = -'^epiupa , £1234 = 1 • (2.1.5) 
The Dirac-like representation of these matrices is: 

where the 2x2 Pauli matrices are: 



Using these conventions the [unrenormalised] Euclidean QCD action is 



S[q,q,A] 



d X 



Nf 



(2.1.8) 



where: = 9^ A" — di,A'J^ — gf°'^'^A''^Al] Nf is the number of quark flavours; ruf are the current-quark 
masses; {A" : a = 1, ... ,8} with [A"^, A''] = 2if"'^'^\'^ are the Gell-Mann matrices for SU{'i) colour; and ^ 
is the covariant gauge fixing parameter. The generating functional follows: 



Z[r],r],J] = / dfi{q,q,A,uj,uj) exp d x qrj + r] q + A"^ 



with sources: 77, 77, J, and a functional integral measure 
d^{q,q,A,u,uj) := 

n n Vq^{x)Vq^{x) n Vuj-{x)Vu\x) \{ VA^^x) exp(-5[g, g, A] - 5^[cu, A]) 



(2.1.9) 



(2.1.10) 



where represents both the flavour and colour index of the quark field, and uj and uj are scalar, 
Grassmann [ghost] fields that are a necessary addition in covariant gauges, and most other gauges 
too. [Without gauge fixing the action is constant along trajectories of gauge-equivalent gluon field 
configurations, which leads to a gauge-orbit volume- divergence in the continuum generating functional.] 
The normalisation 

Z[r/ = 0,r/ = 0, J = 0] = 1 (2.1.11) 

is implicit in the measure. 
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The ghosts only couple directly to the gauge field: 

Sgp,u,A] = J d'x [-d,,ij''d,u;''-gr''d,u^u;'A';^] , (2.1.12) 

and restore unitarity in the subspace of transverse [physical] gauge fields. Practically, the normalisation 
means that ghost fields are unnecessary in the calculation of gauge invariant observables using lattice- 
regularised QCD because the gauge-orbit volume-divergence in the generating functional, associated 
with the uncountable infinity of gauge-equivalent gluon field configurations in the continuum, is rendered 
finite by the simple expedient of only summing over a finite number of configurations. 

It is not necessary to employ the covariant gauge fixing condition in constructing the measure although 
it does yield DSEs with a simple form. Indeed a rigorous definition of the measure may require a gauge 
fixing functional that either completely eliminates Gribov copies or restricts the functional integration 
domain to a subspace without them. Concerned, as they are, with unobservable degrees of freedom, such 
modifications of the measure cannot directly affect the colour-singlet n-point functions that describe 
physical observables. However, they have the capacity to modify the infrared behaviour of coloured n- 
point functions P,p^ and thereby influence the manner in which we intuitively understand observable 
effects. 

The DSEs are derived from the generating functional using the elementary observation that, with 
sensible Euclidean space boundary conditions, the integral of a total derivative is zero: 

= / expy'ci^a;[gr/ + r/g + J^A;:] . (2.1.13) 

Examples of such derivations: the QED gap equation, the equation for the photon vacuum polarisa- 
tion and that for the fermion-fermion scattering matrix [a 4-point function], are presented in Ref. 0, 
Chap. 10. The scattering matrix is important because it lies at the heart of two-body bound state 
studies in quantum fleld theory, being a keystone of the Bethe-Salpeter equation. Section 2.1 of Ref. 
repeats the flrst two derivations and also makes explicit the effects of renormalisation. Herein we simply 
begin with the relevant DSE leaving its derivation as an exercise. That too can be side-stepped if a 
Minkowski space version of the desired DSE is at hand. Then one need merely employ the transcription 
rules: 

Conflguration Space Momentum Space 

1. / d^x'' ^ -I dS"" 1. / d^k"" ^ I / d^k"" 



2. p ^ i-f^ -d^ 2. Ife ^ -i-f^ ■ k^ 

3. 4 ^ -t-f^ -A^ 3. 4 ^ -i-i^ ■ A^ 

4. A^B^' -A^ ■ 4. k^q^' -k^ ■ 

5. x^'^^ ^x^ -d^ 5. fc^x^ ^ -k^ ■ x^ 



These rules are valid in perturbation theory; i.e., the correct Minkowski space integral for a given 
diagram will be obtained by applying these rules to the Euclidean integral: they take account of the 
change of variables and rotation of the contour. However, for the diagrams that represent DSEs, which 
involve dressed n-point functions whose analytic structure is not known a priori, the Minkowski space 
equation obtained using this prescription will have the right appearance but it's solutions may bear no 
relation to the analytic continuation of the solution of the Euclidean equation. 
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2.2 Core Issues 



The use of DSEs as a unifying, phenomenological tool in QCD has grown much since the pubhcation of 
Ref . , and Ref . [|12] reviews apphcations to electromagnetic and hadronic interactions of light mesons 
and the connection to the successful Nambu-Jona-Lasinio |jl2|,|T3|,|l4| and Global-Colour models. In 
most of their widespread applications these two models yield results kindred to those obtained with the 
rainbow-ladder DSE truncation. 



Quark DSE [Gap Equation] The idea that the fermion DSE, or gap equation, can be used to study 



dynamical chiral symmetry breaking [DCSB] perhaps began with a study of the electron propagator [|I^ . 
Since that analysis there have been many studies of dynamical mass generation in QED, and those re- 
viewed in Ref. |^ establish the context: there is some simplicity in dealing with an Abelian theory. 
Continuing research has led to an incipient understanding of the role played by multiplicative renor- 
malisability and gauge covariance in constraining the dressed- fermion-photon vertex fl^, and progress 
in this direction provides intuitive guidance for moving beyond the rainbow truncation. However, the 
studies have also made plain the difficulty in defining the chiral limit of a theory without asymptotic 
freedom. Renormalised, strong coupling, quenched QED yields a scalar self energy for the electron that 
is not positive-definite: damped oscillations appear after the renormalisation point |jl8|. [Quenched 
= bare photon propagator in the gap equation. It is the simplest and most widely explored trunca- 
tion of the gauge sector.] This pathology can be interpreted as a signal that four fermion operators 
~ ['i/j{x)ip{x)]^ have acquired a large anomalous dimension and have thus become relevant operators for 
the range of gauge couplings that support DCSB. That hypothesis has implications for the triviality of 
the theory, and is reviewed and explored in Refs. |]5|JT9[]. Nonperturbative QED remains an instructive 
challenge. 

The study of DCSB in QCD is much simpler primarily because the chiral limit is well-defined. We 
discuss that now using the renormalised quark-DSE: 

Sip)-' = Z2{il-p + mbare) + ^1 9^ D ,,{p - q)^^ {q)Tl{q, p) , (2.2.1) 

where D^^{k) is the renormalised dressed-gluon propagator, T1{q;p) is the renormalised dressed-quark- 
gluon vertex, mbare is the A-dependent current-quark bare mass that appears in the Lagrangian and 
:= d^q/ {2tiY represents mnemonically a translationally-invariant regularisation of the integral, 
with A the regularisation mass-scale. The final stage of any calculation is to remove the regularisation by 
taking the limit A — >■ oo. Using a translationally invariant regularisation makes possible the preservation 
of Ward-Takahashi identities, which is crucial in studying DCSB, for example. One implementation 
well-suited to a nonperturbative solution of the DSE is Pauli-Villars regularisation, which has the 
quark interacting with an additional massive gluon-like vector boson: mass ~ A, that decouples as 
A — i> oo [^. An alternative is a numerical implementation of dimensional regularisation, which. 



although more cumbersome, can provide the necessary check of scheme-independence [^ 



In Eq. (|2.2.1[ ), Zi{('^,A'^) and Z2(C^, A^) are the quark-gluon-vertex and quark wave function renormal- 



isation constants, which depend on the renormalisation point, and the regularisation mass-scale, as 
does the mass renormalisation constant 

ZUe, A') = Z.iC, A')/Z,{e, A'), (2.2.2) 

with the renormalised mass given by 

m(C) := mbare(A)/^™(C', A')- (2.2.3) 
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Although we have suppressed the flavour label, S, and mbare depend on it. However, one can 
always use a flavour-independent renormalisation scheme, which we assume herein, and hence all the 
renormalisation constants are flavour-independent p0|j2^ . 



The solution of Eq. (|2.2.1|) has the form 

sip)-' = t^-pA{p',e)+B{p',e) 



z{p\C 



(2.2.4) 



The functions A{p'^,('^): BijP'X'^) embody the effects of vector and scalar quark-dressing induced by 
the quark's interaction with its own gluon fleld. The ratio: M{p'^,C,'^), is the quark mass function and 
a pole mass would be the solution of 



m: 



pole 



M^p" 



"^pole,C')=0. 



(2.2.5) 



A widely posed conjecture is that conflnement rules out a solution of this equation |2^. We discuss this 
further below. 

Equation ( [2.2.1 ) must be solved subject to a renormalisation [boundary] condition, and because the 
theory is asymptotically free it is practical and useful to impose the requirement that at a large spacelike 



■ p + m{() 



(2.2.6) 



where m{() is the renormalised current-quark mass at the scale (. By "large" here we mean ("^ ^ ^qcd 
so that in quantitative, model studies extensive use can be made of matching with the results of 
perturbation theory. It is the ultraviolet stability of QCD; i.e., the fact that perturbation theory is 
valid at large spacelike momenta, that makes possible a straightforward definition of the chiral limit. It 
also provides the starkest contrast to strong coupling QED. 



Multiplicative renormalisability in gauge theories entails that 



(2.2.7) 



and beginning with Ref. PI this relat ion has been used efficaciously to build realistic Ansdtze for the 
fermion-photon vertex in quenched QED. A systematic approach to such nonperturbative improvements 
is developing [^,^] and these improvements continue to provide intuitive guidance in QED, where they 
complement the perturbative calculation of the vertex They are also useful in exploring model 

dependence in QCD studies. 



At one loop in QCD perturbation theory 



-7f//3i 



^2(C',A') 

and at this order the running strong coupling is 



If 



TT 



2 ' 



- In [CVA^cD 



(2.2.8) 



(2.2.9) 



In Landau gauge: ^ = 0, so = 1 at one loop order. This, plus the fact that Landau gauge is a fixed 
point of the renormalisation group [Eq. ( |2.2.2(J| )], makes it the most useful covariant gauge for model 
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studies. It also underlies the quantitative accuracy of Landau gauge, rainbow truncation estimates of 
the critical coupling in strong QED |lT7| , p8| . In a self consistent solution of Eq. (|2.2.1|) , Z2 7^ 1 even 
in Landau gauge but, at large C^, the ^-dependence is very weak. However, as will become evident, in 
studies of realistic QCD models this dependence becomes significant for < 1-2 GeV^, and is driven 
by the same effect that causes DCSB. 

The dressed-quark mass function: M(p^, ('^) = B{p'^^ ^^), is independent of the renormalisation 

point; i.e., with C 7^ C 

M{p\ e) = M{p\ C') ■■= Mip') , : (2.2.10) 

it is a function only of p'^/^qcbj which is another constraint on models. At one loop order the running 
[or renormalised] mass 



m(C) = MiC' 



m 



7m 5 7m 



12/(33 - 2Nf) 



(2.2.11) 



iln 



QCD 



where rh is the renormalisation point independent current-quark mass, and the mass renormalisation 
constant is, Eq. ( [^.2. 21 ), 



^m(C',A' 



a(A^) 



(2.2.12) 



The mass anomalous dimension, 7^, is independent of the gauge parameter to all orders perturbation 
theory and for two different quark flavours the ratio: '^/i(C)/'^/2(C) = ^hl^f2i which is independent 
of the renormalisation point and of the renormalisation scheme. The chiral limit is unambiguously 
defined by 

chiral limit : m = . (2.2.13) 

In this case there is no perturbative contribution to the scalar piece of the quark self energy; i.e., 
-S(p^?C^) = at every order in perturbation theory and in fact there is no scalar mass-like divergence 
in the calculation of the self energy. This is manifest in the quark DSE, with Eq. (|2.2.1| ) yielding, in 
addition to the perturbative result: B{p^X'^) = 0, a solution M{p^) = B{p^ X'^) / ^^ip^ X"^) that 
is power-law suppressed in the ultraviolet: MijP') ~ 1/p^, guaranteeing convergence of the associated 
integral without subtraction. Dynamical chiral symmetry breaking is 



DCSB : M(p2) ^ when m = . 



(2.2.14) 



As we shall see, in QCD this is possible if and only if the quark condensate is nonzero. The criteria are 
equivalent. 

The solution of the quark DSE depends on the anatomy of the dressed-gluon propagator, which in con- 
cert with the dressed- quark-gluon vertex encodes in Eq. ( p.2.1|) all effects of the quark-quark interaction. 
In a covariant gauge the renormalised propagator is 



k^ 



k^kiy 

1?~ 



(2.2.15) 



where d{k'^, C^) = l/[l + n(A;^, C"^)], with n(fc^, ('^) the renormalised gluon vacuum polarisation for which 
the conventional renormalisation condition is 



^(c^C') = 0; i.e.,ci(c^C') = l 



(2.2.16) 
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For the dressed gluon propagator, multiplicative renormalisability entails 



- ^ -d{e,e)= , (2.2.17) 



and at one loop in perturbation theory 

-71 /A 



a(A2 



7i = iiV;-i(13-30. (2.2.18) 



The gauge parameter is also renormalisation point dependent; i.e., the renormalised theory has a running 
gauge parameter. However, because of Becchi-Rouet-Stora [BRST or gauge] invariance, there is no new 
dynamical information in that: its evolution is completely determined by the gluon wave function 
renormalisation constant 

e(c') = z3-l(c^A')ebare(A). (2.2.19) 

One can express ^(C^) in terms of a renormalisation point invariant gauge parameter: ^, which is an 
overall multiplicative factor in the formula and hence 

Landau Gauge : = ^ ^(C^) = (2.2.20) 

at all orders in perturbation theory; i.e.. Landau gauge is a fixed point of the renormalisation group. 

The renormalised dressed-quark-gluon vertex has the form 

Tt{k,p) = ^T,{k,p); (2.2.21) 

i.e., the colour matrix structure factorises. It is a fully amputated vertex, which means all the ana- 
lytic structure associated with elementary excitations has been eliminated. To discuss this further we 
introduce the notion of a particle-like singularity. It is one of the form: P = k — p, 

(^^^'^e(0'l]- (2-2.22) 

If the vertex possesses such a singularity then it can be expressed in terms of non-negative spectral 
densities, which is impossible if x > 1. x = 1 is the ideal case of an isolated 5-function distribution in 
the spectral densities and hence an isolated free-particle pole, x G (0, 1) corresponds to an accumulation 
at the particle pole of branch points associated with multiparticle production, as occurs with the electron 
propagator in QED because of photon dressing. 

The vertex is a fully amputated 3-point function. Hence, the presence of such a singularity entails 
the existence of a flavour singlet composite (quark-antiquark bound state) with colour octet quantum 
numbers and mass m = h. The bound state amplitude follows immediately from the associated homo- 
geneous Bethe-Salpeter equation, which the singularity allows one to derive. Such an excitation must 
not exist as an asymptotic state, that would violate the observational evidence of confinement, and any 
modelling of r^(A;,p) ought to be consistent with this. 

Expressing the Dirac structure of Vy{k^p) requires 12 independent scalar functions: 

r,(fc,p)=7,Fi(fc,p,C) + ... , (2.2.23) 

which for our purposes it is not necessary to reproduce fully. A pedagogical discussion of the perturbative 
calculation of Vy{k,p) can be found in Ref. while Refs. explore its nonperturbative structure 
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and properties. We only make Fi{k,p,() explicit because the renormalisability of QCD entails that it 
alone is ultraviolet divergent. Defining 



/l(A;^C') ■■=F,{k,-k,C) 
the conventional renormalisation boundary condition is 



which is practical because QCD is asymptotically free. Multiplicative renormalisability entails 



^2.2.24) 



^2.2.25) 



(2.2.26) 



and at one loop order 



-IT /Pi 



(2.2.27) 



Dynamical Chiral Symmetry Breaking At this point each element in the quark DSE, Eq. (2.2.1), is 



defined, with some of their perturbative properties elucidated, and the question is how does that provide 
an understanding of DCSB? It is best answered using an example, in which the mo del- independent 
aspects are made clear. 

The quark DSE is an integral equation and hence its elements must be known at all values of their 
momentum arguments, not just in the perturbative domain but also in the infrared. While the gluon 
propagator and quark-gluon vertex each satisfy their own DSE, that couples the quark DSE to other 
members of the tower of equations and hinders rather than helps in solving the gap equation. Therefore 
here, as with all applications of the gap equation, one employs Ansdtze for the interaction elements 
[D^^{k) and r,^(A;,p)], constrained as much and on as large a domain as possible. This approach has a 
long history in exploring QCD and we illustrate it using the model of Ref. ||2^ . 



The renormalised dressed-ladder truncation of the quark-antiquark scattering kernel [4-point function] 
is 



K{p,q-P)Z = g\e)D,,{p~q) Vl{p+,q+) S{q^) ^ 5(g_)C(g_,p 



tr 



(2.2.28) 



where p± = p ± P/2, g-t = g ± P/2, with P the total momentum of the quark-antiquark pair, 
and although we use it now we have suppressed the ^-dependence of the Schwinger functions. From 
Eqs. ( p.2.16[ - p.2.18|) it follows that for := {p — qY large and spacelike 



^3(c^A^) 



71 //3i 



71 //3i 



DlT{p-q)- (2.2.29) 



Using this and analogous results for the other Schwinger functions then on the kinematic domain for 
which ~ ~ q^ is large and spacelike [^'^(C^) '■= 47ra((^^)] 



K{p, q- P)Z ^ AMQ'') D'^:%p - q) Wl, 5'-^(g+) ^ S*-^(g-) W% 



Tfrcc/ 



tr 



ifrcc / 



because Eqs. (|]23), ( p.2.18|) and (|2.2.27| ) yield 



2 7f , 7i 
/3i Pi 



27r 



1. 



(2.2.30) 



^2.2.31) 
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This is one way of understanding the origin of an often used Ansatz in studies of the gap equation; i.e., 
making the replacement 

9'D,,{k) -> ATraik') D^^^k) (2.2.32) 
in Eq. ( 2.2.1| ), and using the "rainbow truncation:" 



(2.2.33) 



Equation ( p. 2. 32 ) is often described as the "Abehan approximation" because the left- and right-hand- 
sides [r.h.s.] are equal in QED. In QCD, equality between the two sides cannot be obtained easily by 
a selective resummation of diagrams. As reviewed in Ref. Eqs. (5.1-5.8), it can only be achieved 
by enforcing equality between the renormalisation constants for the ghost-gluon vertex and ghost wave 
function: ^1 = ^3. A mutually consistent constraint, which follows formally from Zi = Z^, is to enforce 
the Abelian Ward identity: Zi = Z^- At one- loop this corresponds to neglecting the contribution 
of the 3-gluon vertex to F^, in which case 7r — *■ f^C = 7f- This additional constraint provides the 
basis for extensions of Eq. ( |2.2.33| ); i.e., using Ansdtze for that are consistent with the QED vector 



Ward-Takahashi identity; e.g., Ref. 



Arguments such as these inspire the following Ansatz for the kernel in Eq. (|2.2.1|) | 

Z, g'D,^{p - q)^^,S{q)Tl{q,p) ^ £ g{{p - qf) D^l^p - q)^^^S[q)^^^ , (2.2.34) 



with the ultraviolet behaviour of the "effective coupling:" ^(fc^), fixed by that of the strong running 
coupling. Since it is not possible to calculate Zi nonperturbatively without analysing the DSE for the 
dressed-quark-gluon vertex, this Ansatz absorbs it in the model effective coupling. 



Equation (|2.2.34| ) is a model for the product of the dressed-propagator and dressed-vertex and its 
definition is complete once the behaviour of ^(fc^) in the infrared is specified; i.e., for fc^ < 1-2 GeV^. 
Reference used 



iln 



(1 + F/A 



2 

QCD 



(2.2.35) 



with J^{k^) = [l-exp(-A;V[4m2])]/A;2 and r = e^-l. For Nf = 4, Aq^^ = 0.234 GeV. The qualitative 
features of Eq. ( |2.2.35| ) are plain. The first term is an integrable infrared singularity |^ and the second 
is a finite- width approximation to 5^(fc), normalised such that it has the same / d'^k as the first term. 
In this way the infrared strength is split into the sum of a zero-width and a finite-width piece. The 
last term in Eq. ( |2.2.35| ) is proportional to a{k'^)/k'^ at large spacelike-A;^ and has no singularity on the 
real-fc^ axis. 



There are ostensibly three parameters in Eq. ( p.2.35|) : D, u and m^. However, in Ref. |2^ the au- 
thors fixed u = 0.3GeV(= l/[.66fm]) and = 0.5GeV(= l/[.39fm]), and only varied D and the 
renormalised u = d- and s-current-quark masses in an attempt to obtain a good description of low- 
energy vr- and ii"-meson properties, using a renormalisation point C = 19 GeV that is large enough to 
be in the perturbative domain. [The numerical values of uj and mt are chosen so as to ensure that 
g{k'^) ~ 47ra(/c^) for fc^ > 2 GeV^. Minor variations in u and mt can be compensated by small changes 
in D.] Such a procedure could self-consistently yield D = 0, which would indicate that agreement with 
observable phenomena precludes an infrared enhancement in the effective interaction. However, that 
was not the case and a good fit required 



D = (0.884 GeV)^ 



(2.2.36) 
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with renormalised current-quark masses 

mu,d{C) = 3.74 MeV , m,(C) = 82.5 MeV , 
which are in the ratio 1 : 22, and yielded, in MeV, 



^2.2.37) 
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(2.2.38) 



and other quantities to be described below. An explanation of how this fit was accomplished requires a 
discussion of the homogeneous Bethe-Salpeter equation, which we postpone. Here we instead focus on 
describing the properties of the DSE solution obtained with these parameter values. 



(2.2.39) 



Using Eqs. (|2.2.1[ - ^l273|) and ( p.2.34|) the gap equation can be written 

SipXr' = Z2^7-P + ^4m(C) + S'(p,A), 
with the regularised quark self energy 

S'(P,A) := jy{{p-q)')D^-{p-q)^^,S{q)^^^. (2.2.40) 

When m 7^ the renormalisation condition, Eq. ( |2.2.6| ), is straightforward to implement. Writing 

S'(p, A) := Z7 ■ P {A'ip', A') - l) + B\p', A') , (2.2.41) 
which emphasises that these functions depend on the regularisation mass-scale. A, Eq. ( p.2.6|) entails 
Z^iC, A') = 2- A'(C^ A^) and m(C) = Z2(C', A') ^^^.(A^) + i?'(C^ A^) (2.2.42) 

so that 

A{p\ (') = ! + A{p\ A') - A'(C^ A^) , B{p\ C) = m{C) + B\p\ A') - B'{C\ A') . (2.2.43) 

Multiplicative renormalisability requires that having fixed the solutions at a single renormalisation 
point, (, their form at another point, (, is given by 



S~\p, = ^l -pAip^ + B{p\ C') = lS^fl s-\p, . 



(2.2.44) 



This feature is evident in the solutions obtained in Ref. It means that, in evolving the renormali- 
sation point to C, the "1" in Eqs. {^^2^ is replaced by Z2(C^ A2)/Z2(C^ A^), and the "m(C)" by m(C); 
i.e., the "seeds" in the integral equation evolve according to the QCD renormalisation group. This is 
why Eq. (|2.2.34|) is called a "renormalisation-group-improved rainbow truncation." 

Turning to the chiral hmit, it follows from Eqs. (pIT^) , (^2^, ( p.2.11|) and ( p.2.13|) that for m = 

Z2(C', A^) mbarc(A') = , VA . (2.2.45) 

Hence, as remarked on page ^ there is no subtraction in the equation for B{p'^,(^)j Eq. ( |2.2.43| ) 
becomes 

.2 a2\ _ D'/^2 a2\ 1- D//_2 A 2n 



B{p',C) = B'{p\A') , hm B'{p',A') < oo , 



^2.2.46) 
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Fig. 2.1. Quark mass function obtained as a solution of Eq. ( |2.2.1|) using the model of 
Eqs. ( |2.2.34| ), ( |2.2.35| ), and current-quark masses, fixed at C = 19 GeV: m^^ = 3.7 MeV, mi 



82 MeV, 



0.58 GeV and = 3.8 GeV. The indicated solutions of M'^{p'^) = define 



the Euclidean constituent-quark mass, Mf in Eq. ( |2.2.56|) , which takes the values: 
0.56 GeV, Mf = 0.70 GeV, Mf = 1.3 GeV, = 4.6 GeV. (Adapted from Refs. 



which is only possible if the mass function is at least 1/p^-suppressed. This is not the case in quenched 
strong coupling QED, where the mass function behaves as ~ cos(const. ln[p^/C^])/(p^/C^)^''^ 
and that is the origin of the complications indicated on page |^ |]5|, p!9| , |21| . 



In Fig. |2]T| we present the renormalised dressed-quark mass function, M(p^), obtained by solving 
Eq. ( |2.2.39| ) using the model and parameter values of Ref. [^, Eqs. ( p.2.34| - p.2.37D , and also in the chi- 
ral limit and with typical heavy-quark current-mass values. In the presence of explicit chiral symmetry 
breaking Eq. ( [^.2.11 ) describes the form of M(p^) for > 0(1 GeV^). In the chiral limit, however, the 
ultraviolet behaviour is given by 



largc-p2 2tT 7„ 



(iln 



QCD 



l-7m ' 



(2.2.47) 



where {qqY is the renormalisation-point-independent vacuum quark condensate. This behaviour too is 
characteristic of the QCD renormalisation group and exhibits the power-law suppression anticipated 
on page These results for the large p^ behaviour of the mass function are model independent; i.e., 
they arise only because the DSE truncation is consistent with the QCD renormalisation group at one 
loop. (It has long been known that the truncation defined by Eq. ( |2.2.34|) yields results in agreement 
with the QCD renormalisation group at one loop; e.g., Refs. 



derived in Ref. |^ 



The gauge invariant expression for the renormalisation-point-dependent vacuum quark condensate was 

- ml ■■= UC\ A') iVctro r S°(g, C) , (2.2.48) 

J a 
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where tiD identifies a trace over Dirac indices only and the superscript "0" indicates the quantity was 
calculated in the chiral limit. Substituting Eq. ( |2.2.47| ) into Eq. (|2.2.48|) , recalling that Z4 = in 
Landau gauge and using Eq. ( |2.2.12| ), yields the one- loop expression 



(gg)J=(ilnCVA|cDj {mf ■ 
Employing Eq. (|2.2.11|) , this exemplifies the general result that 



(2.2.49) 



(2.2.50) 



i.e., that this product is renormalisation point invariant and, importantly, it shows that the behaviour 
expressed in Eq. ( p. 2. 471) is exactly that required for consistency with the gauge invariant expression 
for the quark condensate. A model, such as Ref. in which the scalar projection of the chiral limit 
dressed-quark propagator falls faster than up to In-corrections, is only consistent with this quark 

condensate vanishing, and it is this condensate that appears in the current algebra expression for the 
pion mass 



42 



Equation (|2.2.47|) provides a reliable means of calculating the quark condensate because corrections are 
suppressed by powers of K^q-q/C,'^ . Analysing the asymptotic form of the numerical solution one finds 



- (gg)0 = (0.227 GeV)^ 
Using Eq. ( |2.2.49| ) one can define a one-loop evolved condensate 



- (In [1/Aqcd])"™ = (0.241 GeV)^ . 



C=l GeV 



(2.2.51) 



(2.2.52) 



This can be directly compared with the value of the quark condensate employed in contemporary 
phenomenological studies Q: (0.236 ± 0.008 GeV)^. The authors of Ref. [||] noted that mcreasmg 
uj ^ l.buj in Qik"^) increases the calculated value in Eq. ( p.2.52|) by ~ 10%; i.e., the magnitude of the 
condensate is correlated with the degree of infrared enhancement /strength in the effective interaction. 
That is unsurprising because it has long been known that there is a critical coupling for DCSB; i.e., the 
kernel in the gap equation must have an integrated strength that exceeds some critical value [40 . This 
is true in all fermion-based studies of DCSB. 



The renormalisation-point-invariant current-quark masses corresponding to the mf[Q in Fig. ^ 



are 

obtained in the following way: using Eq. ( 2.2.48| ), direct calculation from the chiral limit numerical 
solution gives 

{m)U^OeY = -(0.275 GeV)^ (2.2.53) 
and hence from the values of mj = rufiC) listed in Fig. ^ and Eqs. ( pXSOD , (|2:23l| ), in MeV, 

m„,rf = 6.60, = 147, = 1 030 , rhf, = 6 760 , (2.2.54) 

from which also follow one- loop evolved values in analogy with Eq. ( |2.2.52| ): 

ml GeV ^ 5 5 ^ ^IGeV ^ ^gg ^ ^1 GeV ^ ggg ^ ^1 GeV = 5 70O . (2.2.55) 



Figure |2T^ highlights a number of qualitative aspects of the quark mass function. One is the difference 
in the ultraviolet between the behaviour of M(p^) in the chiral limit and in the presence of explicit 
chiral symmetry breaking. In the infrared, however, the u, d-quark mass function and the chiral limit 
solution are almost indistinguishable. The chiral limit solution is nonzero only because of the nonpertur- 
bative DCSB mechanism whereas the u, d-quark mass function is purely perturbative at > 20 GeV^. 
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Hence the evolution to coincidence of the chiral-hmit and m, (i-quark mass functions makes clear the 
transition from the perturbative to the nonperturbative domain. It is on this nonperturbative domain 
that 74(p^,(^^) differs significantly from one. [This behaviour and that of the light-quark mass func- 



tion depicted in Fig. |2.1| have recently been confirmed in numerical simulations of lattice-QCD 
A concomitant observation is that the DCSB mechanism has a significant effect on the propagation 
characteristics of u, d, s-quarks. However, as evident in the figure, that is not the case for the 6-quark. 
Its large current-quark mass almost entirely suppresses momentum-dependent dressing, so that M{,(p^) 
is nearly constant on a substantial domain. This is true to a lesser extent for the c-quark. 



To quantify the effect of the DCSB mechanism on massive quark propagation characteristics the authors 
of Ref. P7| , |5B| introduced a single measure: := M^/mj, where Mj is the Euclidean constituent- 
quark mass, defined as the solution of 

(Mf )2 - M2(p2 = (M/)2, C') = . (2.2.56) 
In this exemplifying model the Euclidean constituent-quark mass takes the values listed in the figure. 



which have magnitudes and ratios consistent with contemporary phenomenology; e.g., Refs. and 

(2.2.57) 
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These values are representative and definitive: for light-quarks Cg=u,d,s ~ 10-100, while for heavy-quarks 
CQ=c,b ~ 1, and highlight the existence of a mass-scale characteristic of DCSB: M^. The propagation 
characteristics of a flavour with < are significantly altered by the DCSB mechanism, while 
for flavours with ^ momentum-dependent dressing is almost irrelevant. It is apparent and 
unsurprising that ~ 0.2GeV~ Aqcd- This feature of the dressed-quark mass function provides 
the foundation for a constituent-quark-like approximation in the treatment of heavy-meson decays and 



transition form factors p8 |. 



To recapitulate. The quark DSE describes the phenomena of DCSB and the concomitant dynamical 
generation of a momentum dependent quark mass function, with the renormalisation-group-improved 
rainbow-truncation yielding model-independent results for the momentum dependence of the mass func- 
tion in the ultraviolet. For light quarks, defined by ^ 1, the magnitude of the mass function in the 
infrared; i.e., for < 1-2 GeV^, is determined by the behaviour of the effective quark-quark interaction 
on the same domain, and so is the value of the vacuum quark condensate. An infrared enhancement 
in this effective interaction is required to describe observable light-meson phenomena. We have only 
provided a single illustration but it is supported by, e.g., Refs. p7| , |48| , ^ , |50| , |5ll , |52| and the observation 
that with insufficient infrared integrated strength in the kernel of Eq. ( |2.2.1| ) the quark condensate 
vanishes [0,0,^,^j5^, which is a poor starting point for light-hadron phenomenology. The question 



now arises, where does this strength come from? 



Gluon DSE Guidance here comes from studies of the DSE satisfied by the dressed-gluon propagator. 



which is depicted in Fig. p7^ . However, as we now describe, these studies are inconclusive. 



Early analyses |^ used the ghost-free axial gauge: n ■ A" = 0, > 1, in which case the second 



equation in the figure is absent and two independent scalar functions: Fi, F2, are required to fully 
specify the dressed-gluon propagator, cf. the covariant gauge expression in Eq. (|2.2.15|) , which requires 
only one function. In the absence of interactions: -Fi(fc^) = — l//c^, F2{k'^) = 0. These studies employed 
an Ansatz for the three-gluon vertex that doesn't possess a particle-like singularity and neglected the 
coupling to the quark DSE. They also assumed F2 = 0, even nonperturbatively, and ignored it in solving 
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Fig. 2.2. From top to bottom, depictions of the DSEs for the gluon (spring), ghost (dashed- 
hne) and quark (sohd-hne) 2-point functions. Following convention, a filled circle denotes a 
fully dressed propagator and an open circle, a one-particle irreducible vertex; e.g., the open 
circle in the first line represents the dressed-three-gluon vertex. The figure illustrates the 
interrelation between elements in the tower of DSEs: the gluon propagator appears in the 
DSE for the quark and ghost propagator; the ghost and quark propagator in the DSE for 
the gluon, etc. (Adapted from Ref. p8| .) 



the DSE. The analysis then yielded 



Fi(A;')'~°^; (2.2.58) 



i.e., a marked infrared enhancement that can yield an area law for the Wilson loop |^| and hence 
confinement, and DCSB as described above without fine-tuning. This effect is driven by the gluon 
vacuum polarisation, diagram three in the first line of Fig. (p.2|) . A similar result was obtained in 
Ref. |61[]. However, a possible flaw in these analyses was identified in Ref. |]62[, which argued from 



properties of the spectral density in ghost-free gauges that F2 cannot be zero but acts to cancel the 
enhancement in Fi. [Preserving F2 yields a coupled system of equations for the gluon propagator that 
is at least as complicated as that obtained in covariant gauges, which perhaps outweighs the apparent 
benefit of eliminating ghost fields in the first place.] 



There have also been analyses of the gluon DSE using Landau gauge and those of Refs. |]63| , |6^ , |65| , |66| , |67 



are unanimous in arriving at the covariant gauge analogue of Eq. ( p.2.58|) , again driven by the gluon 



vacuum polarisation diagram. In these studies Ansdtze were used for the dressed-three-gluon vertex, all 
of which were free of particle-like singularities. However, these studies too have weaknesses: based on an 
anticipated dominance of the gluon-vacuum polarisation, truncations were implemented so that only the 
third and fifth diagrams on the r.h.s. of the first equation in Fig. ( |2.2|) were retained. In covariant gauges 
there is a priori no reason to neglect the ghost loop contribution, diagram six, although perturbatively 



its contribution is estimated to be small [37 
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Another class of Landau gauge studies are described in Refs. p8| , |69| , which propose solving the DSEs 
via rational polynomial Ansdtze for the one-particle irreducible components of the Schwinger functions 
appearing in Fig. |2.2| ; i.e., the self energies and vertices. This method attempts to preserve aspects of 
the organising principle of perturbation theory in truncating the DSEs. In concrete calculations, for 
simplicity, only the first, third and sixth diagrams on the r.h.s. of the first equation in Fig. ^.2| survive, 
the last [fermion] equation is neglected and the leading order solution of the ghost equation has the 
appearance of the massless free propagator: ~ The analysis suggests that a consistent solution 

for the dressed-gluon propagator is one that vanishes at fc^ = 0; i.e., in Eq. ( |2.2.15| ) 



k^ + -f^ 



(2.2.59) 



However, the associated polynomial Ansatz for the dressed-three-gluon vertex exhibits particle-like 
singularities and this is characteristic of the method. The question of how this can be made consistent 
with the absence of coloured bound states in the strong interaction spectrum is currently unanswered. 



Proponents of the result in Eq. ( p.2.59|) claim support from studies P,[lO| of "complete" gauge fixing; 
i.e., in the outcome of attempts to construct a Fadde'ev- Popov-like determinant that eliminates Gribov 
copies or ensures that the functional integration domain for the gauge field is restricted to a subspace 
without them. Fixing a so-called "minimal Landau gauge," which enforces a constraint of integrating 
only over gauge field configurations inside the Gribov horizon; i.e., on the simplest domain for which 
the Faddeev-Popov operator is invertible, the dressed-gluon 2-point function is shown to vanish at 
k"^ = 0. However, the approach advocated in Refs. [^,0 makes no use of the additional ghost-like 
fields necessary to restrict the integration domain. 

Thus far in this discussion of the gluon DSE we have reported nothing qualitatively new and a more 
detailed review of the studies described can be found in Ref. [|], Sec. 5.1. What about contemporary 
studies? 



The direct approach to solving the Landau gauge gluon DSE, pioneered in Refs. p3| , |6^j65| , |66| , |67| , has 
been revived by two groups: A, Refs. |]58|j70| , [7T| ; and B, Refs. [[f2| , [f3t[f^ , [f5|] , with the significant new 
feature that nonperturbative effects in the ghost sector are admitted; i.e., a nonperturbative solution 
of the DSE for the ghost propagator is sought in the form 



G^\k) = -S' 



ab 



w{k^ 



[without interactions, wik"^) = 1]. 



(2.2.60) 



These studies analyse a truncated gluon-ghost DSE system, retaining only the third and sixth loop 
diagrams in the first equation of Fig. ( p^.2| ) and also the second equation. Superficially this is the same 
complex of equations as studied in Refs. ||68|j69[| . However, the procedure for solving it is different, 
arguably less systematic but also less restrictive. The difference between the groups is that A employ 
Ansdtze for the dressed-ghost-gluon and dressed-three-gluon vertices constructed so as to satisfy the 
relevant Slavnov- Taylor identities while B simply use the bare, perturbative vertices. Nevertheless they 
agree in the conclusion that in this case the infrared behaviour of the gluon DSE's solution is determined 
by the ghost loop alone: it overwhelms the gluon vacuum polarisation contribution. That is emphasised 
in Ref. |]73|, which eliminates every loop diagram in truncating the first equation of Fig. |2.2| except the 



ghost loop and still recovers the behaviour of Ref. . That behaviour is 



w(A;2 



{k 



1 

2Yk ' 



dik"^) ~ {k 



2\2k 



for e < a: 



QCD ■> 



with 0.8 <K<1. 



(2.2.61) 
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Exact evaluation of the angular integrals that arise when solving the integral equations gives the integer 
valued upper bound 
at fc^ 



K 



1 73 



This corresponds to a dressed-gluon 2-point function that vanishes 
0, ahhough the suppression is very sudden with the propagator not peaking until k"^ ~ 



where 



A2 



{d{e)/e 



t2 — A2 
~ QCD 



100/A 
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QCD 1 



(2.2.62) 



i.e., it is very much enhanced over the free propagator. [See, e.g., Ref. Fig. 12.] k = 1 also yields 
a dressed-ghost propagator that exhibits a dipole enhancement analogous to that of Eq. ( p.2.58| ). A 
[renormalisation group invariant] strong running coupling consistent with this truncations is: 

and its value at /c^ = is fixed by the numerical solutions: 



(2.2.63) 





A 


B 


a{k^ = 0) 


9.5 


~ 4 or 12 



(2.2.64) 

[NB. Group A approximates the angular integrals and uses vertex Ansdtze. Group B uses bare vertices 
and in Ref. |j72| approximates the angular integrals to obtain a{0) ~ 12, while in Ref. |^3[ the integrals 
are evaluated exactly, which yields a{0) = |7r ^ 4.2.] 

The qualitative common feature is that the Grassmannian ghost loops act to suppress the dressed-gluon 
propagator in the infrared. That may also be said of Refs. |^8|,^. [Indications that the quark loop, 
diagram seven in Fig. |2.2| , acts to oppose an enhancement of the type in Eq. ( |2.2.58| ) may here, with 
hindsight, be viewed as suggestive.] One aspect of ghost fields is that they enter because of gauge fixing 
via the Fadde'ev-Popov determinant. Hence, while none of the groups introduce the additional Fadde'ev- 
Popov contributions advocated in Refs. P,p!0|, they nevertheless do admit ghost contributions, and in 



their solution the number of ghost fields does not have a qualitative impact. Reference also obtains 
a dressed-propagator for the Fadde'ev-Popov fields with a fc^ = dipole singularity. It contributes to 
the action via the term employed to restrict the gauge field integration domain, in which capacity the 
dipole singularity can plausibly drive an area law for Wilson loops. 

Schwinger functions are the primary object of study in numerical simulations of lattice-QCD and 
Refs. [0 report contemporary estimates of the lattice Landau gauge dressed-gluon 2-point function. 
They are consistent with a finite although not necessarily vanishing value of (i(fc^ = 0). However, sim- 
ulations of the dressed-ghost 2-point function find no evidence of a dipole singularity, with the ghost 
propagator behaving as if w{k'^) = 1 in the smallest momentum bins [^. [NB. Since the quantitative 
results from groups A and B differ and exhibit marked sensitivity to details of the numerical analysis, 
any agreement between the DSE results for w{k'^) or d{k'^) and the lattice data on some subdomain can 
be regarded as fortuitous.] 



The behaviour in Eqs. (|2.2.61| ) also entails the presence of particle-like singularities in extant Ansdtze for 
the dressed-ghost-gluon, dressed-three-gluon and dressed-quark-gluon vertices that are consistent with 
the relevant Slavnov- Taylor identities, [k = 1 corresponds to an ideal simple pole singularity.] Hence 
while this behaviour may be consistent with the confinement of elementary excitations, as currently 
elucidated it also predicts the existence of coloured bound states in the strong interaction spectrum. 
Furthermore, while it does yield a strong running coupling with a{k'^ = 0) > 1, that makes DCSB 
dependent on fine tuning [^, and quantitative calculations based on the present numerical solutions 
give a quark condensate only ~ 5% of the value in Eq. ( |2.2.52| ) [ ]75| . Notwithstanding these remarks, 

^,^j7^J7^ are laudable. They have focused 



the studies of Refs. [p8|[70t[71| and subsequently Refs. [^ 
attention on a previously unsuspected qualitative sensitivity to truncations in the gauge sector 
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To recapitulate. It is clear from Refs. [|^,0,|^,^J5^ that DCSB requires the effective interaction in the 



quark DSE to be strongly enhanced at k"^ ~ ^qcd- [Remember too that modern lattice simulations ^5 



confirm the pattern of behaviour exhibited by quark DSE solutions obtained with such an enhanced 
interaction.] Studies of QCD's gauge sector indicate that gluon-gluon and/or gluon-ghost dynamics 
can generate such an enhancement. However, the qualitative nature of the mechanism and its strength 
remains unclear: is it the gluon vacuum polarisation or that of the ghost that is the driving force? It 
is a contemporary challenge to explore and understand this. 

Finally, in discussing aspects of the gauge sector one might consider whether instanton configurations 
play a role? Instantons are solutions of the classical equation of motion for the Euclidean gauge 
field. As such they form a set of measure zero in the gauge field integration space. Nonetheless, 
they can form the basis for a semi-classical approximation to the gauge field action and models based 



on this notion have been phenomenologically successful ||T9[. In this context we note that the DSEs in 
Fig. ( p.2|) are derived nonperturbatively and their self-consistent solution includes the effects of all field 
configurations. Hence instanton-like configurations may contribute to the form of the solution. However, 
a successful description of observable phenomena does not require that their contribution be quantified 
in a particular truncation. Nevertheless, one might estimate that a dilute liquid of instantons, each 
with radius p ~ 1/(0.6 GeV), could significantly effect the propagation characteristics of gluons only 
on the domain of intermediate momenta: fc^ ~ (0.6 GeV)^. Therefore they cannot qualitatively affect 
the infrared aspects discussed in this subsection. They also make no contribution in the perturbative 
domain: fc^ > 1-2 GeV^, where the perturbative matching inherent in the DSEs is a strength that 
makes possible a unification of infrared and ultraviolet phenomena, such as in the behaviour of bound 



state elastic and transition form factors; e.g., Refs. ||80|,|81|,|82|,83,84,85 



Confinement Confinement is the failure to directly observe coloured excitations in a detector: nei- 
ther quarks nor gluons nor coloured composites. The contemporary hypothesis is stronger; i.e., coloured 
excitations cannot propagate to a detector. To ensure this it is sufficient that coloured n-point func- 
tions violate the axiom of reflection positivity [^], which is guaranteed if the Fourier transform of the 
momentum-space ra-point Schwinger function is not a positive-definite function of its arguments. Re- 
flection positivity is one of a set of five axioms that must be satisfied if the given n-point function is to 
have a continuation to Minkowski space and hence an association with a physical, observable state. If 
an Hamiltonian exists for the theory but a given n-point function violates reflection positivity then the 
space of observable states, which is spanned by the eigenstates of the Hamiltonian, does not contain 
anything corresponding to the excitation(s) described by that Schwinger function. [The violation of 
reflection positivity is not a necessary condition for confinement 
massless two-dimensional QED 



A text-book counterexample is 
but in this case confinement of electric charge and DCSB both arise 



as a peculiar consequence of the number of dimensions.] 

The free boson propagator does not violate reflection positivity: 



A(x) :-- 



(2^ 



Ak-x 



^2 _|_ ^2 Jo 



diJiiix) 



m 



An'^x 



Ki{mx) . 



(2.2.65) 



Here x := (x-xY^'^ > 0, Ji is an oscillatory Bessel function of the first kind and Ki is the monotonically 
decreasing, strictly convex-up, non-negative modified Bessel function of the second kind. The same is 
true of the free fermion propagator: 



S{x) 



d^k 



ik-x 



m — ij ■ k 

^2 _|_ ^2 



(m — 7 ■ d)A{x) 



m 



Att'^x 



Kiimx) H K2(mx] 
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which is also positive definite. The spatially averaged Schwinger function is a particularly insightful 
tool |53J8^. Consider the fermion and let T = X4 represent Euclidean "time," then 



/I 
d'xtiDSix^T) = - / 
TT JO 



di- 



m 



cos(£T) 



1 -mT 

2 ^ 



^2.2.67) 



Hence the free fermion's mass can be easily obtained from the large T behaviour of the spatial average: 

mT = - lim \nas{T) . (2.2.68) 

T — ^00 

[The boson analogy is obvious.] This is just the approach used to determine bound state masses in 
simulations of lattice-QCD. 

For contrast, consider the dressed-gluon 2-point function in Eq. ( |2.2.59| ): 



D{x) :-- 



Ak-x 



A;4 + y 4vr2x 



diJiiix) 



7 



d 



£4 + 



Att'^x \ dz 



ker(2;) 



(2.2.69) 



2=7X 



where ker(2;) is the oscillatory Thomson function. D{x) is not positive definite and hence a dressed- 
gluon 2-point function that vanishes at A;^ = violates the axiom of refiection positivity and is therefore 
not observable; i.e., the excitation it describes is confined. At asymptotically large Euclidean distances 



D(x) 



OC 



1/2 



X 



3/2 



cos(^7x + f ) + sin(^7x + f 



%/2 



^2.2.70) 



Comparing this with Eq. ( 2.2.65| ) one identifies a mass as the coefficient in the exponential: mo = 7/v^2. 
[NB. At large x, Ki{x) oc exp{—x)/y/x.] By an obvious analogy, the coefficient in the oscillatory term 
is the lifetime ||68|j69|l : r = l/m£,. Both the mass and lifetime are tied to the dynamically generated 
mass-scale 7, which, using 

z 1 1 

' ' (2.2.71) 



z'^ + 'j^ z + i'j^ ~ z — I'j^ 

is just the displacement of the complex conjugate poles from the real-fc^ axis. It is a general result that 
the Fourier transform of a real function with complex conjugate poles is not positive definite. Hence 
the existence of such poles in a n-point Schwinger function is a sufficient condition for the violation of 
refiection positivity and thus for confinement. The spatially averaged Schwinger function is also useful 
here. 



D{T) := / d'xD{x,T) = - 



TC Jo 



i^ + Y 



1 
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cos(£T) = ^e'^"'^ cos(^7r + f ) 



and, generalising Eq. (|2.2.68|) , one can define a T-dependent mass: 

m(T) T := - \nD{T) = ln(27) + 1 7 T - In [cos(^7T + f 



^2.2.72) 



;2.2.73) 



It exhibits periodic singularities whose frequency is proportional to the dynamical mass-scale that is 
responsible for the violation of refiection positivity. If a dressed-fermion 2-point function has complex 
conjugate poles it too will be characterised by a T-dependent mass that exhibits such behaviour. 

This refiection positivity criterion has been employed to very good effect in three dimensional QED . 
First, some background. QED3 is confining in the quenched truncation That is evident in the 

classical potential 



V{r) :-- 



dx-^ 



d^k 
(2^ 



^ik-x+ik^X'j, 



P 27r 



In(e^r) , 



2 I 2 

•^1 "^2 ' 



(2.2.74) 
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and one finds 91 



which describes the interaction between two static sources. [NB. has the dimensions of mass in 
QED3.] It is a logarithmically growing potential, showing that the energy required to separate two 
charges is infinite. Furthermore, V{r) is just a one-dimensional average of the spatial gauge-boson 2- 
point Schwinger function and it is not positive definite, which indicates that the photon is also confined. 

If now, however, the photon vacuum polarisation tensor is evaluated at order using Nf massless 
fermions then, using the notation of Eq. ( |2.2.15| ), the photon propagator is characterised by 

^ = , from Uik-) = ^,a = NfC^S , (2.2.75) 

Vir) = ^ [Ho(ar) - iVo(ar)] , (2.2.76) 

where Ho(x) is a Struve function and No{x) a Neumann function, both of which are related to Bessel 
functions. In this case V{r) is positive definite, with the limiting cases 

r(r) ~°-ln(«r), \/(r) =°° |^ , (2.2.77) 

and confinement is lost in QED3. That is easy to understand: pairs of massless fermions cost no energy 
to produce and can propagate to infinity so they are very effective at screening the interaction. 

With (i(A;^) = 1/[1 + n(/i;^)] and sensible, physical constraints on the form of n(A;^), such as boundedness 
and vanishing in the ultraviolet, one can show that |^ 



1 

V(r) ™ — In(eV) + const. + h(r) , (2.2.78) 

^ ' 27r 1 + n(0) ^ ' ^ ^ ' 

where h{r) falls-off at least as quickly as 1/r. Hence, the existence of a confining potential in QED3 just 
depends on the value of the vacuum polarisation at the origin. In the quenched truncation, 11(0) = 
and the theory is logarithmically confining. With massless fermions, 1/[1 -|- 11(0)] = and confinement 
is absent. Finally, when the vacuum polarisation is evaluated from a loop of massive fermions, whether 
that mass is obtained dynamically via the gap equation or simply introduced as an external parameter, 
one obtains 11(0) < 00 and hence a confining theory. 



In Ref. the QED3 gap equation is solved for all four cases and the fermion propagator analysed. The 
results are summarised by Fig. |2.3| . In the quenched theory, Eq. ( ^.2.74[ ), the dressed- fermion 2-point 
function exhibits exactly those periodic singularities that, via Eq. ( |2.2.73| ), are indicative of complex 
conjugate poles. Hence this feature of the 2-point function, tied to the violation of reflection positivity, 
is a clear signal of confinement in the theory. That is emphasised further by a comparison with the 
theory that is unquenched via massless fermions in the vacuum polarisation, Eq. ( p.2.75|) . As we have 



described, that theory is not confining and in this case (Js{T) has the noninteracting, unconfined free 
particle form in Eq. ( p.2.67|) . The difference could not be more stark. The remaining two cases exhibit 



the periodic singularities that signal confinement, just as they should based on Eq. ( [^.2. 78] ). 



At this point we note that any concern that the presence of complex conjugate singularities in coloured 
rz-point functions leads to a violation of causality is misguided. Microscopic causality only constrains 
the commutativity of operators, and products thereof, that represent elements in the space of observable 
particle states; i.e., the space spanned by eigenstates of the Hamiltonian. Since Schwinger functions 
that violate reflection positivity do not have a continuation into that space there can be no question 
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Fig. 2.3. A(T) := -as{T) from Eq. ( p.2.67|) for QED3 with 2 flavours of fermion. 
O: confining, quenched theory; +, massless fermions used to evaluate the photon vacuum 
polarisation tensor; x, as before but with fixed-mass fermions; □, again but fermions with 
a momentum-dependent mass function. (Adapted from Ref. 



of violating causality. It is only required that iS-matrix elements that describe colour- singlet to colour- 
singlet transitions should satisfy the axioms, including reflection positivity. 

The violation of reflection positivity by coloured n-point functions is a sufficient condition for con- 
finement. However, it is not necessary, as the example of planar, two-dimensional QCD shows [Q. 
There the fermion two-point function exhibits particle-like singularities but the colour singlet meson 
bound state amplitudes, obtained from a Bethe-Salpeter equation, vanish at momenta coincident with 
the constituent-fermion mass shell. This excludes the pinch singularities that would otherwise lead to 
bound state break-up and liberation of the constituents. It is a realisation of confinement via a failure 
of the cluster decomposition property [CDP] PJ9^. The CDP is a requirement that the difference 
between the vacuum expectation value of a product of fields and all products of vacuum expectation 
values of subsets of these fields must vanish faster than any power. [This is modified slightly in theories, 
like QED, with a massless, asymptotic state: the photon.] It can be understood as a statement about 
charge screening and its failure means that, irrespective of the separation between sources, the inter- 
action between them is never negligible. That is an appealing, intuitive representation of confinement. 
Failure of the CDP is an implicit basis for confinement in the bulk of QCD potential models with; e.g., 
Refs. [p5| , p6| providing contemporary illustrations. 



Confinement is a more contentious issue than DCSB and its origin and realisation less-well under- 
stood. However, in this subsection we have described a perspective that is common to many authors, 
and an interested reader will find variations and more expansive discussions of various points in; e.g., 
Refs. i0j2|j63j6|j6i, and complementary perspectives in; e.g., Refs. p7|,|98|,p9i[l00 i|l0l| , |102|J103^ . It is, 
of course, because confinement is poorly understood that its study and modeUing are important. To 
presume otherwise is a misapprehension. 
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2.3 Phenomenological Applications 

When Ref. ^ was written the only apphcation of DSEs to observable phenomena consisted in the 
oft-repeated calculation of well-known quantities, such as the pion mass and decay constant. Refer- 
ences [^ , |104|| represent a departure from that, and are progenitors of the wide-ranging application 



of DSEs to observables newly accessible at the current generation of experimental facilities. Many of 
these applications are reviewed in Refs. [pl] , |105|JT06| and herein we only describe three recent, significant 
developments. 



Light Mesons The model used to illustrate renormalisation and DCSB in Sec. 2T2 has been applied 



to the calculation of vector meson masses and decay constants |52], and to elucidate the role of vector 



mesons in connection with the electromagnetic pion form factor [KT?]. These studies are important 



because, in concert with Ref. pO|], they complete a DSE description of those light-mesons in the strong 



interaction spectrum that are most often produced in reactions involving hadrons. In so doing they 
illustrate the efficacy of the renormalisation-group-improved rainbow-ladder truncation for flavour non- 
singlet pseudoscalar and vector mesons composed of light-quarks {u, d, s), and thereby that of the 



systematic, Ward-Takahashi identity preserving truncation scheme introduced in Ref. |108 



The renormalised homogeneous Bethe-Salpeter equation for a bound state of a dressed-quark and 
dressed-antiquark with total momentum P is 

[Tnik; P)]tu = r P)]rsKi:{q, k; P) , P) ■■= S{q+) Tniq; P) S{q+) , (2.3.1) 

Jq 

with: VH{k]P) the Bethe-Salpeter amplitude [BSA], where H specifies the fiavour structure of the 
meson; S{p) := diag[5'„(p), Sdip), Ss{p)]; q+ = q + iJpP, g_ = g - (1 - ■r]p)P; and r,. . .,u represent 
colour-, Dirac- and fiavour-matrix indices, [rjp G [0, 1] is the momentum partitioning parameter. It 
appears in Poincare covariant treatments because, in general, the definition of the relative momentum 
is arbitrary. Physical observables, such as the mass, must be independent of rjp but that is only possible 
if the Bethe-Salpeter amplitude depends on it. rjp = 1/2 for charge-conjugation eigenstates.] 



In Eq. ( p^.3.1| ), Kl^{q,k;P) is the renormalised, fully-amputated quark-antiquark scattering kernel. 



which, as we have seen, also appears implicitly in Eq. (|2.2.1|) because it is the kernel of the inho- 
mogeneous DSE satisfied by T^{q;p). K^^^q, k; P) is a 4-point Schwinger function, obtained as the 
sum of a countable infinity of skeleton diagrams. It is two-particle- irreducible, with respect to the 
quark-antiquark pair of lines and does not contain quark-antiquark to single gauge-boson annihilation 
diagrams, such as would describe the leptonic decay of a pseudoscalar meson. [A connection between the 
fully-amputated quark-antiquark scattering amplitude: M = K + K{SS)K + . . . , and the Wilson loop 
is discussed in Ref. ||100|| .] The complexity of Kl^{q, k; P) is one reason why quantitative studies of the 
quark DSE currently employ Ansdtze for D^y{k) and Ty{k,p). However, as illustrated by Ref. [|^, the 
complexity of Kl^{q, k; P) does not prevent one from analysing aspects of QCD in a model independent 
manner and proving general results that provide useful constraints on model studies. 



Equation (|2.3.1| ) is an eigenvalue problem and solutions exist only for particular, separated values of 



P^. The eigenvector associated with each eigenvalue: TH{k;P), the BSA, is a one-particle-irreducible, 
fully-amputated quark-meson vertex. In the fiavour non-singlet pseudoscalar channels the solutions 
having the lowest eigenvalues correspond to the vr- and if-mesons, while in the vector channels they 
correspond to the u-, p- and (p mesons. 



Following Ref. ||108|| , the renormalised inhomogeneous BSE for the axial- vector vertex, consistent with 
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the renormalisation-group- improved quark DSE, Eqs. ( |2.2.39| ) and ( |2.2.40| ), is 



T[^{k- P) = Z, iA'757M - f^iik - qf) D^f^k - q)^l^S{q^) T[^{q- P) S{q^) ^7/3 , (2.3.2) 



where {\\f '■ ^ = • • • )8} are the generators of 5'f/(3)fiavour. It is straightforward to verify that the 
axial- vector Ward-Takahashi identity is satisfied; i.e., 



P.rUk; P) = S-\k+) + iA'^75 S~\k_) - ^r^(fc; P) - zT^ik; P) M, 



(2.3.3) 



where = diag[m„((^), md((^), ms((^)] and the renormahsed pseudoscalar vertex satisfies its own inho- 
mogeneous BSE: 



T\{k- P) = Z, iA'75 - I Q{{k - qf) D^^:%k - g)-7, 5(g+) T\{q- P) S{qJ, 
[NB. The product M(^T''^{k; P) is renormahsation point independent.] 



(2.3.4) 



The pseudoscalar mesons appear as poles in both the axial-vector and pseudoscalar vertices and 
equating pole residues yields the homogeneous BSE 



Tnik; P) + f g{{k - qf) D^^^k - q)^^, 5(g+) r^(g; P) 5(g_) y 7. = . 



(2.3.5) 



As is characteristic of homogeneous equations, the normalisation of the solution is not fixed by this 
equation. The canonical normalisation enforces a requirement that the bound state contribution to the 
fully-amputated quark- ant iquark scattering amplitude: M, have unit residue. In this rainbow-ladder 
truncation that condition is expressed via 



2P„ 



tr 



^Hiq; -P) Tu{q; P) S{q^) + f ^(g; -P) S{q^) T u{q; P) 



dPa 



dPa 



where 



TH{k,~PY ■.= C-^TH{-k,-P)C , 
with C = 7274 the charge conjugation matrix: 

Ci'.c^ = -7m ; [c, 75] = , 

and X* denotes the matrix transpose of X. The general form of a pseudoscalar BSA is 



r^(A;; P) = T% 



lEnik; P) + 7 ■ PFnik; P)+-fkk- PGnik; P) + a^, k^P, Hnik; P) 



(2.3.6) 
(2.3.7) 
(2.3.8) 

(2.3.9) 



where is a matrix that describes the flavour content of the meson; e.g., T'^^ = ^{X];- + iX'jp) and, 
for bound states of constituents with equal current-quark masses, the scalar functions E, F, G and 
H are even under k ■ P ^ —k ■ P. [NB. Since the homogeneous BSE is an eigenvalue problem, 
Efjik; P) = Enik"^, k ■ P|P^); i.e., P^ is not a variable, instead it labels the solution. The same is true 
of each function.! 



Equation ( |2.3.5| ) also describes vector mesons, as can be shown by considering the inhomogeneous 
equation for the renormahsed vector vertex. In general, twelve independent scalar functions are required 
to express the Dirac structure of a vector vertex. However, a vector meson bound state is transverse: 



P,T^{k\k.P\P' 



-m 



H) 







(2.3.10) 
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where tuh is the bound state's mass, and this constraint reduces to eight the number of independent 
scalar functions. One therefore has 

r;^{k;P) = T''j2 0lFi{k-P), (2.3.11) 



1=1 



with eight ort honor mahsed matrix covariants 



01 



01 
01 
01 



T 

2fcT^ . p , 



[1 + (A; • Pf 



T 

7^! 7 



P 



01 
01 



V2 (7; 



T ^, _ £,T 



where: 7 J : = 



k 



T ._ 



^^.F7-P + fcj7-i') , 

72 (7^7'^ -A;- 7^^ -kil) , 

«^fcj7^-^7-^, 

(2.3.12) 

+ k ■ P Pn- With this decomposition 



2i 



7^ + 7 ■ PP^, P ■ P = -1; and ■ = 1 
the magnitudes of the invariant functions, P;, are a direct measure of the relative importance of a given 
Dirac covariant in the BSA; e.g., one expects Pi to be the function with the greatest magnitude for 
J^*-^ = 1 bound states. 



To calculate the meson masses, one first solves Eqs. ( |2.2.39| ) and ( [^.2. 401 ) for the renormalised dressed- 
quark propagator. This numerical solution for S{j)) is used in the pseudoscalar BSE, Eq. (|2.3.5| ) under 
the substitution of Eq. ( p.3.9|) , which is a coupled set of four homogeneous equations, one set for each 
meson; and the vector BSE, Eq. ( |2.3.5| ) with Eq. ( p. 3.111) , a coupled set of eight equations for each 
meson. Solving the equations is a challenging numerical exercise, requiring careful attention to detail, 
and two complementary methods were both used in Refs. pO|j5^ . While the numerical methods were 
identical, the authors of Ref. [^] used a simplified version of the effective interaction: 



fc2 



-Dk\-'''- 



4tt ■ 



iln 



l + P/A^ 



QCD 



■J^ik' 



(2.3.13) 



and varied the single parameter D along with the current-quark masses: m„ = rhd and rhg, in order 
to reproduce the observed values of m,r, f^K and /,r- All other calculated results are predictions in the 
sense that they are unconstrained. (NB. The Poincare invariant four- dimensional BSE is solved directly, 
eschewing the commonly used artefice of a three-dimensional reduction, which introduces spurious effects 
when imposing compatibility with Goldstone's theorem and also leads to a misinterpretation of a model's 
parameters ||109||.) 



Before reporting the results it is necessary to introduce the formulae for the meson decay constants: 
fn, which completely describe the strong interaction contribution to a meson's weak or electromagnetic 
decay. Following Ref. the pseudoscalar meson decay constant is given by 
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/hP, := (0| Qir^)%j,Q\H{P)) = Z,ti {r'')\,^^Sik^)TH{k;P)S{k^), (2.3.14) 

where here Q = column(M, rf, s). The factor of Z2 on the r.h.s. ensures that fn is gauge invariant, 
and independent of the renormalisation point and regularisation mass-scale; i.e., that Jh is truly an 
observable. Equation ( f^.3.14D is the pseudovector projection of the unamputated Bethe-Salpeter wave 
function, x(/c;P), calculated at the origin in configuration space. As such, it is one field theoretical 
generalisation of the "wave function at the origin," which describes the decay of bound states in quantum 
mechanics. The analogous expression for vector mesons is 



^/^,mHeJ(P) 



H\t 



1,Q\H{P)) 

A 



^JnmH = \Z,iT {T'')\,S{k^)T^{k;P)S{k^), (2.3.15) 
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Table 2.1. Masses and decay constants [in GeV] of light vector and flavour nonsinglet pseu- 
doscalar mesons calculated using the renormalisation-group-improved rainbow-ladder trun- 
cation. The underlined quantities were fltted. The "Obs." value of the vacuum quark 
condensate is the global estimate of Ref. and the masses are taken from Ref. |0 



as 



are /jr, fx- The vector meson decay constants are discussed in connection with Eqs. (|2.3.18 
2.3.201) . For = m^, the rainbow-ladder truncation gives m^^ = mp. This degeneracy 
is lifted by meson-loop self-energy contributions, such as p ^ titi —>■ p [^ , |llCI| , pTT| . The 



root-mean-square error over predicted quantities is just 3.6 %. (Adapted from Ref. [p^ 









rriK 


nip 






U 


fK 


fp 


fK* 


f(t> 


Obs. 


0.236 


0.139 


0.496 


0.770 


0.892 


1.020 


0.130 


0.160 


0.216 


0.225 


0.238 


Gale. 


0.242 


0.139 


0.496 


0.747 


0.956 


1.088 


0.130 


0.154 


0.197 


0.246 


0.255 



where e^(-P) is the vector meson's polarisation vector: P ■ e^{P) = 0. 



A best-flt is obtained |52] with 

D = (1.12GeV)^ (2.3.16) 

which is a 60% increase over Eq. (|2.2.36| ), as expected because the single Gaussian term in Eq. ( |2.3.13| ) 
must here replace the sum of the flrst two terms in Eq. (|2.2.35|) ; and renormalised current-quark masses 



m. 



.^G^v = 5.5 MeV, 



m 



iGeV ^ 124 MeV 



(2.3.17) 



which are little changed from the values used in Ref. ||2^, Eq. ( |2.2.55| ). [NB. uj and rrit are unchanged. 
See the discussion preceding Eq. ( p.2.36|) .] The results are presented in Table and are characterised 
by a root-mean-square error over predicted quantities of just 3.6 %. We emphasise that this is obtained 
with a one-parameter model of the effective interaction. 

Experimental values of the pseudoscalar meson decay constants are obtained directly via observation of 
their prominent (3-decaj mode. However, for the vector mesons this mode is not easily accessible and 
to proceed we note that the rainbow-ladder truncation predicts ideal flavour mixing. Using this and 
isospin symmetry, one can relate the fp matrix element to that describing p — > e~^e~ decay: 



m 
9p 



with Qe 
0.32 keV |3l 



diag[2/3, — 1/3, — 1/3]; i.e., the quark's electromagnetic charge matrix. FpO^e+e 
Qp = 5.03 ± 0.12 and hence fp = 216 ± 5 MeV. For the 0-meson 



m 







9^ 



(2.3.18) 
= 6.77 ± 

(2.3.19) 



and hence F^ 



1.37 ± 0.05 keV ^ 

^T—>K*Ut T—^pUr 



= 12.9 ±0.2 or = 238 ±4 MeV. fx* follows from |2l 
0.051 ^ /x* = 1.042 /p. (2.3.20) 



A number of other important observations are recorded in Refs. [^,^. First: The calculated values 
of observable quantities are independent of the momentum partitioning parameter: r^p, when all of the 
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Dirac covariants, and their complete momentum dependence, are retained; i.e., Poincare invariance is 
manifest. Second: For pseudoscalar mesons the leading 75-covariant is dominant but the pseudovector 
components also play an important role; e.g., fx is < 30% smaller without them. For the vector mesons, 
while Fi is dominant, -F2...5 are also important; e.g., nip is > 20% larger without them. Third: Flavour 
nonsinglet pseudoscalar mesons obey a mass formula P0| , P2|PT2| , exact in QCD, 



fnmH 



H 5 



H 



tr^[Mar^,(r^)n] 



where 



H 



-2i 



c 



1 

7^ 



(2.3.21^ 



(2.3.22) 



is the gauge-invariant and cutoff-independent residue of the pion pole in the pseudoscalar vertex. As a 
residue, it is an analogue of fn and describes the pseudoscalar projection of the unamputated Bethe- 
Salpeter wave function calculated at the origin in configuration space. For small current-quark masses. 



Eq. ( 2.3.21 ) yields the "Gell-Mann-Oakes-Renner" relation as a corollary. However, it is also valid 
for heavy-quarks and predicts ||38| , |113|| itlh oc mq in the heavy-meson domain, which is verified in the 
strong interaction spectrum. Fourth: The behaviour at large fc^ [ultraviolet relative momenta] is model 
independent, determined as it is by the one- loop improved strong running coupling. Fifth: For the 
pseudoscalar mesons, the asymptotic behaviour of the subdominant pseudovector amplitudes [Fh, Gh] 
is crucial for convergence of the integral describing fn- For the vector mesons, the same is true of -F2...5- 
This is also a model-independent result because of the fourth point. 

This subsection illustrates the reliability of the rainbow-ladder truncation for light vector and flavour 
nonsinglet pseudoscalar mesons 
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It 



That is not an accident but rather, as elucidated in Ref 
is the result of cancellations between vertex corrections and crossed-box contributions at each higher 
order in the quark-antiquark scattering kernel. There are two other classes of light-meson: scalar and 
axial-vector 



A separable model |114 



that expresses characteristics of the rainbow-ladder truncation 
has been employed successfully in calculating the masses and decay constants of M,(i-quark axial- vector 
mesons [p.l5|| . Hence, while a more sophisticated study is still lacking, and is indeed required, this 
suggests that the truncation can provide a good approximation in this sector; i.e., that the conspiratorial 
cancellations are also effective here. 

In the scalar channel, however, the rainbow-ladder truncation is not certain to provide a reliable ap- 

This is entangled with the 



proximation because the cancellations described above do not occur ||11( 
phenomenological difficulties encountered in understanding the composition of scalar resonances below 
1.4 GeV [p^ , |117|JTT8| , |119|| . For the isoscalar-scalar meson the problem is exacerbated by the presence of 
timelike gluon exchange contributions to the kernel, which are the analogue of those diagrams expected 
to generate the rj-r]' mass splitting in BSE studies | |120| ]. If the rainbow-ladder truncation is employed 



one obtains ideal flavour mixing and degenerate isospin partners, and; e.g.. 



Ref. 


{u/d)i=Q^i us ss 


calculated mass in cev [47] 

m 

114] 
121] 


0.59 0.90 1.20 
0.67 

0.71 1.18 1.54 
0.59 


averaged mass 


0.64 ±0.06 1.04 ±0.20 1.37 ±0.24 



(2.3.23) 



( 0"*"^ mesons containing at least one s-quark were not considered in Refs. [^|121|] .) Each model rep- 
resented in this compilation was constrained to accurately describe vr- and i^'-meson observables, and 
the standard deviation about their averaged vector meson masses is a factor of 2-5 smaller than that 
exhibited in the last row here. 
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In Eq. ( p .3.231 ) we do not compare directly with observed masses because of the uncertainty in identifying 
the members of the scalar nonet. We only note that: 1) the us channel is least affected by those 
corrections to the rainbow-ladder truncation that can alter the feature of ideal flavour mixing and 
hence it may be appropriate to identify the us scalar with the Kq{1430), in which case the mass is 
underestimated by < 30 %; and 2) an analysis of tttt data | |118| identifies an isoscalar-scalar with 

0.46 GeV , ^0++^^^^ ~ 0.22-0.47 GeV . (2.3.24) 



This supports ideal flavour mixing but, because Ta-/m„ is large, suggests an accurate calculation of this 
state's mass will require a Bethe-Salpeter kernel that explicitly includes couplings to the important tttt 
mode. In contrast, that coupling can be handled perturbatively in the u-p sector P7| , |110| , pn|] . If this 
identification is correct then the mass estimate in Eq. ( p.3.23|) is < 40 % too large. 

The discussion here makes plain that the constituent-quark-like rainbow-ladder scalar bound states 
are significantly altered by corrections to the BSE's kernel. That is good because it is consistent 
with observation: understanding the scalar meson nonet is a complex problem. (NB. The difficulties 

ra- 



in Ref. 



Using the DSEs this complexity is expressed 



to be anticipated are illustrated; e.g., 
in unanswered questions. For example, why do timelike gluon exchange contributions to the kernel 
in the isoscalar-scalar channel not force a deviation from ideal flavour mixing and, returning to the 
pseudoscalar sector, what effect do they play in the rj-rj' mass splitting? 

There are other contemporary questions. For example, which improvements to the rainbow-ladder 
kernel are necessary in order to study bound states containing at least one heavy-quark? The extent 
to which the cancellations elucidated in Ref. ||108|| persist as the current-quark mass evolves to values 
larger than is not known. Even though the rainbow-ladder truncation can provide an acceptable 



estimate of heavy-meson masses; e.g., Ref. I^^l; improvements are necessary and potentials derived from 
dual-QCD models; e.g., Refs. |97|j99(| , have been applied more exhaustively ||123|| . (Heavy-heavy-mesons 
are also amenable to study via heavy-quark expansions ||124||). Other composites admitted by QCD can 



also be considered. The rainbow-ladder truncation has recently been applied to "exotic" light-mesons, 
predicting a J^'" = 1 ^ meson with a mass of ~ 1.4 - 1.5 GeV ||125|| . ["exotic" because such a J^*" value 
is unobtainable in the qq constituent quark model.] However, there is a dearth of DSE applications 
to the glueball spectrum, which has been explored using other methods; e.g., Refs. [ P^jl26|JT?7| . A 
gluon-sector analogue of the rainbow-ladder truncation is an obvious starting point and such studies 
must precede any exploration of hybrid quark-gluon states | |128| . 



That applications in all these areas are actively being pursued and contemplated is an indication of a 
healthy discipline. 



Electromagnetic Pion Form Factor and Vector Dominance We have seen that the renormalisation 
group improved rainbow-ladder truncation of the quark- ant iquark scattering kernel: K, provides a 
good understanding of colour singlet, mesonic spectral functions. That is a key success since describing 
these correlation functions is a core problem in QCD ||79|. However, it is only a single step and one must 
proceed from this foundation to the study of scattering observables, which delve deeper into hadron 
structure. 



The best such observables to study are elastic electromagnetic form factors, because the probe is well 
understood, and the simplest "target" for a theorist is the pion, as long as the theoretical framework 
accurately describes DCSB. The electromagnetic pion form factor is a much studied observable but here, 
to be concrete, we focus on the application of DSEs to this problem, and in that connection Ref. ||129[ is 



a pilot. As our exemplar we choose Ref. ||107|| because it is a direct application of the effective interaction 
described in the previous subsection, Eq. ( |2.3.13| ), and it is the most complete study to date. 
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In the isospin-symmetric limit the renormahsed impulse approximation to the pion's electromagnetic 
form factor is 



(2.3.25) 



2 tr / f ^(A;; -p2)5(A;++) iQeTVk++, T^{k - q/2; pi) S{k^ 

Jk 



k + api/2 + Pq/2 and p2 ■= Pi + q- Here, r^(A;;P) is the pion BSA, which has the form 
in Eq. ( |2.3.9D , and S{k) = diag[Su=d{k), Su=d{k)]. No renormahsation constants appear explicitly 
in Eq. (|2.3.25| ) because the renormahsed dressed-quark- photon vertex: F^, satisfies the vector Ward- 
Takahashi identity: 

{Pi-P2),triipi,p2) = S-\pi) - S-\p2) . (2.3.26) 
Importantly, this also ensures current conservation: 

{Pl -p2)f^A,,{pi,P2) =0, 

and, using Eq. ( p.3.6| ), the correct normalisation of the form factor: 

F{q' = 0) = l; 



(2.3.27) 



(2.3.28) 

i.e., combining the rainbow-ladder truncation of the quark-antiquark scattering kernel with the impulse 
approximation yields a minimal, consistent approximation ||130|| . [NB. As the vr*^ is a charge conjugation 
eigenstate, F^o(g^) = 0, Vg^. The impulse approximation yields this result.] 



The only quantity in Eq. ( 2.3.25| ) not already known is F^. It is the solution of an inhomogeneous BSE, 
which in rainbow-ladder truncation is 



Tlik; P) = ^2 7m - /' Q{{k - q?) D^f^k - g)^7a 5(g+) Vl{k- P) S{q 



-7/3 



(2.3.29) 



It is straightforward to verify that Eq. ( |2.3.26| ) is satisfied, and owing to this the general solution of 
Eq. (2.3.21j) involves only eight independent scalar functions and can be expressed as 



ri{k;P) = rl^{k;P) + j:olF;'{k-p), 

1=1 

where the matrices 0|j are given in Eq. ( |2.3.12| ) and |^| 

Ff (A:; P) = ^S^(A;2 , kl) 7^ + ik+ + k_)^ [1^7 ■ {k+ + k_) Aa(A;2 , kl) + Asikl, kl 
,2 ,2\ iri7/L2\ , i7^l2m a /l2 7^2\ _ ^ik+) — F{k'i) 



(2.3.30) 



^Fiki, ki) = i [F{ki) + F{ki)] , Ap{ki, ki 



kl-k^_ 



(2.3.31) 



(2.3.32) 



with F = A,B; i.e., the scalar functions in the dressed-quark propagator. T^^{k;P) saturates the 
vector Ward-Takahashi identity, and the remaining terms are transverse and hence do not contribute 
to the r.h.s. of Eq. ( |2.3.26|) . 



The importance of determining the dressed-quark-photon vertex from Eq. ( |2.3.29| ) was recognised in 
Ref. ||131j| , where a solution was obtained using the simple model of Q{k'^) introduced in Ref. [p^ . 
As is readily anticipated, the dressed-vertex exhibits isolated simple poles at timelike values of P^. 
Each pole corresponds to a 1" ^ bound state, P^ = — mass-squared, and its matrix- valued residue is 
proportional to the bound state amplitude. For in the neighbourhood of any one of these poles the 
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Q' [GeV'] 



Fig. 2.4. Pion form factor calculated directly from the interaction in Eq. ( |2.3.13| ): solid 
line. The data are: open circles, Ref. ||132|| ; squares, Ref. ||133|| ; filled circles, Ref. ||134|| . No 



parameters were varied to obtain the result. The other curves depict results from simplified 
calculations: dotted line, retaining only the dominant vector covariants in Eq. ( p.3.30D , F^ ^^] 
dashed line, retaining only T^'^ calculated with the numerical solution of the quark DSE; 
dash-dot line, TJ^ik^P) = 7^, which violates the Ward-Takahashi identity. (Adapted from 
Ref. 1107 



behaviour of the pion form factor is primarily determined by the manifestation of that bound state in 
the spectral density. Vector meson dominance, in any of its forms, is an Ansatz to be used for 
extrapolating outside of these neighbourhoods. A direct solution of Eq. ( |2.3.29| ) obviates the need for 
such an expedient and also any need to fabricate an interpretation of an off-shell bound state. 



Using the interaction of Eq. (|2.3.13| ), with its single parameter fixed as discussed in connection with 
Eq. ( p.3.16| ), and solving numerically for the renormalised dressed-quark propagator, pion BSA and 
dressed-quark-photon vertex, Ref. ||107|| obtains the pion form factor depicted in Fig. |2]^ with 



-6 



dy 



0.68 fm cf. 



„Obs 



0.663 ±0.006. 



(2.3.33) 



y=0 



The complete calculation describes the data for both spacelike and timelike momenta, plainly exhibiting 
the evolution to a simple pole corresponding to the p-meson. Here the p is described by a simple pole 
on the real-P^ axis because the rainbow-ladder truncation excludes the vrvr contribution in the kernel. 
However, that can be included perturbatively | 87| , 110 ,111] and estimates show it yields no-more-than 
a 15 % increase in ||135|| . [NB. As illustrated by Refs. ||136| , |137|J138|| , the extension to i^'-meson form 



factors is straightforward but with the interesting new feature that, unlike the neutral pion's elastic 
form factor, Fko{Q'^) 7^ 0, WQ"^ > 0, because the neutral kaons are not charge- conjugation eigenstates.] 



Proponents of vector meson dominance Ansdtze | |139[ | have historically claimed support in the accuracy 
of the estimate ~ 6/m^ and in this connection one can ask what alternative does the direct calculation 
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describe? To address this Ref. [[107|| observed that, for G [— m^, 0.2 GeV^], an interpolation of the 
solution of Eq. ( p.3.29|) is provided by 

^lik; Q) ^ rf (A;; Q) - ^ — ^ ^ F[{e\Q' = -ml) , (2.3.34) 

where ^ 

Ff{e\Q^ = -ml) :=- f dx F[\k'^ ,ikmpx\Q^ = -ml) (2.3.35) 

are the leading Chebyshev moments of the five dominant scalar functions in the p-meson BSA and, fol- 
lowing Ref. 1^ , 1/ Qpis the residue of the p-meson pole in the photon vacuum polarisation. Substituting 
Eq. d^Xg) into Eq. ([2X^ ) yields 

F.m ^ F^'^m - ^ Fp^M) -7^-, , (2.3.36) 

where QpTm -Fp7r7r(Q^) is the impulse approximation to the pTrvr vertex, FpT^T^{Q'^ = —ml) = 1. Implicit in 
Eq. ( p.3.34| ) is a clear but not unique definition of an off-shell p-qq correlation, which yields Eq. ( |2.3.36| ) 
as a calculable approximation to the form factor wherein 

(rPj:=±^FpUO). (2.3.37) 

mj, Qp 

This is simply the vector meson dominance result corrected for nonuniversality of the strong and electro- 
magnetic p-meson couplings; i.e., Qp ^ gpmr, and, via -Fp^7r(0), for the nonpointlike nature of the off-shell 
p-qq correlation. Fp^.^^^Q'^) was not calculated in Ref. |p.071] and therefore here we make an estimate. 
Typically -Fp^^(O) ~ 0.5 ||110|| and using experimental values for the observables [including ^fp^^ = 6.05] 
we find 

«)2~0.5r^, (2.3.38) 
which is a significant suppression with respect to the "naive" vector meson dominance estimate. 

This illustration demonstrates that the pion charge radius is indeed influenced by the p-pole's contribu- 
tion to the 1~ ~ spectral density. However, that is unsurprising: the bound state poles are a significant 
feature of the dressed-vertex. More important is a realisation that the separation in Eq. ( p.3.34|) is 
completely arbitrary and hence so is the fraction of the charge radius attributed to the "off-shell p- 
meson." One can shift any amount of strength between the two terms and yet maintain an accurate 
interpolation of the dressed-vertex. For example. 



+ 



+ ml Q'^ + ml Q"^ + ml 



1 _ Q-i^{l+Q^/ml) 



(2.3.39) 



is an apparently nugatory rearrangement. However, only the first term has a pole at + ml = 0, and 
substituting Eq. ( p.3.39|) into Eq. (|2.3.34| ) and repeating the analysis one arrives at 

(O^ = e--^^F (0) 0.2 r^. (2.3.40) 

ml Qp 

It is obvious now that can be made arbitrarily small while preserving r'^^^. [NB. Eq. ( |2.3.40| ) 
and the procedure leading to it are no more contrived than Eq. ( |2.3.37|) because Eq. ( p.3.39|) can be 
interpreted as "unfreezing" the vector meson BSA; i.e., of allowing the bound state correlation to be 
suppressed off-shell: Ff{k^\Q^ = -ml) exp(-c^[l + gVm^]) Ff{k^\Q^ = -ml).] 



32 



For large spacelike the calculation of Ft^{Q^) directly from Q{k'^) is a computational challenge and for 
practical reasons the study in Ref. ||107|| was restricted to < 1 GeV. The effective interaction doesn't 



appear explicitly in the impulse approximation, Eq. ( |2.3.25| ), only the dressed-quark propagator, pion 
Bethe-Salpeter amplitude and dressed-quark-photon vertex. Hence the numerical analysis is markedly 
simplified if algebraic approximations to these Schwinger functions are employed. 



This was the approach adopted in Ref. ||130 | 
cations; e.g., Refs. |38|,|8^ JTU5| , [T06i[TTTi|l40| , |l41 



which has since been used in a wide range of appli- 
It is efficacious because the vector and axial-vector 



Ward-Takahashi identities can be used to motivate Ansdtze for and F^ that are expressed solely in 
terms of S. In the present context, algebraic models for S and F^r have been developed [^] that encode 
the important qualitative aspects of the DSE and BSE solutions in Ref. These forms, along with 
F^ determined via Eq. ( p.3.31| ), make possible a calculation of F^(Q^), WQ^ > 0, and an analytic anal- 
ysis of the asymptotic behaviour. One finds that the pion's pseudoscalar covariant can alone provide 
a quantitative description of -F,r(Q^) for < 5GeV^. However, beyond this point the pseudovector 
covariants become important. It is these terms in Eq. ( |2.3.9|) that ensure 



Q^F^(Q^) = const.. 



(2.3.41^ 



up to calculable ln(Q/AQCD)^-corrections, and this behaviour is unmistakable for > lOGeV^. The 
anomalous dimension: p, is determined by that of F„ and Gt^, which is a model independent result. [NB. 
If the pseudovector covariants are neglected, -^^(Q^) = const.. Further, the dressed-quark propagator 
obtained in Ref. necessarily yields results inconsistent with Eq. ( |2.3.41| ).] This application is an 
emphatic demonstration of the DSEs' ability to provide a single and simultaneous description of the 
infrared and ultraviolet aspects of observables. 



That is also much in evidence in the study of anomalous processes. As first observed ||142|| in connection 
with processes like K^K~ — *■ 7r"'"7r°7r~, the systematically truncated DSEs yield the anomalies of current 
algebra without any model dependence; i.e., the anomalies remain a feature of the global aspects of 
DCSB [|143|| . For the 7r° — > 77 process this feature was verified in Refs. |]82| , |130| , p!4^ and for the 771* tttt 
transition form factor, in Ref. [ |145| |. Consequently the DSEs provide a single framework wherein the 
value of such transition form factors is fixed and model-independent at the soft-pion, zero momentum 
transfer point, and their evolution is calculable on the entire range of momentum transfer, reproducing 
the ultraviolet behaviour anticipated from perturbative QCD. These features have been elucidated and 
exploited in Refs. PlRPl. 



Describing Baryons Hitherto we have described the application of DSEs to the study of mesonic 
observables, which requires and illustrates a contemporary understanding of the two-body problem in 
quantum field theory. Baryons, as a three-body problem, pose a greater challenge and historically they 
have been described using constituent-quark-like models, which remain an important contemporary tool; 
e.g., Refs. pB|jl46| |. We identify a beginning of progress with a direct assault on the baryon problem in 
a realisation [p.47| , |148|| that field theoretical models of the strong interaction admit the construction of 



a meson-diquark auxiliary-field effective action and thereby a description of baryons as loosely-bound 
quark-diquark composites. 



A head-on DSE approach begins with a relativistic Fadde'ev equation, which can be derived ||149|| by 



exploiting the fact that single gluon exchange between two quarks is attractive in the colour-antitriplet 
channel, whether or not the gluon is dressed. Indeed, using a rainbow-ladder truncation of the quark- 
quark scattering kernel, one obtains nonpointlike, colour-antitriplet diquark [quark-quark] bound states 
in the strong interaction spectrum ||150|| . However, as demonstrated in Ref. [|108|| , this is a flaw of the 
rainbow-ladder truncation: higher order terms in the kernel ensure that the quark-quark scattering 
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p'^p+p' 



p.=^p-p 





P[=(l-ri )P+p 



Pb=(l-Tl)P-p' Pb=(l-ri)P+P 



Fig. 2.5. A quark-diquark Fadde'ev equation for a baryon with total momentum P: single 
line, quark propagator; double line, diquark propagator; quark-diquark Bethe-Salpeter 

amplitude. [For nonpointlike diquarks, the filled circle represents a diquark Bethe-Salpeter 
amplitude.] Here binding is effected by a participant quark leaving the diquark and joining 
with the dormant quark to form another diquark; i.e., via iterated quark exchange. (Adapted 



from Ref. |T54|. 



matrix does not exhibit the simple poles that correspond to asymptotic states. Nevertheless, studies 
with improved kernels [|151|] do support a physical interpretation of the spurious rainbow-ladder diquark 
masses. Denoting the mass by m^g, then Iqq := l/rriqq represents the range over which a true diquark 
correlation in this channel can persist inside a baryon. In this sense they are "pseudo-particle" masses 
and can be used to estimate which two-body correlations should be retained in solving the Fadde'ev 
equation. (NB. Gluon mediated interactions are repulsive in the colour-sextet quark-quark channel, 
just as they are in the colour-octet meson channel ||150||.) 



Reference ||114|] tabulates calculated values of these pseudo-particle masses, from which we extract: 





{ud)o+ 


{us)o+ {uu)i+ 


{us)i+ 




{uu)i- 


{us)i- 


(ss)i- 


niqq (GeV) 


0.74 


0.88 0.95 


1.05 


1.13 


1.47 


1.53 


1.64 



(2.3.42) 



The mass ordering is characteristic and model- independent, and indicates that an accurate study of 
the nucleon should retain the scalar and pseudovector correlations: {ud)Q+, {uu)i+, because for these 
fT^qq ^^N- Thls cxpcctatiou is verified in calculations, where one finds [p!5| , |152| , p!5^ that including the 
{uu)i+ correlation yields a nucleon whose mass is a welcome ~ 33 % less-than that of a quark+scalar- 
diquark-only nucleon. We note that m(^ud)g+/^(uu)-^+ = 0.78 cf. 0.76 = mj^/rri/^ and hence one might 
anticipate that the presence of diquark correlations in baryons is likely to provide a straightforward 
explanation of the A^- A mass-splitting. Lattice estimates, where available ||155|| , agree with these results. 



Reference | |154| | is an extensive study of the octet and decuplet baryon spectrum based on a quark- 



diquark Fadde'ev equation. It represents the nucleon as a composite of a quark and pointlike diquark, 
which are bound together by a repeated exchange of roles between the dormant and diquark-participant 
quarks, as depicted in Fig. and employs the confined-particle representation of the dressed-quark 
and -diquark propagators advocated in Ref. ||114|| . Seven parameters appear in the model and, un- 



fortunately, variations are permitted in some of them that allow the calculated results to override 
intuitive expectations. For example, the A^-A mass-splitting is fitted by adjusting the relative strength 
of the quark-diquark couplings in the scalar and pseudovector channels whilst simultaneously enforc- 
ing m(^ud)Q+ = ^{uu)j^+- Bethe-Salpeter equation studies show that this is erroneous; i.e., these masses 
cannot be equal, and they and the couplings are not independent. Hence this means of generating 
the A^-A mass-splitting is unlikely to be completely correct. Nevertheless, the importance of the study 
is a demonstration that an accurate description of the spectrum is possible, and indeed no calculated 



34 



mass is more-than 1% away from its experimental value. Improvements are now possible and they will 
build on the lessons Ref. ||154|| provides; e.g, replacing undetermined parameters with values obtained 
in precursor calculations. This is an important frontier. 

Even before such thorough relativistic bound state calculations, the notion that diquark correlations in 
baryons could be significant found support in the analysis of scattering observables ||156|| . As remarked 
earlier, these observables delve deep into hadron structure and so are a perennial focus of theory and 
experiment. With a growing understanding of the nature of diquark correlations, detailed nucleon 
models are now being applied in the nonperturbative evaluation of nucleon structure functions; e.g., 
Ref. ||157|| , and nucleon form factors ||158| , |159| , |16CI| , |161| , |162|| . As an exemplar, we review a calculation 
of the nucleon's scalar form factor, which also yields the nucleon a-term ||162| |. This is an important 
application because it illustrates the only method that allows an unambiguous off-shell extrapolation 
in the estimation of meson-nucleon form factors and thereby provides an analogue for the discussion of 
vector meson dominance, Eqs. ( p.3.34[ - p.3.40|) . 



Underlying Ref. ||162|| is the observation that nucleon propagation is described by a 6-point Schwinger 
function 

GIUR - R') ■■= {^liR)^iiR')) , (2.3.43) 



where 



KiR) 
KiR) 



3 

/ n ^"^^i ^(^i ~ R'^ "i' ^i'^ ^)^abc Qaixu tti, Ti) qb{x2; ^2, Qcixs; Q^s, T3) , (2.3.44) 

i=l 
3 

J n ^"^^i ^(^i ~ R'^ ^i'^ ^y^abc qa{Xu "l, n) qb{x2] 02, ^2) qc{xz] 0^3, n), (2.3.45) 



with: qa{xi] ai, ti), etc., quark Grassmann variables; ttj, a quark and nucleon Dirac subscripts; Xj, r the 
isospin analogues; and Sahc ensuring colour neutrality. In these expressions '?/'(xj — i?; Oj, rj; a, r) describes 
the distribution of quarks in the nucleon and; e.g., it can represent a nucleon Fadde'ev amplitude. The 
electromagnetic interaction of this nucleon is described by the current 



J^{E! - Ro, Ro-R) = -{^{R') q{Ro)tQe7MRo) ^(^)) • 



(2.3.46) 



To proceed, Ref. [|162|| writes 
ip^Xi - R; ai,Ti; a,T) 



i=l 



R) 



(2vr)^ 

and employs a product Ansatz for the nucleon amplitude 

iIj{Pu ai, Ti] a, r) = ^^"'^ 5^0,.^ ipipi + ^2,^3) ^{Pi + P2) ^a^a2iPi^P2) , 



^2.3.47) 



(2.3.48) 



where ^'(^15^2) is a Bethe-Salpeter-like amplitude characterising the relative-momentum dependence 
of the correlation between diquark and quark [~ $ in Fig. |2.5|| , A{K) describes the pseudo-particle 
propagation characteristics of the diquark, and 



(2.3.49) 



represents the momentum-dependence, and spin and isospin character of the diquark correlation; i.e., 
it corresponds to a Bethe-Salpeter-like amplitude for what here is a nonpointlike diquark. Complete 
antisymmetrisation is not explicit in this product Ansatz but that is effected in ^ via the contraction 
with the Grassmann elements, Eq. ( [2.3.44 ). 
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Fig. 2.6. Using tlie product Ansatz of Eq. ( p.3.48| ), tlie impulse approximation to the 
electromagnetic current, Eq. ( p.3.46 ), requires the calculation of these five contributions, ip: 
ip{(.iA2) in ( [2.3.481 ); V: Bethe-Salpeter-like diquark amplitude in ( [2.3.49| ); dotted line: ^{K), 
diquark propagator in ( |2.3.52|) ; solid hne: S{q), quark propagator in ( |2.3.55| ). The lowest 
three diagrams, which describe the interchange between the dormant quark and the diquark 
participants, effect the antisymmetrisation of the nucleon's Fadde'ev amplitude. Current 



conservation follows because the photon-quark vertex is dressed, given in ( 2.3.31 



The impulse approximation to the electromagnetic current in this model can now be obtained directly 
from Eq. ( [2.3.461 ). This 8-point Schwinger function expresses a product of eight Grassmann variables 
and, via an analogue of Wick's theorem, that can be reduced to a sum of products of four 2-point dressed- 
quark Schwinger functions and the dressed-quark-photon vertex, with momenta and indices correlated 
via Eq. ( [2.3.48[ ). Using the 1^2 particle exchange symmetry exhibited explicitly by Eq. ( p.3.48[ ) one 
arrives at the result depicted in Fig. ^]6[ 



The product Ansatzis completely specified once explicit forms for the functions are given and Ref. | ]162 | 
employs 



^(^1,4) 
r(gi,g2) 
AW 

Hv) 



X 
1 



2 HK^u^l) , 



1 -e 



y 



1^1 



ill 



1^2, 



\(12 



(2.3.50) 
(2.3.51) 
(2.3.52) 

(2.3.53) 



whose parameters were fixed by fitting the proton's charge form factor on Q'^ E [0, 3] GeV^ |[160[| : 









l/cjr 


l/rriA 


in GeV 


0.20 1.4 0.63 ' 


in fm 


0.99 


0.14 


0.31 



^2.3.54) 
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The two tables demonstrate the internal consistency of the model, dr := l/ur is a measure of the mean 
separation between the quarks constituting the scalar diquark and := l/u^ is the analogue for the 
quark-diquark separation, dr < d^ is necessary if the quark-quark clustering interpretation is to be valid. 
^{ud)^+ =1/ '^A is a measure of the range over which the diquark persists and that must be significantly 
less than the nucleon's diameter. [M^ and A/r are the nucleon and {ud) diquark normalisation constants, 
which are defined and calculated analogously to Eq. ( p.3.6| ), and ensure composite electric charges of 1 
for the proton and 1/3 for the diquark.] 

As is plain in Fig. [^.6| , the calculation also involves the dressed-quark propagator and, based on the 
success of applications such as Refs. |]38| , |8^ , p!05| , |106| , pTT|J14CI| , p^l| , the following algebraic parametrisation 
is employed in Refs. | ]160| , |162[ : 



S{p) 
as{x) 

avix) 



2 m J^(2(x + m^))+ J^{bix) J^{hx) [bo + &2^(ex)] 



X + m 



m 2 



(2.3.55) 
(2.3.56) 

(2.3.57) 



X 



m/A, crs{x) = Xasip"^) and av{x) = X^avip^)- The mass-scale, A = 0.566 GeV, and 



parameter values 



hi 



0.00897 0.131 2.90 0.603 0.185 



:2.3.58) 



were fixed in a least-squares fit to light-meson observables ||137|| . [e = 10~^ in ( |2.3.56| ) acts only to 
decouple the large- and intermediate-p^ domains.] This simple form captures the essential features of 
direct solutions of Eq. ( |2.2.1| ). It represents the dressed-quark propagator as an entire function, which 



is inspired by the algebraic solutions in Refs. ||163| , p.64|| and ensures confinement via the means described 
in Sec. B]^; and exhibits DCSB with 



GeV^ 



A^ 



47r2 hi 63 



In 



^^QCD 



(0.221 GeV)= 



;2.3.59) 



which is calculated directly from Eqs. ( p^.2.12|) , ( p^.2.48|) after noting that Eqs. (|2.3.56| ), (|2.3.57| ) yield 
the chiral limit quark mass function of Eq. ( p.2.47| ) with 7^ = 1. This is a general feature; i.e., 
the parametrisation exhibits asymptotic freedom at large-p^ omitting only the additional \\i{p'^ / K^qj^)- 
suppression, which is a useful but not necessary simplification. As we see here, this omission introduces 
model artefacts that are easily identified and accounted for. 



Reference ||160|| obtains a good description of the proton's charge and magnetic form factors, and also 
the neutron's magnetic form factor. The nonpointlike nature of the scalar diquark correlation makes 
possible a magnetic moment ratio of |yU„//ip| = 0.55, cf. 0.68 experimentally. This ratio is always less- 
than 0.5 when the correlation is pointlike. The neutron's charge form factor is poorly described. The 
charge-radius- squared is negative, consistent with the data, but is 60% too large. That defect results 
primarily from neglecting the contribution of the axial- vector diquark. 



The nucleon's scalar form factor is 

(T(g2) u{P')u{P) := (P'l m{uu + dd) \P) , 
where the nucleon spinors satisfy: 

7 ■ Pu{P) = iMu{P) , u{P) -fP = iMu{P) 



(2.3.60) 



(2.3.61) 
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with the nucleoli mass M = 0.94 GeV, q = (P'-P) and R = (P' + P). The vrA^ a-term is just aiq"^ = 0), 
which is the in-nucleon expectation value of the explicit chiral symmetry breaking term in the QCD 
Lagrangian. The general form of a fermion- scalar vertex is 

\{q,P) = h+t^-qf2 + t7-Rf3 + tcr,.R,q, A , /. = R') (2.3.62) 

since q ■ R = for elastic processes. However, using Eqs. ( |2.3.61| ) the scalar current simplifies: 

Jl{P',P) := u{P')A-^_{q,P)u{P) = s{q')u{P')u{P), (2.3.63) 



h-2Mh + q^U. 



;2.3.64) 



The impulse approximation to J-^(P',P) is also given by the five diagrams in Fig. ^]6| but with the 
dressed-quark-photon vertex replaced by the dressed-quark-scalar vertex, which is the solution of an 
inhomogeneous BSE analogous to Eq. (|2.3.29|) . In Ref. ||162|| , to hasten an exemplifying result, the scalar 
vertex equation was solved using the Goldstone-theorem-preserving separable model of Ref. ||114|| . In 
that model the BSE assumes the form 



d^q 



A(fc-g)7M%+)ri(g;Q)^(g-)7M 



(2.3.65) 



with the interaction 

A(/t -q) = G{k^)G{q^) + k ■ qF{e)F{q^) , (2.3.66) 

where Fiji^)^ G{k'^) are regularised forms of the A, B obtained from Eqs. ( p. 3. 561) , ( |2.3.57|) . [NB. 
The regularisation ensures convergence of necessary integrals in the separable model.] The solution of 
Eq. (|2X65| ) is 



F-, {k;Q) = l + h{Q'')G{k'')+ it2iQ^)F{e 



k-Q^-Q 



+ ih{Q^)F{k'')^-k. 



(2.3.67) 



The cr-term is only sensitive to the vertex at = 0, where the solution reduces to 

F^ {k; Q)|q2=o = 1 + ti(0) G{e) + t3(0) F{e) Z7 ■ A: , (2.3.68) 

with ti(0) = 0.242 GeV, ^3(0) = -0.0140 GeV. Calculation shows that at A;^ = the ti-term is 6- 
times larger than the bare term; i.e., it is dominant in the infrared. That is to be expected because 
it represents the effect of the nonperturbative DCSB mechanism in the solution. This and the other 
tj-terms vanish as fc^ 00, which is just a manifestation of asymptotic freedom. 

What makes the inhomogeneous scalar BSE important here is that it has a solution for all and that 
solution exhibits a pole at the cr-meson mass; i.e., in the neighbourhood of {—Q^) = m'^ = (0.715 GeV)^ 

2 

Tl 771 

T-^_{k; Q) = regular + " \ Y ,(k- Q) , (2.3.69) 

LJ + ui^ 

where reg'w/ar indicates terms that are regular in this neighbourhood and T(j{k; Q) is the canonically nor- 
malised cr-meson Bethe-Salpeter amplitude, whose form is exactly that of (F-| ( fc; Q) — l) in Eq. (|2.3.67|) . 
This is qualitatively identical to the behaviour of the solution of Eq. (|2.3.29| ) discussed above in con- 
nection with F-^i^Q"^) and elucidated in Refs. ||10?| , |131|| . The simple pole appears in each of the functions 
ti{Q'^) and a pole fit yields, with m the current-quark mass, 

mn^ = 3.3 MeV. (2.3.70) 
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rio-m^ is the analogue of the residue of the vr-pole in the pseudoscalar vertex, Eq. ( p.3.22| ): — ^/2 {qq)^/ fw, 
and its flow under the renormalisation group is identical, mn^j is renormalisation point independent 
and its value can be compared with 



—m 



^_ = 3.6MeV; 



(2.3.71) 



i.e., the magnitude of n^- is typical of effects driven by dynamical chiral symmetry breaking. The 
cr-meson-gg coupling is defined on-shell 



9aqq ■= T^{0;Q)\q2 



12.6, 



and its magnitude can be placed in context via a comparison with g^^gg 
analogously. 

The expectation value in Eq. (|2.3.60| ) is obtained with 

Tmik;Q) = mT-^{k;Q) 



(2.3.72) 

11.8, which is defined 



(2.3.73) 



as the probe vertex in Fig. |2.6| , and using the solution for F-^ described here the calculated a-term 
is M 

(x/Mjv = 0.015 . (2.3.74) 

No parameters were varied to obtain this result, which may be compared with a recent lattice com- 
putation ||165|| : ct/Mat = 0.019 ± 0.05, calculated using an extrapolation in the current-quark mass. 
Alternative extrapolation methods can lead to larger values; e.g., ct/Mtv = 0.047-0.059 ||166|| , which are 
also suggested by some phenomenological analyses P^ , |167| |. [The a-term is not directly accessible via 
experiment but a value is theoretically inferred by extrapolating vrA^ scattering data using dispersion 
relations ||168| ]: a/M^ = 0.047-0.076.] Simple estimates indicate that the result in Eq. ( |2.3.74| ) increases 
with decreasing rricr, and a reduction in ma is a likely consequence of using an improved kernel in the 
inhomogeneous scalar BSE. Thus Eq. (|2.3.74| ) is an excellent first estimate. 

The nucleon's scalar form factor is depicted in Fig. |2.7| , where the evolution to the cr-meson pole is 
evident. Fitting {t = —q^) 



(^{t) = 9(7 NN 



t G [0.1,0.5] GeV^ 



(2.3.75) 
27.3. A 



which isolates the residue associated with Fm(fc; Q), one obtains the on-shell coupling: q^nn 
direct calculation using the solution of the homogeneous Bethe-Salpeter equation yields QaNN = 27.7, 
in agreement within Monte-Carlo errors. Equation (|2.3.75| ) alone overestimates the magnitude of the 
calculated a{t) everywhere except in the neighbourhood of the pole. 

As the lowest-mass pole-solution of a rainbow-ladder BSE, Eq. ( |2.3.65| ), the a-meson described here 
is not obviously related to the scalar meson introduced in phenomenological nucleon-nucleon poten- 
tials [ |169|J170|] and meson exchange models ||1 71|] to mock-up two-pion exchange. However, a coupling 
relevant to such models can be estimated by introducing gcr{t): 



a{t) =: gait) 



mnn 



1 — t/ 



where a fit to the calculated a(t) in Fig. yields 

gait) = 1.61 + 2.61- 



2\10 



Aa = 1.56 GeV. 



(2.3.76) 



(2.3.77) 



39 



0.020 



0.015- 



> 

X 0.010 



0.005 - 



0.000 




Fig. 2.7. cr(g^) as calculated in Ref. [162] using the solution of the inhomogeneous scalar 
BSE, Eq. ( p^.3.65D . The rapid increase with decreasing is associated with the evolution 
to the (T-meson pole. A similar feature was also encountered in the calculation of the pion's 
electromagnetic form factor, see Fig. p.4| on page |3l| and the associated discussion. On this 
scale, cr(g^) calculated without the t2,3(<5^) contributions is indistinguishable from the full 
calculation. (Adapted from Ref. ||162|| .) 



with the large exponent merely reflecting the rapid evolution from bound state to continuum dominance 
of the vertex in the spacelike region. ga{t) describes the t-dependent nucleon to scalar-quark-antiquark- 

0.5 GeV, 



correlation coupling strength and at the mock-a-mass: m; 



2-K 



9.3 



(2.3.78) 



which may be compared with a phenomenologically inferred value: g^j = 10 ||172|| . This gives a mi- 
croscopic interpretation of the mock-a. [It is important to note that g^lAmD = 5.2 and hence this 
comparison is meaningful on a relevant phenomenological domain.] Further, g^lq"^ — > oo) = 1.61 so that 
(T(g^) is well approximated by a single monopole for > 1 GeV^. However, the residue is very different 
from the on-shell value. The scalar radius of the nucleon is obtained from 



6 (icr(g^ 



which may be compared with an inferred value 



= (0.89 fm)^ 
: (1.2fm)2. 



(2.3.79) 



Reference ||162|| calculates a range of non-electromagnetic nucleon form factors and provides a uniformly 
good description, yielding meson-nucleon form factors that are "soft" and couphngs that generally 
agree well with those employed in meson-exchange models. Where there are discrepancies with these 
models or with experiment, a plausible cause and means for its amelioration was readily identified. 
The study demonstrates that it is realistic to hope for useful constraints on meson-exchange models 
from well-moulded models of hadron structure. The analysis of the nucleon a-term is particularly 
important because it illustrates the only method that allows an unambiguous off-shell extrapolation in 
the estimation of meson-nucleon form factors: one must employ solutions of the inhomogeneous BSEs 
to describe correlations in a channel when the total momentum is not in the neighbourhood of a bound 
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state pole. Improvements to the simple quark+scalar-diquark nucleoli model are being developed and 
their application to a wide range of hadronic phenomena is being explored. These studies hold the 
promise of providing a much-needed standard model of the nucleon. 



3 Nonzero Density and Temperature 

Identifying the onset and properties of a quark gluon plasma [QGP] is the primary objective of studies 
at nonzero temperature and density. Such a plasma is expected to be characterised by light quarks and 
gluons propagating freely over distances ~ 10-times larger than the proton. The exploration of QCD in 



this domain requires a knowledge of equilibrium statistical field theory, for which Sec. 6 of Ref. |105| is 



a primer and Refs. ||173[ provide a graduate-level introduction. Here the application of DSEs provides 



our primary medium and we will describe the necessary and related concepts. 

Section previews the application of DSEs at zero density and temperature. It shows that the framework 
is efficacious on this domain and well constrained. Both are important because little is known about the 
domain of high density and/or temperature, and much of what we do know comes from the extrapolation 
of models. Since the DSEs provide a nonperturbative framework that admits the simultaneous study 
of DCSB and confinement they are well suited to explore the phase transition that yields the QGP. 
Furthermore, as they also accurately describe bound states, they can be used to explore the response of 
hadron properties to extremes of density and temperature; i.e., to elucidate signals of QGP formation. 
These features also mean that the DSEs can address all the phenomena accessible in lattice-QCD 
simulations and therefore they can be used to confirm the results of such studies. Importantly too 
the lattice simulations can be used to constrain the model-dependent elements of DSE studies, making 
possible a reliable DSE extrapolation into the domain of nonzero chemical potential QCD, for example, 
which is inaccessible in contemporary lattice-QCD. 



3.1 In-medium Essentials 

Equation of State The T ^ analogue of the partition function in Eq. ( |2.1.9| ) is 

Z^^^T)[V,V,J] = J dnp{q,q,A,uj,u) exp^ dr J £x [qv + V Q + , (3.1.1) 
with X4 ^ T, f3 := 1/T and the measure 

diJ,i3{q,q,A,u;,uj):= (3.1.2) 
n n^^0(^' r)Vq^{x, t)1[Vlj''{x, t)Vuj\x, t)X{VAI{x, t) exp(-5^[g, g. A] - S^p, uj, A]) , 

XTe[0,/3] <i> a n 



where Splq, q, A], S^[u, u, A] are the actions of Eqs. ( p.l.8|) , ( p.l.l2|) but for the obvious bounding of the 



r-integral. The frame specified by separating the spacetime integral into a product of a compact integral 
over Euclidean time and the volume integral is not prescribed by the theory. This separation breaks 
the theory's original 0(4) symmetry to 0(3) and is effected by introducing an arbitrary, normalised, 
spacelike vector, u^, with u-p = defining the T = 0-hyperplane. Here we have employed a conventional 
choice for the heat bath vector: u = (0, 1). [NB. The zero temperature theory is recovered by taking the 
limit T — i> 0; i.e., /? — > 00, in Eqs. ( p.l.l| ), ( |3.1.2| ).] The chemical potential is introduced as a Lagrange 
multiplier 

Sp[q,q,A]^S^^fi[q,q,A]:=Sf3[q,q,A]-J^ dr J d^x n qj^q . (3.1.3) 
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With /i 7^ 0, S^^i3[q,q, A] is not Hermitian and hence (i/i^,/3 is not a probabihty measure. That makes 
numerical simulations of lattice-regularised /i 7^ QCD very difficult; in fact, there is currently no sat- 
isfactory algorithm ||174|| . The functional integral in Eq. ( |3 . 1 . 1| ) is evaluated on the space of Grassmann 
fields for which q{x, 0) = —q{x, i.e., over fields satisfying antiperiodic boundary conditions, and over 
gauge fields satisfying periodic boundary conditions. These boundary conditions are an essential aspect 
of the difference between fermion and boson statistics. 

The partition function is particularly important in equilibrium statistical field theory. As usual, it yields 
the DSEs in a straightforward manner, however, more than that, all the thermodynamic functions are 
obtained from the partition function. The thermodynamic potential and pressure densities are 



-^(/i,T) =p(/i,T) 



/3V 



InZ 



(3.1.4) 



where f3V is the four-volume normalising factor. The stable phase of the system is that in which the 
pressure is maximal or equivalently the thermodynamic potential energy is minimised. The expression 
for the pressure is the Equation of State [EOS] and from that one immediately obtains the baryon 
number and entropy densities: 



and also the energy density: 



d_ 

djji 



T fixed 



d_ 
df 



fi fixed 



T) = -p(/i, T) + fipifi, T) + T s(/i, T) . 



(3.1.5) 



(3.1.6) 



The calculation of these quantities in QCD is a contemporary focus; e.g., Refs. ||175| , |176| , |177 



A simple estimate of the pressure due to dressed-quarks is obtained via the "steepest descent" approx- 
imation, which yields 

(3.1.7) 



Pe(/x, T) = ^ {TrLn [(3S-'] - iTr [S S 
where S is the solution of the gap equation with S the associated self energy, and "Tr" and "Ln'' 



are extensions of "tr" and "In" to matrix- valued functions. Equation ( |3.1.7D is just the auxiliary field 



effective action ||178|| , which yields the free fermion pressure in the absence of interactions; i.e., when 
S = 0. At this level of truncation the total pressure receives an additive contribution from dressed- 
gluons: 



PA(/i,T) 



pv 



^TrLn 



(3.1.8) 



where D^,^ is the T 7^ dressed-gluon 2-point function, and this yields the free-gluon pressure in the 
absence of interactions. 



Propagators As just described, the introduction of temperature to Euclidean QCD breaks the original 



0(4) symmetry to 0(3). In this case the most general form of the dressed-quark propagator is ||1 79 
S{p,u)~'^ = i'j-pa{p^, [u-pf^) + h{p^, [u-pf) + n ■uu-pc{p'^, [u ■ pf) + ia^vp^Uud{p^ , [u-pf) , (3.1.9) 



where = (0, 1), and the surviving 0(3) invariance of Eq. ( p.l.9|) is apparent upon inspection. Even 
at T 7^ 0, the chiral-limit dressed-quark propagator should satisfy 



V{a) S{p, u)-^ V{a) = S{p, u)-^ , V{a) = exp (^5 ^X^a^) 



(3.1.10) 
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if chiral symmetry is realised in the Wigner-Weyl mode. Hence b = = d in this case. Furthermore, 
since the term involving d does not appear in the free fermion action of thermal field theory, there are 
no perturbative divergences in this function's evaluation and d is therefore power-law suppressed in the 
ultraviolet, [cf. The discussion of the chiral limit fermion mass function on page ^] These features, 
along with the vector exchange nature of the interaction, underly the fact that d{p'^,u -p^) only plays a 
minor role in the solution of the QCD gap equation. Hereafter we neglect it and assume the following 
as the most general form of the dressed-quark 2-point function: 

^^^'"""^ = t^-pA{P,Col) + B{p\^l)+t^,Co,C{p\ul) ' ^^-^-^^^ 
= -i^-paA{p'^,tbl) + aB{p'^,i^l)-ilitbkCrc{p'^,ujl), (3.1.12) 

where ujk = ojk + with /i the chemical potential, and ujk = {2k + 1) vrT, k ^ Z, are the fermion 

Matsubara frequencies, which ensure q{x,0) = —q{x,/3). The scalar functions: = A, B, C, are 
complex and satisfy 

J^ip^ulr=Tip',u\_,). (3.1.13) 

[NB. The functions A, C cannot be neglected. For T = 0, A = C and A - 1 ^ 0; e.g., Refs. 
and the momentum dependence of A can conspire with that of B to ensure confinement of dressed- 
quarks [|164|| . Furthermore, as will become clear, the (p^, ci)^)-dependence of A, C has a marked effect 
on the behaviour of bulk thermodynamic quantities.] 

The dressed-quark 2-point function satisfies a gap equation 

^(PajJ^^ = ^^^7 ■P+ ^2 (^574^^/0 +"^bare) + S'(p^J , (3.1.14) 

which can be derived in the usual way from the partition function, Eq. ( p.l.l|) . Here we have introduced 
a shorthand notation: j9^^ = {p,u}k), and the regularised self energy is 

S'(Pc.J = -pT^'^iPu^J + hiUJkT.'c{Pu^J + T.'b{Pc.J , (3.1.15) 
^'APujJ = itrP^ / jg^D^^{p-q,uJk-uJi)ji,S{q^^)T^{q^^]p^J, (3.1.16) 

Jl,q 

where the renormalisation factors are Va '■= —{Z^ /\p\'^)ij ■ p, Vb '■= Zi, Vc '■= —{Zi/uJk)i'j4,, and 
:= T J2iZ-oo {'^'^Y with representing a translationally invariant regularisation of the 

three-dimensional integral and A is the regularisation mass-scale. Since introducing /i, T does not 
generate qualitatively new divergences, the regularisation and renormalisation proceeds just as in the 
absence of the medium so that the renormalised self energies are 

HVu^k-X) = e^ + S^(p^,;A)-E^(C-;A), (3.1.17) 

with C the renormalisation point, (Cj^)^ := ~ ^O' = ^ = ^c, and C,b = rriR^Q. [cf. Eqs. (p.2.39| - 



The regularised self energy is expressed in terms of the renormalised dressed-gluon 2-point function and 
dressed-quark-gluon 3-point function. At T 7^ in Landau gauge the complete expression of the former 
requires two 0(3)-scalar functions [cf. Eq. (|2.2.15| )] 

g^D^,{p^^) = P;:^{pn,)Ap{pn,) + P^^ipnJAcipnJ , (3.1.18) 
where p^^. := (p, Qk), = 2A;7rT is the Matsubara frequency for bosons, and 

{0, /i and/or z/ = 4, 

S,-'f, = =1,2,3, 

P'M = S,^-^-P^APnJ. (3.1.20) 
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In the absence of interactions: Ap^pj^J = Ac^pfi^) = The QCD analogue of a Debye [or electric 

screening] mass appears as a T-dependent contribution to Ap, the gluon's longitudinal polarisation 
function. The complete expression of the renormalised dressed- quark-gluon 3-point function requires at 
least fifty-four scalar functions [cf. twelve for T = 0] but, since the theory is renormalisable, ultraviolet 
divergences are only encountered in evaluating two of them. It follows from this that sensible, character- 
preserving truncations are possible. 



Order Parameters In order to demarcate the QGP phase it is necessary to identify order parameters; 
i.e., find operators Xi whose expectation values are nonzero in the normal phase: (Xj) 7^ 0, but vanish 
in the QGP: (Xj) = 0. This is often a difficult task. Having identified such operators then, if (Xj) is 
discontinuous at the transition point we will describe the transition as first order; otherwise we classify 
it as second order. 



In Sec. p.2| we described the phenomenon of DCSB and its connection with the purely nonperturbative 
generation of a scalar term in the dressed-quark self energy when the current-quark mass is zero. One 
characteristic of this effect is the large magnitude, Eq. (|2.2.51|) , of the vacuum quark condensate that 



appears in Eqs. ( 2.2.47 ), ( |2.2.48 ), which has many observable consequences; e.g., it is the primary cause 



of what is naively an unexpectedly large vr-p mass splitting. DCSB is a defining feature of the normal 
phase of QCD. If, with increasing fi and/or T, the condensate vanishes at some particular value of 
these parameters then a new phase of QCD has been reached. It is a phase in which chiral symmetry 
is realised explicitly and, as weakly interacting systems cannot support condensate formation, this is 
plainly a good defining property of a QGP. The discussion of the chiral limit, Eqs. ( |2.2.45[ - p.2.52D , also 
makes plain that it is not possible to have a nonzero quark condensate unless the dressed-quark mass 
function, M{p^^), is nonzero too. Therefore the magnitude of this function at any single point is also 
an order parameter appropriate for characterising the chiral symmetry restoring aspect of the QGP 
transition. In fact, in DSE studies it is by far the most direct, since it appears in all the integrands 
that describe other quantities. 

We also discussed an arguably more significant aspect of cold, sparse QCD, namely confinement. It is 
easy to build a model that exhibits DCSB without confinement; e.g., Refs. p!^Jl3|JI^jl80|| . However, 
building a covariant model with confinement and an uncomplicated realisation of all aspects of DCSB 



is a challenging task. As emphasised by the study of QED3 in Ref. p8[, that challenge is met with 
a modicum of success using the DSE truncation scheme explored in Ref. | 108 | and the realisation of 
confinement via a violation of refiection positivity by coloured n-point functions. This also exposes a 
simple deconfinement order parameter appropriate to both light- and heavy-quarks, whose application 



is illustrated in Fig. |2.3| . To make it plain we return to D{T) defined in Eq. (|2.2.72| ). When calculated 



from a Schwinger function with complex conjugate poles, D{T) has at least one zero. Denote its position 
by Tq^ and define 

no:=l/T^\ (3.1.21) 

If, with increasing /i and/or T, one observes kq approaching zero then the [first] zero in D{T) is moving 
to T = 00. This evolution of corresponds to the limiting case when there is no zero at all, which is 
just the free, unconfined particle described by (JsiT) in Eqs. (|2.2.671 ). Thus kq is a good order parameter 



for deconfinement ||181|| : it is nonzero in the confined phase, and remains nonzero until the Schwinger 



function evolves into one that respects refiection positivity and corresponds to an asymptotic state. It 
is clear that deconfinement is also a good defining property of a QGP. 



Critical Behaviour In second order transitions the length scale associated with correlations in the 
system diverges as the order parameter vanishes and a range of critical exponents can be defined that 
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characterise the behaviour of macroscopic properties at the transition point. In the context of QCD 
this can be elucidated via the free energy, which here we write as f{t,h), where t := T/T^ — 1 is the 
reduced temperature, with Tc a putative critical temperature, and h := jSnifj is the source of explicit 
chiral symmetry breaking measured in units of the temperature. [We omit the chemical potential for 
now so that / = e = {T"^ /V){d In Z(^T)/dT).] h plays a role analogous to that of an external magnetic 
field applied to a ferromagnet. Since correlation lengths diverge it follows that for t,h ^ the free 
energy is a generalised homogeneous function; i.e., 

fit,h) = ]-f{tby\hF'^). (3.1.22) 



As a consequence the "magnetisation" behaves as follows: 

, M{t,h) = W'^~^ M{tW\hW'^) . (3.1.23) 

t fixed 

[NB. Using the partition function one finds easily that M(t, h) is just the vacuum quark condensate. 
In these formulae it can be replaced by any equivalent order parameter.] The scaling parameter: 6, is 
arbitrary and along the trajectory = 1 one has 

M(t, h) = M(sgn(t), h Itl-y^/y^) ■ i.e., M{t, 0) oc {tl'^ , /3 := . (3.1.24) 

Vt 

Alternatively, along the trajectory hby^ = 1 

M{t, h) = h^^-y^'^ly^ M{t h-y''y\ l) ; i.e., M(0, h) oc h^'^ , 5 := . (3.1.25) 



M{t, h) 



df{t,h) 
dh 



Equations ( p. 1.241) , (|3.1.25| ) quantify the behaviour to be expected of an order parameter at a second 
order transition. They also introduce two critical exponents and, using the renormalisation group, 
scaling laws can be derived that relate all other such exponents to (3 and 6 ||182||. It is widely conjectured 



that the values of these exponents are fully determined by the dimension of space and the nature of 
the order parameter. This is the notion of universality; i.e., that the critical exponents are independent 
of a theory's microscopic details and hence all theories can be grouped into a much smaller number of 
universality classes according to the values of their critical exponents. If this is the case, the behaviour 
of a complicated theory near criticality is completely determined by the behaviour of a simpler theory 
in the same universality class. In mean-field theories 



P 



MF 



1 

2 ' 



c-MF 



3.0, 



(3.1.26) 



The success of the nonlinear cx-model in describing long- wavelength pion dynamics underlies a conjec- 
ture ||183|| that chiral symmetry restoration at T 7^ in 2-flavour QCD is a second order transition with 
the theory lying in the universality class characterised by the 3-dimensional, = 4 Heisenberg magnet 



[0(4) model] whose critical exponents are ||184|| : 

= 0.38 ± 0.01 , 5^ = 4.82 ± 0.05 , 



(3.1.27) 



As an alternative to Eqs. ( p.l.24|) , (|3.1.25|) , the critical exponents can be determined by studying the 
pseudocritical behaviour of the chiral and thermal susceptibilities: 



Xhit,h) :-- 



dM{t, h) 



Xt{t,h) 



t fixed 



dM{t, h) 



(3.1.28) 



h fixed 
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which are smooth functions of t with maxima at the pseudocritical points: tp^, and t^^, where tp^, oc 
fi^KPi) oc tp^. Since (35 > 0, the pseudocritical points approach the critical point: t = 0, as h ^ 0^ and 



xT ■= Xhif;,, h) oc h-'^ , := 1 - 1 , xr ■= Xtitl,, h) oc h-'^ ^ Zt:=^{l-/3). 
[The appendix of Ref. [ 185|| provides a demonstration.] Mean field behaviour corresponds to 



3 ' 



MF 



(3.1.29) 



(3.1.30) 



Analysing the susceptibilities can provide more accurate results when numerical noise is anticipated, 
which is why the method has been used in lattice-QCD analyses [ |175|| . 



3.2 Quark Gluon Plasma Phase 

Temperature We are now in a position to describe the application of DSEs to demarcating the QGP 
and exploring its properties. For the moment we persist with yU = and explore the temperature induced 
transition because that is the domain on which the most complete studies exist [in all approaches]. 



The properties of the class of rainbow-ladder models are elucidated in Ref. ||186|| . This class is defined 
by Eq. (|2.2.33| ) with the mutually consistent T 7^ constraints: Zi = Z2, Zf = Z^. [NB. The rainbow 



truncation is the leading term in a 1/iVc-expansion of Tj,{q^^;p^^).] The form of D^^{pq^) has no bearing 
on whether or not a model is in the class. To be concrete, Ref. ||186|| considered models defined by 



^F{pn,) = Vipn,;mg), AGipnJ=-D{pn,;0), (3.2.1) 
V{pn,;mg) := 2Tr^ D ^5okS%p} + VMiPn,;mg) , (3.2.2) 

with three choices of Vyi^po^^^mg). The first term in V{pQ^]mg) is a simple T 7^ generalisation of 
the distribution in Eq. ( ^.2.351 ). It represents the long-range piece of the effective interaction and the 
T^yiip^k] t^g) are chosen so as not to pollute that. (NB. does not contain a simple T-dependent mass. 
Such "magnetic" screening is more complicated than electric screening because of infrared divergences 
in QCD ||173|| .) The three choices are: 

A) I^M:=A(pQfc; ^g) = 0, with D := rj^ /2 and i] = 1.06 fixed [^] by fitting vr- and p-meson masses 
at T = 0. m = 12 MeV yields m^r = 140 MeV. T>a yields an ultraviolet finite model and hence the 
renormalisation point and cutoff can be removed simultaneously. Nevertheless it exhibits many features 
in common with more sophisticated Ansdtze 



B) The model of Ref. [|T81|], obtained with D := (8/9) mf and 



1 _ e-*"fc/(4m?) 

VM:=B{Pn,;rng) = fn^ , (3.2.3) 



s^^ ■= pI^ + ml, where, following Ref. ml = 47r^-f^{Nj3 + Nf/6) = (8/3) vr^T^ [Nf = 3]. The 



mass-scale rrit = 0.69 GeV = 1/0.29 fm marks the boundary between the perturbative and nonpertur- 
bative domains, and was also fixed [H8| by fitting vr- and p-meson properties at T = 0. At C = 9.5 GeV, 



mji = 1.1 MeV yields = 140 MeV. This model adds a Coulomb-like short-range interaction to that 
of Ref. [ p.76| ], thereby marginally improving its ultraviolet behaviour. 

C) A minimal T 7^ extension of the effective interaction in Eq. (|2.2.35| ) with the only modification 
being k"^ — > s^^,, = (16/5) vr^T^, in the last two terms. [Recall that this model's parameters were 
fitted with Nf =4]. As described above, Vq ensures that the model preserves the one-loop renormalisa- 
tion group behaviour of QCD. To further explore model-dependence, two parameter sets were employed 
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Table 3.1. Critical temperature for chiral symmetry restoration and critical exponents char- 
acterising the second-order transition in the four exemplary models. Averaging over the 
independent models yields Tc = 154 ± 24 MeV. (Adapted from Ref. 



186] 



Model 



A 



B 



T,{MeV) 

Zt 



169 
0.666 
0.335 



174 
0.67 ±0.01 
0.33 ±0.02 



120 

0.667 ±0.001 
0.333 ±0.001 



152 

0.669 ±0.005 
0.33 ±0.01 



in Ref. |[186|| , which differ only in the value of uj: lU = 0.6 rrit, 2^ = 1-2 nit. They are equivalent, yielding 
the same values for pion observables: with lU, a renormalisation point invariant current-quark mass of 
rh = 6.6 MeV yields = 140 MeV; while with 21-1^, m = 5.7 MeV effects that. 

Two chiral order parameters were employed in analysing chiral symmetry restoration 

^ :=B(p = 0,0.0), ^c:=^^C^^. (3.2.4) 

C{p = 0,^0) 

They should be equivalent and the onset of that equivalence determines the /i-domain on which 
Eqs. ( |3.1.29|) are valid. These being bona fide order parameters is useful and particularly important 



because it means that the lowest Matsubara frequency completely determines the character of the chiral 
phase transition, as conjectured in Ref. ||188|| . Furthermore, it was numerically verified in Ref. ||186 



that in the chiral limit and for t ~ 0: oc {qq) oc X{t, 0); i.e., these quantities are all bona fide order 
parameters. 



The results for the critical temperature and exponents are presented in Table RTTl where the quoted crit- 



ical temperatures are easily determined using either Eqs. ( ^X24D , (^X25D or tjj, oc /i^/^^"^) oc t^,. The 



models are all constrained by observables in two-flavour QCD and the critical temperature agrees with 



the current estimate from simulations of lattice-QCD with two dynamical light-quark flavours ||1 75|| . 
(Quenched lattice-QCD exhibits a first order [deconfining] transition at a much higher critical tempera- 
ture: T^^'^'^'^^ = 270(5) MeV ||175| , |189| |.) Studies using simple versions of the Nambu-Jona-Lasinio model 
typically yield higher values: Tc > 200 MeV |T^. However, that is likely a model artefact. Modifica- 



tions that make possible a relaxation of the overly restrictive connection between the model's mass-scale 
and the dressed-quark mass; such as l/A^c-[or meson-loop-] corrections or replacing the contact inter- 
action by a separable form, can easily effect a reduction of the critical temperature by < 25%; e.g., 
Refs. [|190| , p.91|jl92|J193|| . This sort of reduction brings it in-line with the estimate in Table |3.1| . 

The behaviour of the pseudocritical maxima in the chiral and thermal susceptibilities is illustrated in 
Fig. |3.1| . It is plain that the results exhibit curvature and hence the scaling relations of Eq. ( |3.1.29| ) are 



not satisfied on the entire domain. The domain on which they are valid may be established by defining 
a "local" critical exponent: 

_ _ InxT^hix^ ^ .3 2.5) 

* In hi — In /ij+i ' 

where {hi, xf^) and (/ij+i, xf+i) are adjacent data pairs, h lies in the scaling region when Zi is independent 



of the order parameter. Applying this to the model yields the results in Fig. |3]^, which indicate 
clearly that the scaling relations are not valid until 

logio(VM<-7, (3.2.6) 
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Fig. 3.1. /i-dependence of the maxima of tlie cliiral (circles) and thermal (squares) suscep- 
tibilities calculated with the effective interaction of Eq. ( p.2.2|) using Vc with u = 1.2 mt. 
The filled symbols are obtained from X and the open from Xq, Eq. ( p.2.4| ). The slope of the 
straight lines is given in Table SJ. and they are drawn through the two smallest /i-values, 
which is explained in the discussion of Eq. ( |3.2.5| ) . (Adapted from Ref . 



where hu := mn/Tc is defined with the current-quark mass that gives = 140 MeV in this model. 
Applying the same method to the renormalisation-group improved models, T>c, shows that in these 
cases the scaling relations are only valid for 



lu: logio(/i//i„) < -5, 2^^: \og^o{h/hu) < -6; 
i.e., for current-quark masses six orders of magnitude smaller than those of real u,(i-quarks. 



(3.2.7) 



While these are model studies the result is also likely to be true in QCD; i.e., while the critical temper- 
ature is relatively easy to determine, very small current-quark masses may be necessary to accurately 
calculate the critical exponents from the chiral and thermal susceptibilities. If that is the case, the 
calculation of these exponents via numerical simulations of lattice-QCD will not be feasible. The 
lattice-to-lattice variation of the critical exponents described in Refs. ||175| , [189|| , with even a first order 
transition being possible, could be a signal of this. 



It is clear from Table |3.1|that each of the models considered in Ref. ||186|| is mean field in nature. As the 



long-range part of the effective interaction is identical in each case and the correlation length diverges as 
t —>■ 0, that may not be surprising. The models are distinguished by the feature that they describe the 
long-wavelength dynamics of QCD very well, with mesons represented as nonpointlike quark-antiquark 
composites. Coulomb gauge models ||187|| also describe mesons as composites and exhibit mean field 
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Fig. 3.2. (circles) and z\ (squares) from Eq. ( |3.2.5| ) calculated with the effective inter- 
action of Eq. ( |3.2.2| ) using Vb in Eq. (|3.2.3|) : filled symbols - A", open symbols - Xc- The 
dashed lines are the mean field values: Zh = 2/3, Zt = 1/3. (Adapted from Ref. ||186||.) 



critical exponents. The long-range part of the interaction there corresponds directly to the regularised 
Fourier amplitude of a linearly rising potential and therefore it is in no simple way equivalent to the 
distribution in Eq. (|2.2.35|) . Separable models with ||191j| : T^ij)^^ — qui^'mg) oc g HpI) g {\(^) , where g^pl) 
is a non-increasing function of its argument, can also provide a good description of long- wavelength 
pion dynamics. They too describe mesons as composites and exhibit an explicit fermion substructure, 
and they also have mean field critical exponents, as do models with g = g{u>l + p'^) [ |194|| . [NB. The 
Nambu-Jona-Lasinio model can be viewed as a separable model obtained with gdp]) = 6 {I — IpI/sA), 
where 3A is a cutoff parameter. Also, in the rainbow-truncation DSE model of Ref. [p.95|| , which employs 
an infrared-finite effective interaction, the critical exponent obtained from X is 1/2. However, a different 
exponent is reported to be obtained from {qq). We expect that is erroneous, as these must be equivalent 
order parameters.] 

All these examples are members of the class of rainbow truncation models. Therefore the results 
indicate that chiral symmetry restoration in this class is a mean field transition. This is an essential 
consequence of the fermion substructure in the rainbow gap equation: non-trivial critical exponents 
require the appearance of an infrared divergence in this equation. However, in rainbow truncation 
the self energy is infrared-regulated by the zeroth fermion Matsubara frequency: uq = ttT ^ in 
the vicinity of the transition. For the same reason chiral symmetry restoration at finite-T in QCD 
will be a mean field transition unless 1/iVc-corrections to the vertex are large for T ~ Tc. There are, 
however, examples in which that is certainly true ||196|| . Infrared divergences might be anticipated 



if bosonic modes dominate the gap equation near T^, because the zeroth boson Matsubara frequency 
vanishes. It is therefore important to note that pseudoscalar meson-like fluctuations are nonperturbative 
l/A'^c-corrections to the dressed- vertex. They do not contribute to the dressed-gluon 2-point function. 
Therefore a true determination of the critical exponents must await the accurate incorporation of these 
l/A'^c-corrections. 



We have seen that chiral symmetry restoration in 2-flavour QCD occurs at T ^ 150 MeV via a second 
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Fig. 3.3. Chiral and confinement order parameters calculated in Ref. | 181 | using T>b 
of Eq. ( p.2.3| ) but with an erroneous value of = Svr^T^, which is why the coincident 
deconfinement and chiral symmetry restoring transitions here occur at Tc ~ 155 MeV, cf. 
Table The solid lines are fits: X oc (— t)^/^ oc kq. 



order transition whose critical exponents are not yet reliably determined. What of deconfinement? In 
QCD with two dynamical flavours the Wilson-line/Polyakov-loop is not strictly a good order parameter 
for deconfinement because of flux tube breaking mediated by quark- ant iquark pair formation. However, 
in simulations the Polyakov loop's susceptibility develops a peak at the same temperature as the chiral 
susceptibility and this has been interpreted as evidence of coincident deconfinement and chiral symmetry 
restoration 



Ko in Eq. ( |3.1.21| ) provides an order parameter that can be used to explore deconfinement via the gap 
equation. It was exploited in Ref. ||181|| where the deconfinement transition in model "B)" was studied. 
In this case the Schwinger function analysed is 



dpp sin(pr) asip^ 



(3.2.8) 



where as is calculated in the chiral limit. Figure |3.3| depicts the behaviour of ko(T) obtained from 
(Ts{r) and clearly the deconfinement order parameter is nonzero in the normal phase. ko{T) vanishes 
at the same temperature as does A". Hence, in this model, the deconfinement and chiral symmetry 



restoring transitions are coincident [to better-than 1%] and analysis that parallels Ref. [ 186 | shows that 
deconfinement is also a mean field transition. (This improves over the critical exponents estimated in 
Ref. ||181|| .) Whether the transitions are always coincident is unknown, however, that is not unlikely 



given the role the scalar piece of the dressed-quark self energy plays in determining both transitions. 

The rainbow-ladder class of models fails to exhibit any flavour dependence in the transition to a QGP, 
a flaw it shares; e.g., with random matrix models [|188|| . It is anticipated that with three light quark 



flavours the chiral symmetry restoring transition is first order ||183|| but becomes second order if the mass 
of the "s-quark" is raised above a presently-undetermined critical value. This can be expected because 
for an infinitely heavy s-quark one recovers the two-flavour theory with its second order transition. (A 
summary of lattice results is presented in Ref. [|189|| .) A coupling between flavours is necessary to achieve 
such a pattern of symmetry breaking and, as seen above, that appears only when l/A^^c-corrections to 
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the vertex are incorporated. 



Chemical Potential The introduction of a chemical potential moves us to a domain in which the only 
available nonperturbative information comes from models. For yU 7^ the dressed-quark self energies in 
general acquire an imaginary part driven by the magnitude of /i; e.g., Refs. | p.76|J195| , [197|| . This effect 



is not observed in explorations of the Nambu-Jona-Lasinio model or indeed in any model in which the 
interaction is energy-independent; i.e., instantaneous. That follows from Eq. ( |3.1.13| ): an instantaneous 
interaction cannot yield an energy-dependent solution and hence M = M* ^ M G IR. It is an unrealistic 
limitation since the interaction in the normal [confined] phase of QCD is not instantaneous. 

An extensive exploration of the phase structure and thermodynamics of two-flavour QCD with /i 7^ 7^ 
T is presented Ref. ||176|| , which employs the "Da model. "Da is an infrared-dominant model that does not 
represent well the behaviour of D^y{p, Vt^) away from ^ and that introduces model-dependent 

artefacts. However, there is significant merit in its algebraic simplicity and, since the artefacts are easily 
identified, the model remains useful as a means of elucidating many of the qualitative features of more 
sophisticated Ansatze. Therefore we use it here as an exemplar. 

In this model the gap equation is 

5(p^J-i = 5o(p.J-^ + ir/^.5(p^j7. (3.2.9) 

and the simplicity is now apparent: the model allows the reduction of an integral equation to an 
algebraic equation, which is always a useful step when the goal is to develop an intuitive understanding 
of complicated phenomena. Substituting Eq. (|3.1.11|) yields 

2 2 
rj m 



+ mB^ + (ApI^ -rj^ - m^) B"^ - m (2 rf + + 4p^J B , (3.2.10) 
A{p.,) = C{p^J= (3.2.11) 



and it is now apparent that neglecting A, C, as in Refs. ||195| , |198|| , is a poor approximation in the 
presence of DCSB; i.e., on the domain where B{j)^^) ^ m, which must at least introduce quantitative 
errors. 



In the chiral limit: m = 0, Eq. ( |3.2.10|) reduces to a quadratic equation for B{p^^,), which has two 



qualitatively distinct solutions. The "Nambu-Goldstone" solution, for which 



5(p^J = <; V''' '^(P-J<'r (32.12) 

( , otherwise 

^^"'^^ ~ \ i (1 + y^l + 2r/Vp2 J , otherwise, ""^ 

describes a phase in which: 1) chiral symmetry is dynamically broken, because one has a nonzero quark 
mass-function, B{pi^^), in the absence of a current-quark mass; and 2) the dressed-quarks are confined, 
because the dressed-quark 2-point Schwinger function violates reflection positivity. The alternative 
"Wigner" solution, for which 

Bip^,) ^ , Cip^,) = i (1 + vTTv7<) , (3.2.14) 
characterises a phase in which chiral symmetry is not broken and the dressed-quarks are not confined. 
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To explore the phase transition one can consider the relative stability of the confined and deconfined 
phases, which is measured by their difference in pressure. Using Eq. ( p.l.7|) alone, because the gluon 
contributions cancel, that difference is 



\p\'A' + ulC' + B 
which in the particular case considered here yields 




(3.2.15) 

+ \p\H^A-(rA)+^U(rc-ac)] , (3.2.16) 



In 



(3.2.17) 



with: f = T/r], fi = ^m/t]; uf^^^ < ^ + p,^^ = of^^^ - ujf, pi = {y,ui + ip). B{fi,T) defines a 

(3.2.18) 



(/i, T)-dependent "bag constant" ||199|| and in this model 

i3(0,0) = (0.102 r^)^ = (0.109 GeV)^ 



which can be compared with the value ~ (0.145 GeV)^ commonly used in bag-like models of hadrons. 
The positive value indicates that the confining phase of the model with DCSB is stable at /i = = T. 
The phase boundary in the (yU, T)-plane is defined by the line i3(/i, T) = , which is evident in Fig. ^A . 
In this model the transition is first order except at = 0, which is clear from the nonzero derivative at 
the phase boundary, and the deconfinement and chiral symmetry restoring transitions are coincident, 
which is a consequence of the character of the solutions of the gap equation, Eqs. ( |3.2.12| - ^.2.14|) . The 
T = transitions are also coincident in the model obtained with T>b, Eq. ( p.2.3| ) [ |197|| . The chiral 
order parameters increase with increasing /i, as they do in all models that preserve the momentum 
dependence of the dressed-quark self energies; e.g, Refs. ||195| , |197| , |198| |, typically reaching a value < 20% 
larger than that at = = T. This is an anticipated result of confinement: as long as /i < /ic, each 
additional quark must be paired with an antiquark thereby increasing the density of condensate pairs. 
It cannot be observed in models that omit the momentum dependence of the dressed-quark self energy, 
such as those with an instantaneous interaction. 



The criterion of maximal pressure was also employed in Ref. | 198 ], wherein the model effective interaction 
saturates in the infrared and exhibits one-loop improved ultraviolet behaviour analogous to that in 
Eq. ( ^.2.351 ). The calculated = critical temperature agrees with the estimate in Table ^TT] and the 
T = critical chemical potential: fic = 422 MeV, is just 13% larger than that estimated in Ref. ||197|| . 
A similar interaction was employed in Ref. [|195|| , wherein (/i = 0,Tc) agrees with Table |3.1| . However, 
the estimate of {fic, T = 0) is 60% larger than in Refs. ||1971 , |198|| and is likely too high. The usual range 
is /ic e [300,400] MeV; e.g., see also Refs. i200| , M - 



As clear in Fig. |3.4| , the chiral phase boundary traces out a curve in the (yU, T)-plane. In Ref. | ]198| |, 
beginning at {fic,T = 0), this curve describes a line of first order critical points that ends with a 
tricritical point at (/itc ~ 0.5 fic, ^tc ~ 0.8 Tc). For T > Ttc the transition is second order. In comparison 
with the Pa model, this is a shift of the tricritical point away from the /i = axis. We expect such 
a relocation to be a feature of all models with a better treatment of the ultraviolet behaviour of the 
effective interaction; e.g., there are preliminary indications that in the separable model of Ref. ||194 
a tricritical point is located at (/itc ~ 0.4/ic, ~ 0.9 T^). The existence of a tricritical point has 



observable consequences ||202|| . However, if it occurs at a value of /i as large as that estimated in these 
models it may be difficult to detect using the relativistic heavy ion collider [RHIC] , which will focus on 
the domain of low baryon number density. 
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Fig. 3.4. T) in Eq. ( p.2.17| ): T) > marks tlie domain: T>, of confinement and 
DCSB, with the line B{fi,T) = demarcating the phase boundary. At /i = 0, the transition 
is second order with mean field critical exponents and = 169 MeV, using rj = 1.06 GeV. 
For all nonzero fi the transition is first order with a tricritical point at [fi = 0,Tc). At 
T = 0, /ic = 300 MeV. The value is 25% larger; i.e., /ic = 375 MeV, in the model defined 
by Vb, which improves the ultraviolet behaviour of the effective interaction [ |1971| , and in a 
two-flavour, ultrarelativistic Fermi gas at this chemical potential the baryon number density 
is 2.9po, Po = 0.16 fm-3. (Adapted from Ref. |T76^ .) 



With B{fi, T) > 0, then, in the exploration of hadronic observables using the rainbow-ladder truncation, 
Pi;ng(P' describes the free energy of the "vacuum" or "ground state," relative to which all excitations 
are measured V (/i, T) G T>] i.e., in the domain of confinement and DCSB. Therefore Psng is not directly 
observable and does not contribute to the thermodynamic pressure. Nevertheless, it does evolve with 
changes in (/i,T), as evident in Fig. and this evolution simply reflects the (/x, T)-dependence of 
the dressed-quark self energies, upon which it depends. [NB. If it did not evolve there could never 
be a phase transition.] The changes in the dressed-quark self energy represent the (/i, T)-dependent 
rearrangement of the vacuum and are indirectly observable; e.g., in the behaviour of hadron masses and 
decay constants. Hence, on the domain T>, they are evident in changes of the only true contributions to 
the pressure; i.e., those of the colour singlet hadronic correlations. At the phase boundary, denoted by 
dV: pY;j^^{fi,T)\gjy = pY;^{fj,,T)\gjj. This is the last vacuum reference point: hereafter dressed-quarks 
and -gluons are the active degrees of freedom. 

The dressed-quarks' contribution to the pressure is 

P,{fi,T) = 9{V) {pj:^{fi,T) - pj:^{fi,T)\g^} , (3.2.19) 

where 0{I)) is a step function, equal to one for {T,fi) G V. This pressure can be calculated numerically 
and is depicted in Fig. 3^ for model T>\. As evident in Eq. ( p. 2. 14 ), nonperturbative, momentum- 
dependent modifications of the dressed-quark's vector self energy persist into the deconfined domain 
and this feature retards the approach of the pressure to the ultrarelativistic limit: 

1 



Pur(p,^) 



127r2 



(3.2.20) 
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Fig. 3.5. Pq{n,T) in Eq. (|3.2.19| ) calculated using the Va and normalised by the pressure 
of an ultrarelativistic Fermi gas, Eq. (|3.2.20|) . This ratio approaches one slowly; e.g., at [fi = 
0, T ~ 2 Tc) or (yU ~ 3 /ic, T = 0), Pq{fi, T) is only one-half of the free particle pressure. The 
retardation is a consequence of the persistence into the deconfined domain of nonperturbative 
effects in the dressed-quark's vector self energy. Note the "mirroring" of the T-dependence 



in the evolution with /x, which we anticipate is a general feature. (Adapted from Ref. ||176 |. 



A qualitatively similar result is observed in numerical simulations of T 7^ lattice-QCD [ |175| ]. Such a 
correspondence is impossible in models where the dressed-quark vector self energy is omitted. 



The pressure in Eq. ( p.2.19|) is incomplete. It is zero in the confined domain where, as remarked 
above, the only contribution is that of colour singlet hadronic correlations. Those contributions can be 



estimated; e.g., using the hadronisation techniques of Ref. ||203|| , but that has not yet been done. In the 
deconfined domain, one must add the contribution of the dressed-gluons. However, hitherto the model 
studies have employed "frozen" gluons; i.e., the structure of the effective interaction does not evolve 
with (/i, T). T 7^ lattice-QCD simulations suggest that to be a good approximation, at least until very 
near the phase boundary [|175|| , and Ref. [p.87|| indicates that changes in the interaction in this vicinity 
have little effect on the properties of the transitions. That is good because improving on the frozen 
approximation is difficult, requiring information about the effective interaction in a new domain. A first 
step is to allow the mass-scale characterising the infrared behaviour of the interaction to drop suddenly 
and significantly near the phase boundary, thereby mocking-up the T-dependence of the string tension. 
Such a study is presently lacking. 



Not withstanding these weaknesses, the pressure depicted in Fig. |3.5| still provides a rudimentary model 
for the dressed-quark pressure that exhibits significant new qualitative features; e.g., the persistence 
into the deconfined domain of nonperturbative effects in the vector self energy yields an EOS much 
softer than that of a bag model. In this capacity it has been used to explore the structure and stability 
of nonstrange quark stars | 204 |. References [176 ,197, 205 1 indicate that a transition to quark matter 



should occur at approximately three-times nuclear matter density, and the comparison between results 
obtained with the EOS depicted in Fig. |3]^ and a bag model EOS shows; e.g., that the prediction for 
the maximum stable radius of a pure quark star is model-independent: 8-10 km, also in agreement 
with estimates based on a colour dielectric model ||206|| . The maximum mass, however, is very sensitive 
to the EOS. The quark core described by Eq (|3.2.19| ) must be surrounded by an hadronic crust and 
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that will also affect the properties of the star. The problem rapidly becomes involved, and determining 
the composition of dense astrophysical objects and those signals that point to the realisation of quark 
matter in their core is an important focus of contemporary research; e.g., Refs. ||194| , p06| , p07| , |208| . 



Superfluidity in Quark Matter In Sec. |2]^ we described aspects of diquark correlations and the man- 
ner in which they help to understand baryon properties. Another aspect follows from the meson-diquark 
auxiliary-field effective action ||147| , |148|| : the steepest descent approximation to the vacuum pressure 
suggests the possibility of diquark condensation; i.e., quark-quark Cooper pairing, and that was first 
explored using a simple version of the Nambu-Jona-Lasinio model ||209| |. A nonzero chemical potential 
promotes Cooper pairing in fermion systems and the possibility that it is exhibited in quark matter was 
considered using the rainbow-ladder truncation of the gap equation |210|. Interest in that possibility 



has been renewed A quark-quark Cooper pair is a composite boson with both electric and colour 



charge, and hence superfluidity in quark matter entails superconductivity and colour superconductivity. 
However, the last feature makes it difficult to identify an order parameter that can characterise a tran- 
sition to the superfiuid phase: the Cooper pair is gauge dependent and an order parameter is ideally 
describable by a gauge-invariant operator. 

As we saw above, at yU = = T QCD exhibits a nonzero quark- ant iquark condensate but it is under- 
mined by increasing fi and T, and there is a domain of the (/i, T)-plane, evident in Fig. |3.4| , for which 
(qq) = 0. Increasing T also opposes Cooper pairing. However, since increasing fi promotes it, there 
may be a (large-/i,low-T)-subdomain in which quark matter exists in a superfiuid phase. While that 
domain may not be accessible at RHIC, such conditions might exist in the core of dense astrophysical 
objects IplOII , which could undergo a transition to superfiuid quark matter as they cool. Unambiguous 



signals of such a superfiuid phase are actively being sought; e.g, Refs. [|194|j208| 



We take Ref. ||151|| as our exemplar, in which it was observed that a direct means of determining whether 



a SUc(A^) gauge theory supports O"*" diquark condensation is to study the gap equation satisfied by 

Here T = 0, for illustrative simplicity and because temperature can only act to destabilise a Cooper 
pair, and, with cj^^j = p4 + ifi, 

D{p, /i) = Z7 ■ pA{p\ ul^) + B{p\ ool^) + a;^ C{p\ ool^) ; (3.2.22) 

i.e., the inverse of the dressed-quark 2-point function, Eq. (|3.1.11|) . In Eq. ( p. 2. 211) , {A^, i = 1 . . .n^, 
= Nc{Nc — l)/2} are the antisymmetric generators of SUc(Ai'c) and C is the charge conjugation 
matrix, Eq. ( p.3.8|) . Using such a gap equation to study superfluidity makes unnecessary a truncated 
bosonisation, which in all but the simplest models is a procedure difficult to improve systematically. 

In addition to the usual colour, Dirac and isospin indices carried by the elements of D{p, fi), the explicit 
matrix structure in Eq. (|3.2.21| ) exhibits a quark bispinor index and is made with reference to 



Q{x) := ^^^^ t%Cct ) ' ^^"^^ ( ^^"^^ ^^"^^ '^'^''^ ) ' ^^-^-^^^ 

where {rj : i = 1,2,3} are Pauli matrices, Eq. ( p.l.7|) , that act on the isospin index. (Only two-flavour 
theories are considered in this exemplar. Additional possibilities open in three-flavour theories pi2[ .) 



As we have seen, a nonzero quark condensate: (qq) ^ 0, is represented in the solution of the gap equation 
by -B(p^, cu^^]) ^ 0. The new but analogous feature here is that diquark condensation is characterised 
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by A*(p, n) ^ 0, for at least one i. That is clear if one considers the quark piece of the QCD Lagrangian 
density: L[q, q]. It is a scalar and hence L[q, g]* = L[q, q]. Therefore L[q, q] oc L[q, q] + L[q, qf, and that 
sum, and hence the action, can be re-expressed in terms of a diagonal matrix using the bispinor fields 
in Eq. (|3.2.23| ): Q dia:g[D,CD^C^Q. It is plain now that a dynamically-generated lower-left element in 
D{p,fi), the inverse of the dressed-bispinor propagator, corresponds to a qq [diquark] correlation. 

The bispinor gap equation can be written in the form 

D{p, /.) = Doip, /i) + ( i^^;^^' ^\ ^J^fP' \ , (3.2.24) 

where in the absence of a diquark source term 

£>o(p, ii) = {ii-p + m)Tl -fiT^, (3.2.25) 

with m the current-quark mass, and the additional Pauli matrices: {tq,^ = 0,1,2,3}, act on the 
bispinor indices. The structure of Sjj(p, /i) specifies the theory and, in practice, also the approximation 
or truncation of it. 

Two colour QCD [QC2D] provides an important and instructive example. In this case = Ar^ in 

Eq. ( |3.2.21| ), with ifc the generators of SUc{2), and it is useful to employ a modified bispinor 

with Q2 the obvious analogue of Q in Eq. ( |3.2.23| ). Now the Lagrangian's fermion-gauge-boson inter- 
action term is simply 

Q2ix) ^^gi.T^Q Q2ix) 4(x) (3.2.27) 

because SUc(2) is pseudoreal; i.e., (— t^c)* Tc = 'fc , and the fundamental and conjugate representations 
are equivalent. (That the interaction term takes this form is easily seen using g]* = L[g, q].) 

Using the pseudoreality of SUc(2) it can be shown that, for = and in the chiral limit, the general 



solution of the bispinor gap equation is |151 



D{p) = i-i -pAip") + V{-iv)M[p') , V(7r) = exp |i75 T^vr^j = V(-7r)-^ , (3.2.28) 

where 7r^=i'-'5 are arbitrary constants and {T^'^'^ = r|(g)r7, = r^®r^, = r^®r)?}, {T\ T^} = 
so that 

S{p) = p2^2(p2)^_yvi2(p2) •= -n-VOv{p)^V{i^)os{p). (3.2.29) 

[tt = (0, 0, 0, 0, — ivr) produces an inverse bispinor propagator with the simple form in Eq. (|3.2.21|) .] 
That the gap equation is satisfied for any constants vr^ signals a vacuum degeneracy: if the interaction 
supports a mass gap, then that gap describes a five-parameter continuum of degenerate condensates: 

(Q2V(7r)Q2)7^0, (3.2.30) 

and there are 5 associated Goldstone bosons: 3 pions, a diquark and an anti-diquark. In this con- 
struction, Eq. ( |3.2.28| ), one has a simple elucidation of a necessary consequence of the Pauli-Giirsey 
symmetry of QC2D. For m 7^ 0, the gap equation requires trj^Q [T*V] = , so that in this case only 
{Q2Q2) 7^ and the Goldstone bosons are now massive but remain degenerate. 
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For /i 7^ the general solution of the gap equation has the form 

and in the absence of a diquark condensate; i.e., for A = 0, 

[f/B(«),£)(p,/i)] =0, UB{a) :=e^°-^®^/, (3.2.32) 



which is a manifestation of baryon number conservation in QC2D. In this case the explicit form of the 
gap equation is complicated but its features and those of its solutions are easily illustrated using the 
model Pa introduced above; e.g., it becomes clear that a chemical potential promotes fermion-pair 
condensation, at the expense of fermion-antifermion pairs, and A G IR, V/x. The rainbow-truncation 
gap equation was solved using this model in Ref. ||151|| , and the relative stability of the quark- and 
diquark-condensed phases measured via the pressure difference 

Sp{n) := Ps[B=o,A] - PmA=o] (/^) > (3.2.33) 

where the pressure here is obtained as an obvious generalisation of Eq. ( p.l.Tf ). Sp{fi) can be expressed 
in terms of /i-dependent bag constants 

Bsifi) :=Ps[i?,A=o](At) -PE[B=o,A=o](/^) > ^A(/i) := J5e[b=o,a] -Ps[i?=o,A=o](/i) , (3.2.34) 

which measure the stability of a quark- or diquark-condensed vacuum relative to that with chiral 
symmetry realised in the Wigner-Weyl mode. [NB. Improving on rainbow-ladder truncation may yield 
quantitative changes of < 20% in the exemplary results that follow, however, the pseudoreality of QC2D 
and the equal dimension of the colour and bispinor spaces, which underly the theory's Pauli-Giirsey 
symmetry, ensure that the entire discussion remains qualitatively unchanged.] 

The (m, /i) = degeneracy of the quark and diquark condensates, Eq. ( p.2.30|) , is manifest in 

Bb{0) = Ba{0) = (0.13 mj=i)^ , (3.2.35) 

where m j=i is the (m, fi) = mass of the model's vector meson, calculated in rainbow-ladder truncation. 
Increasing /i at m = and excluding diquark condensation, chiral symmetry is restored at 

/xf/=° = 0.40 mj=i, (3.2.36) 

where Bb{ij) = 0. However, V/i > 0: 5p{jji) > and Bi^{ji) > , with 

^A(/i£^=°) = (0.20 mj=i)^ > Sa(0). (3.2.37) 

Therefore the vacuum is unstable with respect to diquark condensation for all > 0, and this is always 
dynamically preferred over quark condensation. 

(5 = 0, A ^ 0) in Eq. ( ^X3lD corresponds to tt = (0, 0, 0, 0, ivr) in Eq. ( p.2.30|) ; i.e., ((52«75^qQ2) ^ 0. 
The usual chiral [S't/^(2)] transformations are realised via 

D{p, /i) V{n) D{p, fi) V{n) , 1^(7?) := e'^'^'^ , vf = {n\ vr^, vr^) , (3.2.38) 

and therefore, since the anticommutator {T,T^'^} = 0, a diquark condensate does not break chiral 
symmetry. However, (i? = 0, A 7^ 0) does yield a dressed-bispinor propagator that violates reflection 
positivity and hence the model exhibits confinement to arbitrarily large densities. [NB. Reference [151 
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employs a frozen gluon approximation. In a more realistic analysis, the /i-dependence of t], the mass- 
scale characterising the model, would be significant for fx ~ /^fc^~°, and ^ as /i ^ oo, which would 
ensure deconfinement at large-/!.] Finally, although the yU 7^ theory is invariant under 

Q2 ^ UB{a) Q2 , Q2 ^ Q2 Uei-a) , (3.2.39) 

which is associated with baryon number conservation, the diquark condensate breaks this symmetry: 

(Q2Z75r^Q2) ^ cos(2a) (Q2^75r^Q2) - sin(2a) (Q2Z75r^Q2) • (3.2.40) 

Hence, for (m = 0,yU 7^ 0), one Goldstone mode remains. (These symmetry breaking patterns and the 
concomitant numbers of Goldstone modes in QC2D are also described in Ref. 



213| 



For m 7^ and small values of /i, the gap equation only admits a solution with A = 0; i.e., diquark 
condensation is blocked and this is because the current-quark mass is a source of quark condensation. 
However, with increasing /i a diquark condensate is generated and the model exhibits the following 
minimum chemical potentials for diquark condensation: 



0.013 mj=i ^ = 0.051 mj= 



1 5 



m 



0.13 mj=i ^ = 0.092 mj= 



(3.2.41) 



This retardation of diquark condensation by a nonzero current-quark mass can also be seen in Ref. [214 . 



The exploration of superfluidity in true QCD encounters two differences: the dimension of the colour 
space is greater than that of the bispinor space and the fundamental and conjugate representations of 
the gauge group are not equivalent. The latter is of obvious importance because it entails that the 
quark-quark and quark- ant iquark scattering matrices are qualitatively different. As we saw above in 
connection with Eq. ( p.3.42|) , these differences ensure that colour singlet meson bound states exist but 
[necessarily coloured] diquark bound states do not. 



n' 



3 in QCD and hence in canvassing superfluidity it is necessary to choose a direction for the 



condensate in colour space; e.g., A*A\ = A in Eq. (|3.2.21|) , so that 



where P\ 






A(p,/i)75A2 


-A(p, -/i)75A^ 




-p,fi)C^Pll+CD^ip,fi)C^P^ 



p± + P\\ 



(3.2.42) 



diag(l,l,l), and D\\, D± are deflned via obvious generalisations of 



Eqs. ( p.2.21| ), ( p.2.22| ). (NB. It is this selection of a direction in colour space that opens the possibility 
for colour-flavour locked diquark condensation in a theory with three effect ively-massless quarks; i.e., 
current-quark masses <^ fi [ pi5|[ .) In Eq. ( |3.2.42| ) the evident, demarcated block structure makes explicit 
the bispinor index. Here each block is a 3 x 3 colour matrix and the subscripts: [[, -L, indicate whether 
or not the subspace is accessible via A2. The bispinors associated with this representation are given in 
Eqs. ( |3.2.23| ) and in this case the Lagrangian's quark-gluon interaction term is 











-^7m(A")* , 



(3.2.43) 



It is straightforward to derive the gap equation at arbitrary order in the truncation scheme of Ref. 108 
and it is important to note that because 

D||(p,/i)P|,+Z}^(p,/i)Px = A°{fZ}|,(p,/i) + iDx(p,/i)} + ^A«{Z}||(p,/i)-Dx(p,/i)} (3.2.44) 
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Fig. 3.6. Illustration of the dressed-ladder vertex-corrected gap equation, which is the 
next-to- leading-order in the truncation scheme of Ref. ||108|| . [The ladder truncation is ob- 
tained by dropping the last diagram.] Each bispinor quark-gluon vertex is bare, given by 
Eq. ( p.2.43| ), but the shaded circles mark quark and gluon 2-point functions that are dressed. 
The corresponding truncation in the relevant Bethe-Salpeter equations ensures the absence 
of diquark bound states in the strong interaction spectrum. (Adapted from Ref. ||151| .) 



the interaction: r^5(p, /i)r", necessarily couples the ||- and ^-components. That interplay is discarded 
in models that ignore the vector self energy of quarks, which, as we have repeatedly seen, is a qualita- 
tively important feature of QCD iT6^ , [T76igT^ , gTl] . 



Reference ||151[| explored the possibility of diquark condensation in QCD via the gap equation using the 
model defined by Pa in both the rainbow and a vertex-corrected truncation. The latter is illustrated 
in Fig. |3.6| . For = the rainbow-ladder truncation yields 

2 --2 -1^2^ (gg)0 = (0.11r/)3, ^^(/i = 0) = (0.10 r/)^ (3.2.45) 



m , = m 



P 2 



and momentum-dependent vector self energies, Eq. (|3.2.12|) , which lead to an interaction between the 



and ^-components of D that blocks diquark condensation [pi8||. This is in spite of the fact that 



A'^A (— A*^)* = iA^A*^, which entails that the rainbow-truncation quark-quark scattering kernel is purely 
attractive and strong enough to produce diquark bound states [ |150| . [Remember that in the colour 
singlet meson channel the rainbow-ladder truncation gives the colour coupling A^A"^.] For /i 7^ and 
in the absence of diquark condensation, we saw in connection with Fig. (|3.4|) that the model exhibits 



coincident, first order chiral symmetry restoring and deconfining transitions at 

f-^ttrlo. = 0-28 V = 0.3 GeV . (3.2.46) 



For /X 7^ 0, however, the rainbow-truncation gap equation admits a solution with A{p, fi) ^ and 
B{p,fi) = 0. Sp{fi) in Eq. ( |3.2.33| ) again determines whether the stable ground state is the quark- 
condensed or superfiuid phase. With increasing /i, Bb{ij) decreases, very slowly at first, and B/^{ii) 
increases rapidly from zero. That evolution continues until 

A'-c, rainbow — 0.25?] = 0.89 /i^^ ^j^j^jj^^ , (3.2.47) 

where B/^^ii) becomes greater-than Bb{ij)- This signals a first order transition to the superfiuid ground 
state and at the boundary 

(Q275r2A2Q) s^o.A =(0.65)=^(QQ)^=o. (3.2.48) 

^ "c, rainbow 

The chemical potential at which the switch to the superfiuid ground state occurs is consistent with 



other estimates made using models comparable to the rainbow-truncation class 



20C, 


201 


214 


21C 



the magnitude of the gap at this point pnU| , pCTTiPT3| , pD|l . 
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A question that now arises is: how sensitive is this phenomenon to the nature of the quark-quark 
interaction? As we discussed in connection with Eq. ( |2.3.42|) , the inhomogeneous ladder BSE exhibits 
particle-hke singularities in the O"*" diquark channels and such states do not exist in the strong interaction 
spectrum. Does diquark condensation persist when a truncation of the gap equation is employed that 
does not correspond to a BSE whose solutions exhibit diquark bound states? The vertex corrected gap 



equation depicted in Fig. 2]6| is just such a truncation and it was also studied in Ref. |[L51|| 



In this case there is a A ^ solution even for /i = 0, and using X>a 

ml = (1.1) m^l^dder^ ^gg^ _ ^1 0)3 (QQ)™'^bow^ _ ^l j)4^minbow ^ (3.2.49) 

where the rainbow-ladder results are given in Eqs. ( |3.2.45| ), and 

{Qt^^r^X^Q) = {0A8f {QQ) , Ba = (0.42)^ Bb . (3.2.50) 

The last result shows, unsurprisingly, that the quark- condensed phase is favoured at /i = 0. Precluding 
diquark condensation, the solution of the vertex-corrected gap equation exhibits coincident, first order 
chiral sjnnmetry restoring and deconfinement transitions at 

/.f'^=° = 0.77/.f;t:bow (3.2.51) 



Admitting diquark condensation, however, the /i-dependence of the bag constants again shows there is 
a transition to the superfluid phase, here at 

;,f =0.^ = 0.63 /if '^=° , with (g^75r^A2Q)^=o.63Mf = (O-^l)' {QQ),=o ■ (3.2.52) 

Thus the material step of eliminating diquark bound states leads only to small quantitative changes in 
the quantities characterising the still extant superfluid phase. 

Solving the inhomogeneous BSE for the O"*" diquark vertex in the quark-condensed phase provides 
additional insight ||151|| . At = and zero total momentum: P = 0, the additional [confining] 



contributions to the quark-quark scattering kernel generate an enhancement in the magnitude of the 
scalar functions in the Bethe-Salpeter amplitude. However, as P^ evolves into the timelike region, 
the contributions become repulsive and block the formation of a diquark bound state. Conversely, 
increasing n at any given timelike-P^ yields an enhancement in the magnitude of the scalar functions, 
and as II ^ fj,f'^^^ that enhancement becomes large, which suggests the onset of an instability in 
the quark-condensed vacuum. This "robustness" of scalar diquark condensation is consistent with the 
observations in Ref. | p21[ . However, the studies described herein do not obviate the question of whether 
the diquark condensed phase is stable with-respect-to dinucleon condensation [p22|| , which requires 
further attention. [NB. As remarked above, the inclusion of temperature undermines a putative diquark 



condensate and existing studies |P00| , pi4|| suggest that it will disappear for T > 60-100 MeV. However, 
such temperatures are high relative to that anticipated inside dense astrophysical objects, which may 
therefore provide an environment for detecting quark matter superfluidity.] 

Finally, this discussion illustrates that, in some respects; such as the transition point and magnitude 
of the gap, the phase diagram of QC2D is quantitatively similar to that of QCD. That is a useful 
observation because the simplest superfluid order parameter is gauge invariant in QC2D, and the fermion 
determinant is real and positive, which makes tractable the exploration of superfluidity in QC2D using 



numerical simulations of the lattice theory 1 223 ]. Hence, the results of those studies may provide an 



additional, reliable guide to the nature of quark matter superfluidity. 
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4 Density, Temperature and Hadrons 



4.1 Masses and Widths 



Temperature Hitherto we have canvassed the bulk thermodynamic properties of QCD at nonzero 
(/i, T). However, the terrestrial formation of a QGP will be signalled by changes in the observable 
properties of those colour singlet mesons that reach detectors. In this connection, a primary feature of 
the QGP is chiral symmetry restoration and, since the properties of the pion (mass, decay constant, other 
vertex residues, etc.) are tied to the dynamical breaking of chiral symmetry, an elucidation of the T- 
dependence of these properties is important; particularly since a prodigious number of pions is produced 
in heavy ion collisions. Also important is understanding the T-dependence of the properties of the scalar 
analogues and chiral partners of the pion in the strong interaction spectrum. For example, should the 
mass of a putative light isoscalar-scalar meson ||117|J118| , |119| , |224|| fall below 2171^^, the strong decay into a 



two pion final state can no longer provide its dominant decay mode. In this case electroweak processes 



will be the only open decay channels below Tc and the state will appear as a narrow resonance |p02 
Analogous statements are true of isovector-scalar mesons. 



The effective interaction denoted above as T>c in Eq. ( |3.2.2|) , with uj = 1.2 m*, has been employed [|121 | 



in an exploration of the T-dependence of scalar and pseudoscalar meson properties. As we saw in 
Sec. mesons appear as simple poles in 3-point vertices and these vertices alone provide information 
about the persistence of correlations away from the bound state pole, which can be useful in studying 
the T-evolution of a system with deconfinement. For T ^ and two flavours, the ladder-truncation of 
the inhomogeneous BSE for the isovector O""*" vertex is 



r;.(p.,;Po;C) = ^4ir^75- r lg'D^,ip^^-q^^)^,Siq^^)Tliq^-Po;C)Siq-^)^., (4.1.1) 

Jl,q 



where q^^ = q^^ ±Pq/2, and with Pq = (P, 0) this is the equation for the zeroth Matsubara mode. [This 
is an extension of Eq. ( p.3.4|) . Remember, mji{() rpg(p^^; Pq^; Q is renormalisation point independent.] 



The solution of Eq. ( ^.1.1| ) has the form [hereafter the ^-dependence is often implicit] 

r;,(p^,;P)= (4.1.2) 
iEp,{p^^;P) + 7 ■ PFp,{p^^; P)+^-pp- PGj,(p^,; P) + ^^cu^ p- P G^,{p^^; P)] , 



where the T 7^ analogues of the a^^-like contributions in Eq. ( |2.3.9D are omitted because they play 
a negligible role at T = The breaking of 0(4) symmetry is responsible for making this T 7^ 



amplitude more complicated than its zero-temperature counterpart. For the higher Matsubara frequen- 
cies the form is still more complicated, with three additional terms. The scalar functions in Eq. ( |4.1.2| ) 
exhibit a simple pole at P^ + = so that 

rJC) 

r;s(P^.; P) = KiP^.;P) + regular, (4.1.3) 

— * 

where again regular means terms regular at this pole and r^(p^j.;P) is the bound state pion Bethe- 
Salpeter amplitude, canonically normalised via the obvious T 7^ extension of Eq. ( p.3.6|) . 



The pole position in Eq. ([4.1. 3|) determines the spatial "screening-mass" of the pion's zeroth Matsubara 
mode and its inverse describes the persistence length of that mode at equilibrium in the heat bath. There 
is a screening mass and amplitude for each mode, and each mode's amplitude is canonically normalised. 
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The full T 7^ bound state propagator can be calculated via any polarisation tensor that receives a 
contribution from the bound state but only once all the screening masses have been determined. (For 
the pion: the pseudoscalar or pseudovector polarisations will serve [0.) The propagator so obtained 
is defined only on a discrete set of points along what might be called the imaginary-energy axis and 
the "pole-mass;" i.e., the mass that yields the bound state energy pole for p ~ 0, is obtained only after 
an analytic continuation of the propagator onto the real-energy axis. [The loss of 0(4) invariance for 
T 7^ means that, in general, the pole mass and screening masses are unequal.] That continuation is 
not unique but an unambiguous result is obtained by requiring that it yield a function that is bounded 
at complex- infinity and analytic off the real axis |p25| . From this description it is nonetheless clear that 



the screening masses completely specify the properties of T 7^ bound states. Furthermore, the often 
used calculational expedient of replacing the meson Matsubara frequencies by a continuous variable: 
Qn ^ —iv-, and the identification of the energy scale thus obtained with a pole mass, is seen to be merely 
an artefice. However, since this prescription yields the correct result for free particle propagators, it 
might provide an illustrative guide. 



The residue in Eq. ( [4.1.3| ) is 



5'^ ir^ = Z,iv r ir^75xi(g.,; P) , (4.1.4) 



where XiriQuji', P) '■= S{q^^)TT,{q^^; P)S{q^^) is the unamputated Bethe-Salpeter wave function. [Ref. ||121 



employs the = 92 MeV normalisation, cf. Eq. ( |2.3.22| )] Substituting Eq. ( ^X3| ) into Eq. (|iXTD and 
equating pole residues yields the homogeneous pion BSE, which provides the simplest way to obtain 
the bound state amplitude. As already noted, the pion also appears as a pole in the axial-vector vertex 
and there the residue is the leptonic decay constant 

5'^PU = Z^tr r ir^757Xi(g.pi^)- (4-1.5) 

At T = in this Vq model the calculated chiral limit values are [remember: f.^ = 92 MeV normalisation] 
/° = 0.088 GeV, -(gg)°Gev2 = (0.235 GeV)^, rO(lGeV^) = (0.457 GeV)^. (4.1.6) 



The analogue of Eq. ( |4.1.1|) for the O"*""*" vertex is 

f-A 

n,q 

where a = 0, 1, 2, 3. However, as discussed in connection with Eq. ( p.3.23|) , the combination of rainbow 



r^(p^,;Po;C) = Z.^T^l - r y'D,,{p^,-q^,)^,S{q^;)Tl{q^-Pn^-OS{qZ^)^u, (4.1.7) 

Jl.a 



and ladder truncations is not certain to provide a reliable approximation in the scalar sector. Nev- 



ertheless, in the absence of an improved, phenomenologically efficacious kernel, Ref. ||121|| employed 
Eq. ( |4.1.7| ) in the expectation that it would provide some qualitatively reliable insight, an approach 
justified a posteriori. 

The solution of Eq. (^4.1.7|) has the form 

n(Pc.,;P)= (4.1.8) 

E,{p^^;P) + i7 ■ pGl{pu,^;P) + i^^uj^, G^ip^^; P) + i^ ■ Pp- P F,{p^^; P)] , 



r-1 



2 



where here the requirement that the neutral mesons be charge conjugation eigenstates shifts the p ■ P 
term cf. the O""*" amplitude in Eq. (^4.1.2| ). As already observed in connection with Eq. (|2.3.69| ), the 
scalar functions in Eq. ( [4.1.8|) exhibit a simple pole at P^ + = with residue 

r^r. = Z4tr rir"x^(gc.pP). (4.1.9) 

Jl,q 
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Fig. 4.1. Leading Dirac amplitude for the pseudoscalar (shaded symbols) and scalar (open 
symbols) vertices [E in Eqs. (^4.1.2| ), (|4.1.8| )] evaluated at (p = 0,a;o) and plotted as function 
of in the chiral limit; i.e., E{p = 0, ujq; P^), for a range of temperature values. The bound 



state poles are evident in each case. (Adapted from Ref. ||121 



However, since a V — A current cannot connect a O"*"^ state to the vacuum, the scalar meson does not 
appear as a pole in the vector vertex; i.e, 

6''^PU = Z^tT r ir"7x^(g.,;P) = 0. (4.1.10) 

J l,q 

The homogeneous equation for the scalar bound state amplitude is obtained, as usual, from Eq. ( |4.1.7| ) 
by equating pole residues. 

The T-dependence of the pole positions in the solution of the inhomogeneous BSEs is illustrated in 
Fig. [4.1| , from which it is evident that: 1) at the critical temperature, Tc = 152 MeV in Table |3.1| , 



one 



has degenerate, massless pseudoscalar and scalar bound states; and 2) the bound states persist above 
Tc, becoming increasingly massive with increasing T. These features are also observed in numerical 
simulations of lattice-QCD ||175|| . The bound state amplitudes are obtained from the homogeneous 



BSEs. Above T^, all but the leading Dirac amplitudes: E^^, E„, vanish and the surviving amplitudes 
are pointwise identical. These results indicate that the chiral partners are locally identical above 
they do not just have the same mass. 

The results are easily understood algebraically. The BSE is a set of coupled homogeneous equations 
for the Dirac amplitudes. Below each of the equations for the subleading Dirac amplitudes has an 
"inhomogeneity" whose magnitude is determined by Bq, the scalar piece of the quark self energy, which 
is dynamically generated in the chiral limit. Bq vanishes above eliminating the inhomogeneity and 
allowing a trivial, identically zero solution for each of these amplitudes. Additionally, with i?o = the 
kernels in the equations for the dominant pseudoscalar and scalar amplitudes are identical, and hence 
so are the solutions. It follows from these observations that the Goldberger-Treiman-like relation |^ 

/X(p^.;0) = i?o(p.J (4.1.11) 

is satisfied for all T only because both /° and Bq{p^^) are equivalent order parameters for chiral 
symmetry restoration. 
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T-dependence of the meson masses-squared and 

2 



Fig. 4.2. 

pole-residue matrix elements, Eqs. (|4.1.4|) , ( [4.1.9|) , along with the mass-squared: rnl^^, and 



residue: rgym, calculated in the chirally symmetric i?o = phase. For T < Tc 



which is a signal of the instability of the chirally symmetric phase at low T. For T > Tc, 
= rri^ = m'^^^. Right Panel: T-dependence of the masses and pion decay constant, 
Eq. (4.1.5). = within numerical error. (Adapted from Ref. [121|). 



0.2 



Using the calculated bound state amplitudes and dressed-quark propagators, the T-dependence of the 
matrix elements in Eqs. ( [4.1.4| ), ( |4.1.5| ), ( |4.1.9| ) follows. It is depicted in Fig. [4.2| for the chiral limit. 
Below Tc the scalar meson residue in the scalar vertex, To- in Eq. ( [4.1.9|) , is a little larger than the residue 



of the pseudoscalar meson in the pseudoscalar vertex, r^^ in Eq. ( |4.1.4| ). However, they are nonzero and 



equal above Tc, which is an algebraic consequence of -Bq = and the vanishing of the subleading Dirac 
amplitudes. As a bona fide order parameter for chiral symmetry restoration 



/°ocv^, t = (T/Tc-l)<0, (4.1.12) 



as anticipated from the discussion accompanying Table |3?T| . The same result was obtained in Ref. ||187| 
where it was observed too that 

1 



('"JJ)em 



cx/^ (4.1.13) 



i.e., the pion charge radius diverges at Tc in the chiral limit. This is plausible but Ref. [121] did not 



attempt its verification. Reference |121| did confirm the behaviour of the ratio of pole residues observed 



in Ref. 186 



r°(C) 1 

"^^V^, t<0. (4.1.14) 



The concomitant results: /° = and r° 7^ for t > 0, demonstrate that the pion disappears as a pole 
in the axial-vector vertex ||181|| but persists as a pole in the pseudoscalar vertex. Evident also in Fig. ^]2| 



is that 

m°(xv^, t<0: (4.1.15) 

follows a linear trajectory for t < 0. Such behaviour in the isoscalar-scalar channel might be 
anticipated because this channel has vacuum quantum numbers and hence the bound state is a strong 



interaction analogue of the electroweak Higgs boson ||119|| . 



m^y^ in Fig. [4.2| is the mass obtained when the chirally symmetric solution of the quark DSE is used 



in the BSE. [Bq = is always a solution in the chiral limit.] For t > 0, rn^yjj^iT) is the unique meson 
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Fig. 4.3. T-dependence of the meson masses for large T with m = 0, see Fig. [4.2| . m^r = mg- 
for T > Tc and m/(27rT) — > 1~. This behaviour persists with m 7^ 0, as illustrated in 
Refs. [|194|,p26||, and observed in lattice-QCD simulations ||175||.) (Adapted from Ref. ||121||.) 



mass-squared trajectory. However, for t < 0, rrigy^ < 0; i.e., the solution of the BSE in the Wigner- 
Weyl phase exhibits a tachyonic solution, [cf. The Nambu-Goldstone phase masses: > = 0.] By 
analogy with the cr-model, this tachyonic mass indicates the instability of the Wigner-Weyl phase below 
Tc and translates into the statement that the pressure is not maximal in this phase. Figure |4.3| depicts 
the evolution of this [common] meson mass at large T. As expected in a gas of weakly interacting 
quarks and gluons 

where ujq = ttT is a quark's zeroth Matsubara frequency and "screening mass." Similar behaviour is 
observed for the p-meson mass in Ref. ||226|| and can be demonstrated algebraically using the V\ model 



for the effective interaction |P27 



The results described hitherto were all calculated in the chiral limit. The extension to m 7^ is 
straightforward although calculations with the renormalisation group improved rainbow-ladder trunca- 



tion become more time consuming. That is why simpler models, such as employed in Refs. ||194| , p2(: 
can be useful. For m 7^ 0, chiral symmetry restoration with increasing T is exhibited as a crossover 
rather than a phase transition. The solutions of the inhomogeneous BSEs again exhibit a pole for all T, 
with the bound state amplitudes obtained from the associated homogeneous equations. In this case the 
scalar and pseudoscalar bound states are locally identical for T > |Tc. Figure is the m 7^ analogue 



of Fig. [4.2| . An important result is that the axial- vector Ward-Takahashi identity is satisfied, both above 
and below the chiral transition temperature, which was demonstrated in Ref. ||121|| via Eq. (|2.3.21|) : the 



two sides remain equal VT. The Gell-Mann-Oakes-Renner formula, however, which involves r^{C), fails 



for t > —0.1, as observed too in the separable model study of Ref. ||228 . 



As one can anticipate from Sec. |2.3| , the calculated a and n bound state amplitudes and dressed-quark 
propagator also make possible a study of two-body decays. For example, the impulse approximation to 
the isoscalar-scalar-TTTT coupling is described by the matrix element 

W := (vr(pi)7r(p2)|a(p)) (4.1.17) 
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Fig. 4.4. T-dependence of the meson masses and pole-residue matrix elements, 
Eqs. ( [4.1.4| , |CT75| , [4.1.9D , for m 7^ 0. That the transition has become a crossover is evident 
in the behaviour of fj^. The meson masses become indistinguishable at T ~ 1.2 Tc, a little 
before the local equivalence is manifest, which is unsurprising given that the mass is an inte- 
grated quantity. The small difference between and r^r below Tc is again evident and they 
assume a common value at the same temperature as the masses. (Adapted from Ref. ||121||.) 



0.2 



.A 

2A^ctrD / T„{kuji]p) Su{k++) iT^{ko+-, -pi) Su{k+^) iT^^k^o] -P2) Su{k. 



ka(3 = kuji + {a/2)pi + {/3/2)p2, in terms of which the width is 



2 JaTTTV 



1 — 4m^ / 
16 vrmo- 



(4.1.18) 



The coupling and width obtained from Eq. ( ^.1.17| ) are depicted in Fig. |4.5| , which indicates that both 
vanish at in the chiral limit. Again this can be traced to Bq — * 0. For m 7^ 0, the coupling 
reflects the crossover but that is not observable because the width vanishes just below where the 
isoscalar-scalar meson mass falls below 2m^ and the phase space factor vanishes. [See the right panel of 
Fig. [4.4|.] The evolution 2m.n- may, however, be observable via an enhancement in the tttt 
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cross-section [ p29|| 



Additionally, the particular properties of the tt^ — > 77 decay, which is mediated by the "triangle 
anomaly," [see the paragraph after Eq. ( |2.3.41| ) on page Q make interesting the behaviour of this 
process at T 7^ 0. At T = 0, the anomalous contribution to the divergence of the axial-vector vertex is 
saturated by the pseudoscalar piece, Et,-, of the pion Bethe-Salpeter amplitude [p2| 



T^u{ki, k2) = tr / ^(gi) ^^T-^iE^{q] -P) S{q2) iQe^M'i^ ^12) 'S'(gi2) iQe^v{qi2, qi) , (4.1.19) 

Jl,q 



(fci,0), k2 = (^2,0), P = ki + k2, qi = q^ji - q2 



where here, just to be specific, k\ = 

Q = ki^li + ^2), qi2 = quji — ki + k2- Equation ( |4.1.19| ) involves the dressed-quark-photon vertex: 
which also appeared in the calculation of F^^Q"^) 
reliable numerical solutions of the T 



described in Sec. 2.3 



As we saw. 



+ k2 
V 

quantitatively 
However, this 



A" 



vector vertex equation are now available 
anomalous coupling is insensitive to details and an accurate result requires only that the dressed vertex 
satisfy the vector Ward-Takahashi identity. For T = with real photons, Eq. ( 4.1.19| ) is expressed in 
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Fig. 4.5. T-dependence of the isoscalar-scalar-vrvr coupling and width, both in the chiral 
hmit and for reahstic hght current-quark masses. The phase space factor (1 — (2m^/mo-)^)^^^ 
is 6{Tc — T) in the chiral limit but nontrivial for rh ^ 0, vanishing at T 0.98 Tc; i.e., this 
decay channel closes at a temperature just 2% less-than T^. (Adapted from Ref. ||121|| .) 



terms of one scalar function: 



-/lupcr 



hp T(0) . 



(4.1.20) 



Now, as long as satisfies the Ward-Takahashi identity, Eq. ( p.3.26|) , one finds algebraically in the 
chiral limit 

/°r(0):=(7.o^^ = l/2, (4.1.21) 

completely independent of model details P2| , |130| , P^^ , as required since at T = the anomalies are a 
feature of the global aspects of DCSB ||143|| . This value reproduces the experimental width. [Remember, 
the normalisation in this subsection yields /^r = 92 MeV.] 



The T 7^ calculation requires only a valid extension of Eq. ( |2.3.31| ) and one such is 

«f (g,.,^ , g,., J = J:A{qt^^,qtjil + {qi + q2) [^iG{q^,^,q^^^) + AB{ql^^,qlJ], 
ir4(g^,^,ga.,J = ^ciql^^,ql^Jili+ {uJi, +uJi^) [^iG{q^^^,q^^^) + AB{ql^^,ql^^)], 



G{q^,^,q^,^) = 7 ■ (gi + q2)AA{ql^^ , qlj + 74(^«i + uJi^)Ac{q^^^ , q^^^ 



(4.1.22) 
(4.1.23) 
(4.1.24) 



It is a particular case of the Ansatz in Ref. | J23U| | and satisfies the T 7^ vector Ward-Takahashi identity 

{qu^i, -qu.i^)pi'^i.{qui^,qu.i^) = 5"^(9^,J - S'\q^^^). (4.1.25) 

For T 7^ the tensor structure of Eq. ( [1.1. 201) survives to the extent that, with ki, k2 as defined, it 
ensures one of the photons is longitudinal (a plasmon) and the other transverse, with 



fiA^ki, k2) 



TT 



(fcixfc2).r(o) 



(4.1.26) 



The T-dependence of the anomalous coupling follows from that of T(0), and it and the T-dependence 
of the width are depicted in Fig. In the chiral limit the interesting quantity is: T(0) = fi'°o^^//°, 
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71-0 _^ width: 

agj„m^|T(0)p/(167r^), which is identically zero in the chiral limit because rrin = 0. 



Fig 

(Adapted from Ref. |T2l .^ 



and obvious in the figure is that it vanishes at T^. (It vanishes with a mean field critical exponent ||121|| .) 
Thus, in the chiral limit, the coupling to the dominant decay channel closes for both charged anc? neutral 
pions. These features were anticipated in Ref. | |231|| . Further, the calculated T(0) is monotonically 
decreasing with T, supporting the perturbative 0(T^//^) analysis in Ref. ||232|| . For m ^ both 
the coupling: QnO-yy/fn, and the width exhibit the crossover, with a slight enhancement in the width 
as T ^ Tc due to the increase in m.,^. This is similar to the results of Ref. |P33|| , although the T- 
dependence depicted here is much weaker because the pion mass approaches twice the T 7^ free-quark 
screening- mass from below, never reaching it, Eq. ([4.1.16|) ; i.e., the continuum threshold is not crossed. 

The T-dependence of meson properties illustrated here is robust: it agrees with the results obtained in 
lattice simulations when there is an overlap, and with the results of other models. The local equivalence 
exhibited by isovector chiral partners above Tc might be expected as a general feature of the bound state 
spectrum in the Wigner-Weyl phase. However, an explicit demonstration is numerically challenging; 
e.g., in the p-ai complex the bound state amplitudes have eight independent functions even at T = 0, 
compared with the four in the pseudoscalar and scalar amphtudes at T 7^ 0. 

Adding a third light quark introduces one qualitatively new aspect: the rj-r]' system and the restoration 
of Ua{^) symmetry, which can affect the order of the chiral transition ||183|| . As already observed, it is 
necessary to move beyond the rainbow-ladder truncation before that can be addressed using the DSEs. 
The question has been explored in lattice simulations but the results are not currently conclusive: 
the mass splittings used to characterise the symmetry breaking might become smaller near ||234|| 
but strong, topological arguments can nevertheless be made in favour of the non-restoration of Ua{1) 
symmetry ||235|| . Much remains to be done and improved models can be useful. 



Chemical Potential The relation between chemical potential and baryon number density can only be 



determined after the EOS is known; i.e., via Eq. ( |3.1.5|) . As described in connection with Eq. ( 3.2. 19| ) 
on page ^ the dressed- quarks and -gluons contribute nothing to the EOS in the confined domain, even 
though they dominate it in the QGP, and the only true contributions to the pressure in the confined 
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domain are those of colour singlet bound states. This physical requirement is overlooked in many 
model explorations of the density dependence of meson properties. For example, in applications of the 
Nambu-Jona-Lasinio model the EOS for a free fermion gas is used to define a baryon number density. 
While this is the EOS for that non-confining model it is not a good model for QCD's EOS. 

The application of DSEs in calculating the /i-dependence of hadron properties is rudimentary. How- 
ever, even that is significant given the problem is currently inaccessible in lattice simulations. The 
model obtained with Pa in Eq- (|3.2.2| ) again provides a useful, algebraic exemplar. As described after 
Eq. ( |3.2.18| ), the QGP transition is first order at (/i, T = 0) and the chiral order parameters increase 
with increasing n when the dressed-quark self energy is momentum dependent. Mechanically, the latter 
is an obvious consequence of analyticity in the neighbourhood of the real axis: any function, 0(4) 
invariant at /i = = T, has the expansion 



2 ~,2n T^Or^fi 
,<^i.l — 



dy 



+ ... 



(4.1.27) 



y=Cjl=0 



and since \R[u 



— then, if IR[/(p^, u;^)] decreases with T^, the derivative is negative and 



IR[/(p^, cj^)] must increase with /i^. [Only the real-part is important because the imaginary-part of 
physical quantities vanishes after summing over the Matsubara frequencies. The derivative is zero in 
models with an instantaneous interaction.! 



Equation (|4.1.27| ) is exemplified in the behaviour of the pion's leptonic decay constant, which using the 
algebraic Va model is simply expressed: 



J n 



^ — ~ 



^max A 3 

TE ^(1 + 4 



Atuf 



8 

5 ^^l 



(4.1.28) 



where the notation is specified in connection with Eq. ( |3.2.17|) on page |5^. Its behaviour is depicted in 



Fig. [4.7|, as is that of the pion's mass, m^r is almost insensitive to changes in /i and only increases slowly 



with T, until T is very near the critical temperature, as already seen in Fig. f4.4| , which was calculated 
with the renormalisation-group-improved model. The insensitivity to fi mirrors that to T and is the 
result of compensating changes in r.„{() and f^; i.e., a consequence of the axial-vector Ward-Takahashi 
identity. All these features are also evident in Refs. ||195|jl98| | . 

A first step in exploring the properties of the p-meson is solving the the vector meson BSE. That too 

(4.1.29) 



takes a particularly simple form in the T)\ model [p27|| : 

n2 



There are two solutions: one longitudinal and one transverse with-respect-to P : 
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(4.1.30) 

|p|2 ^ I ^ J ■^P- 

where 6p^ labels the longitudinal and 6p^ the transverse solution. Substituting, one finds that for m = 

m'^p- = 5?7^, independent of /i and T. (4.1.31) 

Even for nonzero current-quark mass, mp_ changes by less than 1% as n and T are increased from zero 
toward their critical values. This insensitivity is just what one would expect for the transverse mode: 
remember, there is no constant mass shift in the transverse polarisation tensor for a gauge-boson. 
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Fig. 4.7. The pion's weak decay constant, Eq. (^4.1 .281) , as a function of /i for a range of 
values of T: increases with /x and decreases with T. [/^(/ic,T = 0) = 1.51/7^(0,0). In 
Ref. 1193, which employs model Pb, fAf^cO) = 1.25/^(0,0).] 
It falls slowly and uniformly with fi [mvr(/ic, 0) 



too 

(Adapted from Ref, 



The pion's mass is plotted 
0.95 m,r(0, 0)] but increases with T. 



For the longitudinal component one obtains in the chiral limit: 



m 



p+ 



\rf + 4(7r2r2 



/i 



(4.1.32) 



where the anticorrelation between the response of mp+ to T and its response to /i is plain. From 
Eq. (|4.1.32| ) and continuity at the second-order phase boundary it follows that 



{2TiTy 



(4.1.33) 



which is analogous to Eq. ( [4.1 .161) . Here, however, the limit is approached from above because mp+ 7^ 



in the chiral limit and increases with T. (NB. It is only because this model exhibits confinement that 
such a result is possible. Studying the p-meson in non-confining dressed-quark-based models requires 

p36i . 



However, 



that some means be employed to suppress or eliminate the p ^ qq threshold; e.g., Ref. 
that artefice is merely an indigent expression of confinement. For the transverse component of the p, 
rn?p- 1 {2ujqY ~^ 1" because of Eq. ( [4.1.31| ); e.g., Ref. [p26|| .) The (/i, T)-dependence of the p-meson mass 
is depicted in Fig. ^78|and; e.g., at T = 0, Mp+(/Xc) ~ 0.65 Mp+(/i = 0). As observed in the introduction. 



though, the connection between /i-dependence and baryon-density-dependence cannot be determined 
until the EOS is calculated. Without it one can only observe that, in a two-flavour free-quark gas, the 
T = critical chemical potential corresponds to 3pQ, see Fig. on page Therefore, at 1-2 po a 
mass reduction less-than this should be anticipated, plausibly no more than 25% ||227 |. 



The T-dependence described above was also observed in the confining model of Refs. ||194| , p26|| , wherein 
too it was found that, because of the T-dependence of m,r, f^p-, the dominant 2'7r-decay mode of the p_ 
meson mode is phase-space blocked for T/Tc > 1.2. [cf. The 2tt mode of the isoscalar-scalar discussed 
in connection with Eq. ( 4.1.18 ).] The T-dependence of the p-meson's dilepton decay width was also 
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Fig. 4.8. Mp_|_ and as a function of /i for T = 0,0.1. On the scale of this figure, 
is insensitive to this variation of T, which does not take approach near to {fi = 0,Tc). 

The current-quark mass is m = 0.011 rj, which for 7] = 1.06 GeV yields Mp+ = 770 MeV and 
= 140 MeV at T = = /i. (Adapted from Ref. p2^.) 



considered: it is suppressed by a factor of 0.9 in the vicinity of (/i = 0,Tc). However, its /i-dependence 
is yet to be explored. 

The anticorrelation, anticipated in Eq. ( fl.l.27| ), between the fi and T dependence of mass-dimension- 
two observables; such as X, /^r, mp+, etc., is apparent. It entails that, in these cases at least, there 



is a trajectory in the (/x, T)-plane along which the observables are constant |p37|| . It also means that 
observables calculated using the rainbow-ladder truncation do not exhibit a /^-scaling law of the type 
conjectured for baryon-number-density in Ref. [ ^38| . 



We anticipate that the mass of the Oi-meson [the chiral partner of the p] will decrease with increasing 
T so that it can evolve to meet the increasing p-mass, with mp+ = at the phase boundary. How- 
ever, since model Pa is defective in not supporting an axial-vector bound state, that remains to be 
verified. (It fails to support a scalar bound state too As remarked in Sec. ^]3|and Ref. ||115|| , much 



remains to be learnt about axial-vector mesons, in which connection confinement is an important ele- 
ment. Determining the /x-dependence of rria^ is particularly interesting given the fi cf. T anticorrelation 
exemplified above. 

In the studies described here, and also in lattice simulations, /i is an intensive thermodynamic param- 
eter whose presence modifies the propagation characteristics of dressed, confined particles, and this 
modification is transmitted to the observable hadrons they comprise. It is clear from the existence of 
a critical quark chemical potential, below which asymptotic quark states cannot be produced, that in 
the confined domain there is no simple proportionality between the quark chemical potential and the 
chemical potential associated with colour singlet baryons. In approaches based on elementary hadronic 
degrees of freedom; e.g., those reviewed in Ref. [p39|| , this consideration is bypassed. Colour singlet 
baryon density is introduced directly via the expedient of in-medium elementary meson and nucleon 
propagators, which are then employed in calculating the myriad nuclear-matter many-body loop in- 
tegrals that contribute to observable processes. The approach has a long history and yields a useful 
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phenomenology. However, it ignores the dressed-quark level effects described above and also questions, 
such as, just what is represented by an off-shell hadron propagator? [cf. The discussions associated 
with Eqs. ( |2.3.34| )-( [23:iq ), page |2|, and ( |2.3.69| )-( |2X79| ), page |8[] We judge that in understanding 
QCD at nonzero baryon density it is important to uncover the nature of the relationship between these 
approaches. 



4.2 Hadronic Signatures of the Quark-Gluon Plasma 

A QGP existed approximately one microsecond after the big bang, and primary goals of current gen- 
eration experiments at CERN and Brookhaven are the terrestrial recreation of the plasma and an 
elucidation of its properties. A number of signals for QGP formation at high temperature and low 
baryon number density have been suggested and here we briefly review three of them. They and others 



are discussed more extensively in Refs. ||240| , p41 



J/'^ Suppression. Reference |242] proposed "...that J/\l' suppression in nuclear collisions should 



provide an unambiguous signature of quark-gluon plasma formation." The reasoning is simple. A cc- 
pair produced in a hard parton-parton interaction will evolve into a J/^^-meson if the colour interaction 
is sufficiently strong to effect binding. That will always be the case unless the cc-pair is produced 
in a heat-bath of deconfined, colour- carrying excitations that [Debye-] screen the cc-attraction; i.e., 
unless the cc-pair is produced in a QGP. QGP formation occurs only for temperatures greater than 
some critical value. Hence the J /'^ production cross-section should evolve smoothly with controllable 
experimental parameters; such as, projectile and target mass numbers, and impact parameter, until a 
QGP is produced, when a dramatic suppression should follow. 

Following this suggestion, the J/\E' production cross-section has been systematically explored in rel- 
ativistic heavy ion collisions at CERN using the Super Proton Synchrotron [SpS]. As illustrated in 
Fig. all results are described by a simple scaling law: 

(^projectile X ^target)" , « = 0.92 ± 0.01 , (4.2.1) 

where A, B are the mass numbers, except those for Pb-Pb collisions. The scaling relation is easily 
understood as a consequence of nuclear absorption ||244|| and is termed "normal" J/ \E'- suppression. While 



the Pb-Pb data agree with this normal pattern for peripheral collisions; i.e., L < 8 fm, that is not the case 
for the most central collisions, which correspond to an energy density range of 2-3 GeV/fm^. Hitherto, 
using standard in-medium hadronic tools, this "anomalous" suppression is inexplicable. However, an 



explanation can be found in the transition to a QGP; e.g., Refs. [ 245 ,246], and this signal has recently 
been claimed ||247|| as ". . . evidence for the creation of a new state of matter in Pb-Pb collisions at the 
CERN SPS." 

Low-mass Dilepton Enhancement. Leptons produced in a relativistic heavy ion collision escape the 
interaction region without attenuation by strong interactions, which means they are a probe of phe- 
nomena extant in the early phase of the collision. Lepton pair production data collected in relativistic 
S-Au and Pb-Au collisions at the CERN SpS exhibits an excess with-respect-to proton-nucleus data in 



the "low-mass" region: 0.25 GeV < Me+e- < 0.70 GeV, which is illustrated in Fig. 4.10 



Since the clear excess appears above the tttt threshold and below the position of the p-meson peak 
in proton-nucleus collisions, the search for an explanation has focused on exploring the in-medium 
properties of the p-meson. Adequate explanations are found in: 1) coUisional broadening; i.e., the p has 
a shorter lifetime in an hadron-rich medium, and an increase in the p-meson's width due to in-medium- 



modified hadron-loop contributions to the p-meson's self-energy ||239| , |249|| ; and 2) a simple reduction 
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Fig. 4.9. Ratios of J/\[^-to-Drell-Yan cross- sect ions as a function of the average nuclear 
path: L, traversed by the cc-pair, which is a mo del- dependent but calculated quantity. The 
straight line is the scaling relation of Eq. ( |4.2.1D . (Adapted from Ref. ||243|| .) 



in the p-meson mass 
[T ~ 0.15 GeV] and baryon density [p 



IS. 



The phenomenological models require large values of temperature 
' 1-2 po] to describe the data. On this domain the p-meson 
screening masses, calculated in Refs. ||226| ^ 
by < 25% and, while Tp^^ 



and described in connection with Fig. O, are increased 

These constraints and 



is roughly unchanged, Fp^^^r is much reduced ||226 



effects are ignored in contemporary analyses, even though such features can affect photon production 

rates 



251 . 



Strangeness Enhancement. The lifetime of a terrestrial QGP can only be of the order: r ~ 5-10 fm, 
which is much too short for weak interactions to be important. Hence strangeness, once produced, can 
only disappear through s-s annihilation. Such events are unlikely unless there is a super-abundance 
of strangeness. Therefore strangeness carrying reaction products in the debris are a good probe of the 
conditions created by a relativistic heavy ion collision. 



The amount of strangeness produced in collisions can be quantified via the ratio ||252 

AA _ 2(S + S) 



A 



{u + u) + {d + d) 



(4.2.2) 



where {s + s), etc., are the mean multiplicities of newly produced valence quark-antiquark pairs at 
primary hadron level before resonance decays. Experimentally, as illustrated in Fig. [4.11| , 



A 



AA 



(4.2.3) 



In addition to this global enhancement, specific enhancements in the yields of K, K, A, A, etc., have 
been observed in CERN experiments [p^54 |. Detailed analyses indicate that quark degrees-of-freedom 
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Fig. 4.10. Recent Pb-beam data collected by the CERES group. Two independent 
data sets: 1995 and 1996, are compared to the dilepton yield expected in proton-nucleus 
collisions. The excess is characterised by an enhancement factor: 3.3, 1996; 2.5, 1995. It is 



an interesting effect but not a direct signal of QGP formation. (Adapted from Ref. |p48|| .) 



are necessary to describe the strangeness enhancement, with the suggestion that it may be indicative 



of QGP formation p4T . 



5 Toward a Kinetic Description 

Hitherto we have described features of cold, sparse QCD: the nature of DCSB, confinement and bound 
state properties, and then the effect that the intensive thermodynamic parameters chemical potential 
and temperature have on these phenomena. In the latter we applied the methods of equilibrium sta- 
tistical field theory and elucidated properties of the QGP phase. The terrestrial creation of this QGP, 
however, is expected to be effected via relativistic heavy ion collisions, which initially yield a quantum 
system far from equilibrium. This system must then evolve to form the plasma, and the study of that 
evolution and the signals that characterise the process are an important contemporary aspect of QGP 
research. 



5.1 Preliminaries 

The energy density in an ideal gas of eight gluons and two flavours of massless quarks is 



(2 X 8) + (2x3x2) 



30 



120 



30 



(5.1.1) 



As we saw in connection with Table |3.1| , the critical temperature for QGP formation is Tc ~ 0.15 GeV 
and at this temperature: egj^u+d = O.SGeV/fm'^. Construction of RHIC at Brookhaven National Lab- 
oratory is complete and it will soon provide counter-circulating, colliding 100 A GeV ^^^Au beams to 
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Fig. 4.11. By the measure in Eq. (|4.2.2| ), global strangeness production is enhanced by 
a factor of ~ 2 in A-B [S-S, S-Ag, Pb-Pb] collisions cf. e~^-e~ , p-p, p-p. The data are from 
Refs. [|53l. (Adapted from Ref. p52| ). 



generate a total centre-of-mass energy of 40 TeV, which corresponds to an initial energy density: e ~ 10- 
lOOGeV/fm^. The Large Hadron Collider [LHC] project at CERN is scheduled for completion in 2005. 
Plans are for it to provide ^°^Pb-^°^Pb collisions with y/s > 2 000 TeV and a consequent initial energy 
density e > 1 000 GeV/fm^. With these energy scales RHIC and LHC should certainly provide the 
conditions necessary for QGP formation. 

Control over the conditions produced in a relativistic heavy ion collision can only be exerted via two 
experimental parameters: the beam/target properties and, to some extent, the impact parameter. 
They can be used together to analyse the debris collected in the detectors. It is in the behaviour of 
this debris, summed over many events, that signals of the evolution and formation of a QGP must be 
found. Predicting just what the signals are requires an understanding of these processes, including and 
perhaps especially their the non-equilibrium aspects. 



Figure 5.1 illustrates the spacetime evolution of a relativistic heavy ion collision. All stages subsequent 



to the formation of a thermalised QGP are adequately addressed using hydrodynamical models |P56| 



However, in this approach the plasma's initial conditions [energy density, temperature, etc.] and EOS 
must be specified. These initial conditions can only be reliably determined by following the complete 
evolution of the system produced in the collisions; i.e., by understanding the pre-equilibrium stage. 
Furthermore, the very existence of a pre-equilibrium phase can lead to signals of QGP formation; e.g., 
plasma oscillations, a disoriented chiral condensate ||257|| , and out-of-equilibrium photon and dilepton 
emission 



The QGP is a hot, equilibrated, many-parton agglomeration, and in recent years two main approaches 
have been employed in describing how such a system might be produced in a relativistic heavy ion 
collision. In the perturbative parton picture |f^59|| , the colliding nuclei are visuahsed as pre-formed clouds 
of quarks and gluons. The collision proceeds via rapid, multiple, short-range parton-parton interactions, 
which generate entropy and transverse energy in a cascade-like process. In the string picture ||260|| , 
after passing through one another, the colliding nuclei are imagined to stretch a high energy-density 
fliix tube between them, which decays via a nonperturbative particle-antiparticle production process. 
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Fig. 5.1. Spacetime evolution of a relativistic heavy ion collision. Two [Lorentz contracted] 
disc-like nuclei collide and pass through each other, transforming a large amount of their 
kinetic energy into the potential energy of a highly excited spacetime volume. This energy is 
dissipated via particle production, which is a non-equilibrium process. The system evolves to 
equilibrium, expanding and passing through the QGP phase, hadronisation and freeze-out, 
where the observable hadronic debris is created. (Adapted from Ref. ||255|.) 



Each approach has its merits and limitations, and their simultaneous pursuit provides complimentary 
results. The analysis of data proceeds via one of the many Monte-Carlo event generators that have 
been developed for both pictures [ p6iy262| , |263|| . 

Flux Tube Model and Schwinger Mechanism. One intuitively appealing, semi-classical picture of con- 
finement is provided by the notion that it is effected via the formation of a small-diameter colour 
flux tube between colour sources; and there is evidence in lattice simulations for the appearance of 
such flux tubes between heavy-quark-antiquark pairs; e.g., Ref. ||264|| . This is a motivation for the 
string-like models just introduced, which have been used to study particle production in e~-e~^, p-p and 
p-nucleus collisions [26C, 265|| . A flux tube yields a linearly-rising, confining quark-antiquark potential: 



Vqq{r) = ar, where the string tension: a, can be estimated in static-quark lattice simulations. An over- 
stretched flux tube can be viewed as a strong background field. As such it destabilises the QCD vacuum, 
which is corrected by particle-antiparticle production via a process akin to the Schwinger mechanism in 
QED ||266|J267|| . We use this particle production mechanism as the primary medium for our discourse. 



The fermion production rate for a constant, homogeneous electric field E in a flux tube is 



dN 



dtdVd^p^ 



\eE\ In 



1 + exp 



27: {m^ +p\) 



\eE\ 



(5.1.2) 



where m is the mass and e the charge of the produced particles, and it is plain from this formula 
that the rate increases with increasing E and is suppressed for large m and/or p^. Particle-antiparticle 
production via this mechanism is analogous to a tunnelling process in quantum mechanics, as illustrated 
in Fig. f).2\ [NB. The Schwinger formula |266| follows immediately by integrating Eq. ( |5.1.2D over the 
transverse momentum: / d'^p±] 

What happens once the particles are produced? Naturally, they are accelerated in the tube by the field. 
This produces currents, which generate an electric field that works to screen the flux-tube [background] 
field. In the absence of other effects, the net field vanishes and particle production stops. However, 
the currents persist, now generating a field and renewed particle production that opposes their own 
existence. That continues until the net current vanishes. But at this point there is again a strong 
electric field . . . and the process repeats itself. This is the "back-reaction" phenomenon and the obvious 
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Fig. 5.2. Left Panel: For E = the vacuum is characterised by a completely filled 
negative-energy Dirac sea and an unoccupied positive-energy continuum, separated by a 
gap: 2eT = 2{rn? +p\Y^'^. Right Panel: Introducing a constant external field: E = CxE, 
which is produced by a potential: A'^ = —Ex, tilts the energy levels. In this case a particle 
in the negative-energy sea will tunnel through the gap with a probability ~ exp(— vre^/eE'). 
Succeeding, it will be accelerated by the field in the — x-direction, while the hole it leaves 
behind will be accelerated in the opposite direction. The energy-level distortion is increased 
with increasing E and hence so is the tunnelling probability. [NB. Ttn ~ ex/eE is a measure 
of the time between tunnelling events. eE can be related to the fiux-tube string-tension.] 




consequence is time-dependent fields and currents; i.e., plasma oscillations characteristic of the system. 
In recent years its affects have been studied in both boson and fermion production [ p68|| and we exemplify 



the process in Sec. |573 . 



There are a number of shortcomings in extant applications of the flux-tube particle-production picture 
to QGP formation: 1) The non-Abelian nature of the chromoelectric fleld is often ignored because the 
QCD analogue of the Schwinger mechanism is poorly understood. Instead the flux tube is represented 
via a classical electromagnetic potential; 2) On the whole, flnite volume effects are neglected, with the 
electric fleld assumed to be homogeneous in space. A more realistic description would account for the 
geometry of the interaction region. Some steps have been taken in this direction, with a consideration of 
effects produced by a cylindrical boundary [p69|| and those in a flnite-length flux tube with a conflning 
transverse potential p7(J|| . Going further, the geometrical representation of a flux tube can be replaced, 
allowing the flux tube profile itself to be a dynamical quantity, whose presence and stability is affected 
by the charged particle currents | ^71|| ; 3) The time-evolution of the system is described using either 
mean-field theory, which retains quantum effects but makes problematic an exploration of the effects 
of collisions, or a transport equation ||272|] , which neglects quantum effects. Below we will describe 
a partial reconciliation of these approaches; and 4) Little attention has been paid to particular non- 
Abelian features in transport equations |p73|| . A study [p74|| of Wong's equation [ p75|| is one step in the 



direction of explicitly including colour algebra effects. 

Quantum field theory can be applied directly to out-of-equilibrium plasma phenomena; for example, 
Refs. ||268| , |276|| . However, kinetic equations provide an appealing alternative because of their intuitive 
character. This approach begins with a transport equation 



dl_ 



-Qp'F. 



dpu 



S{p,q) + C{p,q) 



(5.1.3) 



where: / is the single particle distribution function, which gives the ensemble fraction of particles in a 
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given phase-space cell [p is four-momentum and g is a spacetime four- vector] ; F^y is the field strength 
tensor; and Q is the particles' charge. The source term: S{p,q), describes the production of particle- 
antiparticle pairs and C [p, q) is the collision term, whose origin is intuitively obvious but which is 
difficult to approximate accurately. 

A Boltzmann equation of this form has been widely applied to particle production using a source term 
of the type in Eq. ( p.l.2| ) both with and without a collision term, and in the following we exemplify this. 
However, to provide an introductory overview, Refs. | 272| | employ a classical, Abelian electric field in the 
source term but completely ignore the effect of collisions, C = in Eq. ( p. 1.31 ). The effect of collisions, 
represented via a "relaxation time approximation," [described in connection with Eq. ( |5.4.2|) ] has been 
considered in an hydrodynamical approximation to Eq. ( |5.1.3| ), Refs [ ^77|j278[| , with the influence of 
back-reactions neglected in the former but explored in the latter. A comparison between the transport 



equation and mean-field theory approaches has also been made ||268|| and the results are remarkably 
similar. That, however, is problematic since; e.g., the application of the Schwinger source term in the 
presence of a rapidly changing electric field is a priori unjustified and, although quantum field theory 
with its manifest microscopic time reversal invariance must underly the behaviour all quantum systems, 
experience confirms that systems far from equilibrium exhibit macroscopically irreversible behaviour 
that is amenable to treatment using [inherently time-irreversible] kinetic theory. The nature of the 
connection has recently been established P79| , P8(]| , ^81|| . 

At this point we re-emphasise that flux-tube models describe the nonperturbative production of soft 
partons. The production of hard and semi-hard partons is described by perturbative QCD and that 
mechanism is explored in Refs. ||282|| . Simultaneously incorporating both types of particle production is 
challenging but Ref. | 283 | is a step in that direction. Therein hard and semihard partons are produced 
via "minijet gluons" and provide the initial conditions necessary to solve the transport equations, and 
the subsequent evolution of the plasma is described by a classical but non- Abelian transport equation. 
Collisions are accounted for using a relaxation time approximation but quantum effects in the source 
term are neglected and, as will become apparent, they can be important in strong fields. 



5.2 Quantum Vlasov Equation 

The derivation of a quantum Vlasov equation in Refs. ||279|p80| , |281|| provides a connection between 
the quantum field theoretical and transport equation approaches to particle production and plasma 
evolution, and shows that the particle source term is intrinsically non-local in time; i.e., non-Markovian. 
Therefore calculating the plasma's properties at any given instant requires a complete knowledge of the 
history of the formation process. 

The derivation begins with the Dirac [Klein- Gordon] equation for fermions [bosons] in an external, 
time-dependent, spatially homogeneous vector potential: A^, in Coulomb gauge: Aq = 0, taken to 
define the 2;-axis: A = (0, 0, A{t)). The corresponding electric field is 

E{t) = -A{t) = , (5.2.1) 

also along the 2;-axis. The vacuum instability created by this external field is corrected via particle- 
antiparticle production. Fig. ^.2| , which is a time-dependent process. The transition from the in-state 
to the instantaneous, quasi-particle state at time t is achieved by a time- dependent Bogoliubov trans- 
formation, which effects the diagonalisation of the system's Hamiltonian. By this means one obtains a 
kinetic equation for the single particle distribution function 

f{P,t) = {0\al{t)ap{t)\0), (5.2.2) 
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which is defined as the vacuum expectation value, in the time- dependent basis, of creation and anni- 
hilation operators: a^()f:), ap{t), for particle states at time t with three-momentum P. That equation 
is 

^^iP^^^MM + eE(t)^|^ f_Jt'W^{t')[l±2U{P,t')]cos[x{t\t)], (5.2.3) 

where the lower [upper] sign in Eq. (|5.2.3| ) corresponds to fermion [boson] pair creation. The momen- 
tum is defined as P = {pi,p2, P\\{t)), with the longitudinal [kinetic] momentum P\\(t) = p\\ — eA(t), 
P|l = P3. [NB. eE(t) = P\\it), the particle velocity attained via acceleration by the field E{t).] For 
fermions [|279| , ^8(][] and bosons [ ^8(J| , ^81[ 1 the transition amplitudes are 



where the transverse energy e± = \Jim? + p]_, p± = (^1,^2), and uj{t) = ^e]_ + P'^{t) is the total energy. 
In Eq. (m), 

x(t\t) = 2[e{t) - e{t')] , e{t) = f dt'uit') , (5.2.5) 

is the dynamical phase difference. Equation ( ^.2.3|) is the precise analogue of directly solving QED 



with an external field in mean-field approximation, as done; e.g., in Ref. ||268|| . The physical content is 
therefore equivalent and, in particular, the fundamental quantum character is preserved. The appeal 
of this kinetic equation, however, is that it simplifies: the identification of reliable approximations; 
widespread applications; and numerical studies. 

This quantum Vlasov equation has three qualitatively important features: 1) The source is non- 
Markovian for two reasons - (i) the source term on the r.h.s. requires complete knowledge of the 
distribution function's evolution from t_oo — > t\ and (ii) the integrand is a non-local function of time, 
which is apparent in the coherent phase oscillation term cos x{t\t) and refiects the quantum nature 
of the source term; 2) Particles are produced with a momentum distribution cf. the Schwinger source 
term, which produces particles with zero longitudinal momentum; and 3) the production rate is affected 
by the particles' statistics, as evident in Eqs. (|5.2.4D and also in the sign appearing in the statistical 
factor [1 ± 2/-I-], which leads to different phase space occupation [p84|| . 



These features can have a material impact on the solution of the kinetic equation, and their importance 
depends on the field strength and the time-scales characterising the production process. The first 
time-scale is set by the Compton wavelength of the produced particles: 

Tqu ~ — . (5.2.6) 

Events taking place over times less-than this expose the negative-energy elements in particle wave- 
packets, a core quantum field theoretical feature. This is the time-scale of the rapid oscillations generated 
by the dynamical phase difference in cosx{t' ,t). The high frequencies involved mean that the main 
contributions arise when t ~ t', and therefore a local approximation can be justified for weak field 



plasmas |p81j| . We identified a second time-scale in Fig. p.2| ; i.e., the time taken for a particle to tunnel 
through the barrier. It can also be motivated by considering the transition amplitudes in Eqs. (|5.2.4|) , 
which for a constant electric field can be written in the form; e.g., 

^- = 77 / t-'if, / ^.,2 ' (5-2-7) 

[t - p\\/eEY + [e^/eEy 
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which is a Lorentzian characterised by the time-scale [width] 

n„~^. (5.2.8) 
This is also the time required to accelerate a charged particle to the speed of light in an electric field. 



The times scales in Eqs. ( p.2.6D and ( |5.2.8| ) are comparable when eE ~ e^. For weak fields, Tqu <^ Ttu 



and, in the integrand on the r.h.s. of Eq. ( |5.2.3|) , the cos x{t',t) oscillations occur on a time-scale 
so-much shorter than variations in the other elements that a stationary phase approximation is valid. 
That approximation yields a local source term and this important limit of Eq. ( |5.2.3[ ) was discussed in 



Ref. 281 



Ttu ^ Tqu for strong fields and hence the particle's Compton wavelength extends across the gap. In 
this case wave-like [quantum] interference effects become important in the behaviour of f{p,q), which 
changes at a rate comparable with cosx{t' ,t), so that a stationary phase approximation is not valid. 
Can this strong field scenario be relevant to contemporary relativistic heavy ion collisions? It is easy 
to make an estimate. The field strength in a fiux tube is commensurate with the QCD string tension; 
i.e., eE ~ a ~ 0.4 GeV^ ~ (2.5Aqcd)^, and e± ~ 2-3Aqcd is achievable. Therefore the answer is 
yes: effects arising from the non-Markovian structure of the source term are likely to be exhibited in 
contemporary relativistic heavy ion collisions. 



This was illustrated in Ref. ||285|| , wherein a comparison was made between the complete solution and 
that obtained in the low density limit; i.e., retaining cos x{t',t) in the source term but making the 
replacement [1 ± 2f±\ 1, which is the assumption 

/-t ~ , almost everywhere, (5.2.9) 

so that the source term becomes 

Sl{p,t) = \W±{t) f dt' W±{t') cos x{t',t). (5.2.10) 



Using Eq. (|5.2.iq) , the solution of Eq. (|5X3D is 



fl{m= dt'Sl{p,t'), (5.2.11) 

J —oo 

and is depicted for fermions in Fig. [5.3|. Particle production begins at t = because, with p\\ = 0, the 



transition amplitudes in Eq. ( |5.2.4|) are maximal at this point [the vector potential is zero], and the 
distribution function rapidly approaches a Schwinger-like asymptotic value: 

f-'Ut^oo)-exv(-ne^ /eE) - ""'^'^ 0.35 (strong) 

/ ^t^ooj-exp^ T^ejet.)- ^ 0.0019 (weak) " ^^'^-^^^ 



The qualitative features of the results depicted in Fig. are easily understood. A fermion, once 
produced with a certain momentum, retains it because there are no further interactions [collisions are 
ignored]. The number density is greater for strong fields because they produce more particles and this 
plainly means that the low-density limit will smoothly become invalid with increasing field strength. The 
low-density limit overestimates the fermion distribution function in strong fields because it eliminates 



Pauli blocking. The opposite effect is seen for bosons ||285|| . For very strong fields: /_ 1, again 
because of Pauli blocking, but there is no such bound on the boson distribution function. Finally, 
since ale\> is achievable in contemporary relativistic heavy ion collisions, the low-density limit is 
quantitatively unreliable. 
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Fig. 5.3. Time-evolution in a constant electric field of the complete single fermion 
distribution function cf. that obtained in the low-density limit, Eq. ( |5.2.9D , with P|j = = p± 
and initial condition f{t — oo) = 0. Upper Panel: Strong field, eEo/m? = 3.0; Lower 
Panel: Weak field, cEo/m'^ = 0.5. (Energy unit: m ^ Aqcd- Adapted from Ref. [|285|| .) 



5.3 Back-reactions 

The illustration above assumed a constant electric field. Allowing the more realistic case of a time- 
dependent field introduces another time-scale, namely that characterising the response of the field to 
the system's evolution. This brings us to the phenomenon of back-reactions, which have been explored 
in connection with models in cosmology and recently much in connection with QGP evolution. In both 
cases the particles produced by the strong background field modify that field: in cosmology it is the 
time-dependent gravitational field, which couples via the masses, and in a QGP, it is the chromoelectric 
field affected by the partons' colour charge. 

The effect of feedback is incorporated by solving Maxwell's equation, which for a spatially homogeneous 
but time dependent electric field is just 

m = -j{t). (5.3.1) 

The total field is a sum of two terms: an external field, Egx{t), excited by an external current, jex{t), 
such as might represent a heavy ion collision [that is a model input]; and an internal field, £'j„(t), 
generated by the internal current, jj„(t), which characterises the behaviour of the particles produced. 
The internal current has two components | |277|| : continued spontaneous production of charged particle 
pairs creates a polarisation current, jpoiit), that depends on the particle production rate, S{p,t); and 
the motion of the existing particles in the plasma generates a conduction current, jcondit), that depends 
on their momentum distribution, f{p,t). The internal field is therefore obtained from 



- Ein{t) = jin{t) = jcondit) + jpol{t) 



(5.3.2) 



In mean field approximation the currents can be obtained directly from the constraint of local energy- 
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Fig. 5.4. Time evolution for bosons of the electric field and current for an impulse external 
field, Eq. ( |5.3.7| ), with Aq = 10.0, b = 0.5 and the coupling = 5. The plasma oscillations 
driven by the back-reactions are evident. A complete explanation of the qualitative features 
evident here is given after Eq. ( p. 3. 7] ). (Energy unit: m ^ Aqcd- Adapted from Ref. 



density conservation: e = 0, where 

e(t) = ii?^(t) + 2 1 uip, t) fip, t) , 

and the factor of 2 accounts for antiparticles. For bosons the constraint yields 

. Pii - eA{t) 



-2e 



(2^ 



uj{p,t) 



uj{p,t) at 



and one easily identifies the currents 

P\\ - eA{t) 



jcondit^ 2e 



(27r)3 



uj{p,t) 



fiP,^), jpol{t) 



^u{p,t)S{p,t). 



(5.3.3) 



(5.3.4) 



(5.3.5) 



Equation ( 5.3.2 ) now yields the internal field. This construction has been used extensively to study 
back- reactions; e.g., Refs. 



The expression for the polarisation current exhibits the usual short-range divergence associated with 



charge renormalisation, however, its regularisation and renormalisation is straightforward |p68| , |28(: 
That accomplished. Maxwell's equation assumes the form 



A^{t) = E^{t) = -r{t) 



(5.3.6) 



(27r)3 



P\\{t) 



/.(P,t) + i 2."(^'^) eE-(t)P,Kt)W 



2 W±{P,t) 4u^{P,t) 



where = 2, = 1, and all fields and charges are understood to be fully renormalised. 

The effect of back-reactions can now be illustrated by solving this equation in concert with Eq. (|5.2.3|) . 
[Collisions are still neglected: C = in Eq. ( p. 1.3 ).] A relativistic heavy ion collision can be mimicked 
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Fig. 5.5. fip,t = fixed). The rapid oscillations are a quantum coherence effect, a 
manifestation of interference between effects on the Tqu and time-scales. Employing a 
local [or stationary phase or Markov] approximation to the source term, these effects are 
averaged-out. (Energy unit: m ~ Aqcd- Adapted from Ref. ||281|| .) 



by an impulse profile for the external field [p86|| 



Eex{t) 



An 



sech2(t/6) 



(5.3.7) 



which "switches-on" at t ~ —26 and off at t ~ 26, with a maximum magnitude of -Emax = Aq/I) at 
i = 0. Once this field has vanished only the induced internal field remains to create particles and affect 
their motion. The calculated field and current profiles are depicted in Fig. p.4| , and the qualitative 
features are easily understood. The impulse electric field is evident at t ~ 0. It produces particles and 
accelerates them, and their motion generates an internal current that opposes the impulse field. Shortly 
after the "collision" the current reaches a short-lived plateau, when the total field vanishes and particle 
production halts temporarily. At about this time the external, collision-mimicking field "turns-off." 
Nevertheless, in its absence, the total field grows in magnitude but now acts in the opposite direction, 
decelerating the existing particles, causing new particles to be produced and accelerating them in the 
new direction. The effect of that is clear, after a time the total current must vanish. A pattern is 
now established and, in the absence of other influences such as collisions or radiation, it repeats itself, 
reaching a steady state once the wash from the collision-mimicking impulse conflguration disappears 
completely. The oscillations characteristic of a plasma with fleld-current feedback have now set-in. The 
oscillation period is the new time scale: 



Tpi, the plasma period. 
It decreases with increasing field strength. 



(5.3.8) 



One more feature of the results is the high frequency oscillations evident at the current's peaks and 
troughs. They are not a numerical artefact and become more pronounced with increasing values of 



eE/m 



I.e. 



when Tqu ~ 



Ttu. This makes plain that they are a non-Markovian feature and result from 
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interference between effects on the tunneUing and quantum time-scales. Of course, as illustrated via an 
analogous feature in Fig. p75| , they disappear if a local approximation to the source term is used ||281 
because the stationary phase approximation suppresses such interference effects. 

One additional observation is important here. Our example employed an electric field whose magnitude, 
oc Aq, is large and hence the tunneling time, Ttu, is small, being inversely proportional to eE. The period 
of the plasma oscillations, Tpi, also decreases with increasing Aq but nevertheless, as clear in Fig. 



Ttu Tpi- Thus, in contrast to the effect it has on the production process [ 285 |, the temporal nonlocality 
of the non-Markovian source term is unimportant to the collective plasma oscillation | |286[| . This is the 
reason why kinetic equations with a simple source term of the form in Eq. (|5.1.2| ) are successful in 
describing plasma oscillations [ |268| |. Whether or not these oscillations are observable in relativistic 



heavy ion collisions depends on the effect of dissipative processes, which we now discuss. 



5.4 Collisions and Evolution of the Quark- Gluon Plasma 

Thus far we have illustrated the phenomena of quantum particle production in strong fields and field- 
current feedback. The stationary state is unrealistic because the dissipative processes: "collisions," that 
lead to thermalisation have been neglected. [Collisions can also lead to particle production and effect 
hadronisation but we defer that discussion.] 

In the presence of collisions the kinetic equation takes the form 

df±{p,t) 



dt 



S±{p,t) + C±{p,t). 



A simple and widely-used model |P78| , P871J288|J289| , |290|| for the colhsion term is 

fl'{p,t)-fAP,t) 



C±{p.t) 



where is the "relaxation time" and f^"^ is the thermal equilibrium distribution function 



exp[u;ip,t)/T{t)]Tl 



(5.4.1) 



(5.4.2) 



(5.4.3) 



with T{t) a time-dependent temperature profile, discussed on page The relaxation time is a fourth 
time-scale and it is plain that plasma oscillations can only be observable if Xpi -C t,.] i.e., if collisions 
take place much-less frequently than oscillations. 

In many of the exploratory calculations hitherto undertaken, Tj- is a constant. That might be argued 
to be inadequate because the collision time is supposed to characterise a system's thermalisation, a 
process which is interdependent with time-evolving quantities such as density and temperature. A 
more realistic representation might therefore employ a time-dependent Tr, which is calculated self- 
consistently as the plasma evolves. Reference ||291|| , employing collisions in a gluon- minijet plasma, is 
a step in that direction. However, from another perspective, any sophisticated collision term should be 



derived from, and justified by, an underlying microscopic theory (see; e.g., Refs. |P73| , |279| , [292|| .); and, 
furthermore, a relaxation time approximation of any sort can only be valid under conditions of quasi- 
equilibrium, which cannot be justified a priori in the presence of strong fields. A truly realistic collision 
term will introduce non-Markovian effects in addition to those already present in the source term, and 
exemplary studies exist in connection with: relativistic heavy ion collisions ||293|| ; collective effects in 



nuclear matter | 294 |; nuclear fragmentation | 295 |; and the damping rates of giant dipole resonances 
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Fig. 5.6. Time evolution for bosons of the electric field obtained using different relaxation 
times in the collision term of Eq. ( ^.4.2|) , and with the impulse external field, Eq. (|5.3.7|) , 
where Aq = 7.0, b = 0.5 and the coupling = 4. (Energy unit: m ^ Aqcd- Adapted from 
Ref. 



286| 



These studies make clear that even a binary collision approximation, as characteristic of a Boltzmann- 
Uehling-Uhlenbeck kinetic equation, can be inadequate under extreme conditions; e.g., in the presence 
of strong fields and/or when the particle density is high. These are precisely the conditions relevant 
to QGP formation. The patent complexity justifies the use in exploratory, illustrative studies of the 
simple Tj. = constant relaxation time approximation. Improvements are a contemporary challenge. 

As intuition suggests, collisions effect a damping of the distribution functions' time- dependent struc- 



ture. That is illustrated in Fig. ^.6| , with this particular calculation obtained using an Ansatz for the 
temperature profile ||286 



where T^g = A 



^QCD, 



T 

J. rr 



9T f"^ 

^ eq) ''0 



Tit) = T,, + (T^ - T,,) e-* /*o , (5.4.4) 

10/Aqcd- [The definition and meaning of T{t) is discussed below.] 
As evident in the figure, for large relaxation times the plasma oscillations are unaffected, as might be 
anticipated because this is the coUisionless limit. However, for ~ Xpi, the collision term has a significant 
effect, with both the amplitude and frequency of the plasma oscillations being damped. Finally, there 
is a value of Tr below which oscillations are not possible, just as in the case of an overdamped harmonic 
oscillator, and the system evolves quickly and directly to thermal equilibrium. The time taken by the 
plasma to thermalise depends on the ratio Tpi/rj., being longer for larger values. 



Why use the temperature profile in Eq. (|5.4.4| ) and, indeed, what is temperature in a system far 
from equilibrium? The notion of temperature is introduced via an assumption of quasi-equilibrium 
at each time t, which is only truly valid if the fields are not too strong; i.e., only as long as ^ 
Tpi. This temporally local temperature can be calculated self-consistently with the evolution of the 
distribution functions, an approach which represents a slight improvement over employing Ansdtze, 
such as Eq. ( ^.4.4|) , and, in fact, can provide an a posteriori justification for that expedient. The 
calculation of T{t) can, however, proceed in a number of ways so that the resulting profile is not unique 
but the differences are only small and quantitative. 



To illustrate the definition of a temperature profile we follow Ref. | |290| | . The total energy density in 



the evolving plasma is given in Eq. ( |5.3.3| ), where the second term is the particle contribution: e/(t). 
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Fig. 5.7. Time evolution of the fermion's energy density for RHIC (upper panel) and LHC 
(lower panel) conditions. For RHIC conditions, the result obtained using the low density 
approximation, Eq. ( [5. 2. 101 ), is indistinguishable from the complete solution. That is not the 



case for LHC-like conditions. (Adapted from Ref. [290 . 



and the particle number density is 

n{t) = 2j^Up,t). (5.4.5) 

An intuitive definition of the local temperature is to require that, at each t, the average energy-per- 
particle in the evolving plasma is that of a a quasi-equilibrium gas; i.e.. 



(5.4.6) 



n{t) n^i{t) ' 
where 

e^\t) = 2 coip, t)/5(p, t) , n^^(t) = 2 / ^g, t) . (5.4.7) 

With given in Eq. ( |5.4.3| ), Eq. ( |5.4.6| ) is an implicit equation for T(t), which must be solved simul- 
taneously with Eqs. and ( ^XTP [|9|. 



In Ref. | 290 |, this system of equations was solved with m = Aqcd, e = 2, and Tr = 1/Aqcd ~ 1 fm, and 
using two exemplary impulse field configurations: Eq. ( [5.3.7| ), with b = 0.5/Aqcd ~ 0.5 fm and either 
j^RHic ^ 4 Aqcd or Aj^^'^ = 20 Aqcd- As evident in Fig. |5]3, Aj^^^^ yields a calculated initial energy- 
density characteristic of RHIC conditions [see Sec. |5.1|] while A^^'-'' gives a very much greater value 
characteristic of the initial energy- density expected at LHC. For RHIC conditions the energy-density 
rises rapidly but, after reaching a maximum, decays monotonically. The low density approximation 
is valid here. For LHC conditions, with an initial energy- density twenty-times larger, the situation 
is different: the solution obtained in the low density limit is only a qualitative guide to the plasma's 
evolution and plasma oscillations are evident on observable time-scales. These oscillations retard the 
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time [fm] 

Fig. 5.8. Time evolution of the temporally local quasi-equilibrium temperature for RHIC 
and LHC initial conditions. For comparison we also plot the Ansatz in Eq. ( p.4.4| ), which is 
a fair guide for RHIC-like conditions. (Adapted from Ref. ||290|| .) 



equilibration of the plasma, roughly doubling the time taken cf. RHIC-like conditions. This application 
makes evident an expected terrestrial hierarchy of time-scales: 



Tqu ~ Ttu < Tpl ~ T, . 



(5.4.8) 



In a very intense field the produced particles travel almost at the speed of light and hence the plasma 
oscillation period can be estimated via the ultrarelativistic formula [^97] ]: 



2 V 



m 



+ ^2 



(5.4.9) 



using the input and LHC-like results of Ref. ||290|| , which yield rimax — 12/Aqqj-, from Eq. ( p.4.5|) . That 
this is a reliable guide for LHC-like conditions is evident in Fig. |^!^, remembering that e{t) ~ E'^{t) so 
that the peaks exhibited in the lower panel are separated by ^Tpi. 



The calculated temperature profile is depicted in Fig. |5^. Under RHIC-like conditions, an initial 
temperature of T{t = 0) ~ 0.5 GeV is reached and the temperature decreases monotonically. The 
LHC-like source conditions yield an initial temperature twice as large: T^^'~'{t = 0) ~ 0.9 GeV, and the 
temperature fluctuates in tune with the energy density. 



As described above, the exemplary impulse models in Ref. ||290|| set initial energy density scales appro- 
priate for RHIC and LHC. An improvement over such Ansdtze is to calculate the initial conditions, as 
done; e.g., in Refs. ||283|,p91|] for the case of a gluon minijet plasma. Even with this improvement, how- 



ever, the results obtained for observable quantities; such as, e{t), n{t), the temperature and equilibration 
time, are semi-quantitatively identical and all features are qualitatively unchanged. 

It is a uniform prediction that plasma oscillations are present under LHC-like conditions and they may 
manifest themselves in the dilepton spectrum. The dilepton production rate is (see; e.g., Ref. ||298|| ) 



dt d^x 



d'pid%f{Ei)f{E2)a{M) vu 



(5.4.10) 
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where / is the distribution function of the charged elements in the plasma, a{M) is the cross-section 
for particle-antiparticle annihilation into a dilepton pair with invariant mass M and f 12 is the standard 
phase-space factor [a "relative particle-antiparticle velocity"]. The distribution function responsible for 
the features in Fig. 5/7 itself exhibits structure and, according to Eq. ( ^.4.10 ), that will be transmitted 
to the dilepton spectrum. The oscillations characterising the pre-equilibrium phase of the plasma could 
then be exhibited as repeated dilepton bursts. However, this signal will only be detectable if it is strong 
enough relative to other process, such as Drell-Yan. That quantitative question is currently unanswered. 



5.5 Hadronisation 

The penultimate stage of QGP evolution is hadronisation; i.e., the re-establishment of the conditions 
that prevail at zero density and temperature. That process in not instantaneous and hence here too a 
dynamical approach is necessary, one that provides the means of describing phenomena such as bound 
state formation, and the onset of DCSB and confinement. Consequently the development of a transport 
theory in which such effects are directly accessible is an important current focus. 

Relativistic transport theory and hydrodynamics have long been useful tools in the study of non- 
equilibrium states of matter at nonzero density and temperature. They have been employed extensively 
in relativistic nuclear physics, both in the intermediate energy region |p92| , p99|| , where it may be appro- 
priate to neglect quark and gluon degrees of freedom, and in the high energy domain directly relevant to 
a non-equilibrium QGP ||262| , |300| , p01| , |302|| . One approach [p03|| to the derivation of Vlasov-like transport 
equations can be described as a relativistic generalisation of the Zubarev method ||304|| . Another is the 
contour Green function technique |P05|| , which has been applied formally; e.g., in Refs. ||300| , P01| , |302|| , 
and also in quantitative model studies; e.g., Refs. [|306| , |307i^08| , p09|j310|] . The latter is widely used 



because it admits a systematic definition and exploration of approximations, and will be our focus. 

To be concrete, extant phenomenological applications begin with a model Lagrangian density; e.g., 
that of the Friedberg-Lee or Nambu-Jona-Lasinio model. Such models are chosen because of their 
efficacy in describing DCSB, and have been employed in studies of: the spacetime evolution of the 
dressed-quark mass [ ^06|j307| , |308| , |309| | , which characterises the onset of DCSB; and hadronisation via 
the models' fermion-antifermion-bound-state interaction terms [piCI|| . In the former case, the onset of 
DCSB is explored by solving a quantum Vlasov equation coupled with a distribution-function-dependent 
gap equation, which describes the evolution of the particles' mass and feeds this information back into 
the system. One important qualitative result is a parton mass that increases as the system moves 
toward equilibrium; e.g., Ref. p07|| , which results in a dynamical softening of the EOS. This softening 
can have observable consequences; e.g., slowing the [hydrodynamical] expansion of the plasma |pll| . 
These models, however, do not incorporate confinement and hence cannot describe the replacement of 
coloured degrees of freedom by hadronic matter through the transition. 

An indication that the onset of confinement might significantly affect the evolution of the plasma was 
observed in Ref. ||309|| , wherein a T-dependent dressed-quark mass was introduced by fitting the model's 
quasi-particle energy density to that determined in lattice simulations. A self-consistent solution of the 
coupled Vlasov and gap equations then exhibits confinement; i.e., the partons cannot leave the QGP 
volume because their mass becomes infinite outside this region. The model does not include a parton- 
antiparton-hadron interaction and hence the confining effect also stops the hydrodynamical expansion. 
This study only admitted the importance of the scalar piece of the dressed-quark self energy. However, 
as we saw; e.g., in connection with Fig. |3.5| , the vector self energy is more important in QCD's deconfined 
domain. The effect it has on plasma evolution has not hitherto been explored. 
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The incorporation of such quahtatively robust features of DSE studies into the Vlasov equation is an 
obvious next step [P37|| . It is comphcated for [at least] two reasons: 1) The quasi-particle energy is a 
functional of the scalar and vector dressed-quark self energies, which are both nonzero and momentum 
dependent. The scalar self energy is small in the QGP. [It vanishes in the chiral limit.] However, the 
significant vector self energy remains | |17(j|| ; 2) In the confined domain the dressed-parton 2-point func- 
tions violate reflection positivity and hence a single particle distribution function for these excitations 
cannot be defined. That is as it should be because hadrons are the relevant degrees of freedom in this 
phase: a kinetic theory based on quarks and gluons is only appropriate in the deconfined QGP phase. 
As described in connection with Fig. p.3| , DSE models can adequately represent the transition from 
confined excitation to propagating quasiparticle ||181|| , an attribute exploited in Ref. ||176|| in calculat- 
ing the thermodynamic functions describing a QGP. These observations suggest that a truly realistic 
description of the quark and gluon distribution functions in a QGP will see them vanish in the vicinity 
of the phase boundary. 

We close this section with an illustration of the concepts introduced. In a QGP the evolution of the 
dressed-quark distribution function can be described using a Vlasov equation with the quasi-particle 
energy 

E,{p,x,t) =po[l + Ec{p,x,t)] = ^\p,\^ + B{p,x,ty, (5.5.1) 

where = p[l + x, t)] is the rescaled three-momentum and B{p, x, t) = [m + T^b{p, t)]. Using 

Eqs. ( p.l.lTp , Eq. ( |5.5.1[ ) can be rewritten 

E,{p, X, t) = C{p, X, t) E{p, X, t) = C{p, X, t) ^\p\''J{p,x,ty + M{p,x,t)\ (5.5.2) 

where I = A/C and M. = B / C . As our exemplar we follow Ref. p37|| and employ an instantaneous 
variant of the Pa model in Sec. 13.21: 



'Dip, 



3'K^r]6'\p). 



(5.5.3) 



With this interaction and using rainbow truncation, which is akin to a mean-field approximation, the 
Matsubara sum can be evaluated algebraically because the self energies are po-independent, and the 
gap equation assumes the form 



B{p,x,t) = m + r] ^^fl^^fl f^ [p, X, t) , C(p, x,t) = l + \ri ^^fl^^^^} f^p, x, t) 



{p, X, t) ' 



E^ (p, f , t) ' 



(5.5.4) 



with A = C, so that X = 1, and = f /C. [NB. / appears explicitly here because, with a quasiparticle 
pole, the Matsurbara sum can be evaluated algebraically.] Again, in spite of the model's simplicity, 
the solution exhibits qualitative features that are characteristic of a realistic dressed-quark 2-point 
function; e.g., momentum-dependent scalar and vector self energies, and the persistence of this aspect 
of the solutions in the deconfined domain | |237| 1 . In this class of models, for a coUisionless plasma, 
satisfies 

= dtfM^,t) (5.5.5) 
^ E{p X t) ^ + ^' t)^p-^(P^ ^' t)] ■ V^/,(p, X, t) - M{p, f, t)V^M{p, X, t) ■ Vp/,(p, x,t)} . 
It is important to note that the dressed-quark mass function is momentum-dependent so that in general 

VpM{p, X, t) ■ V. fM X, t) ^ . (5.5.6) 



[NB. The Vlasov equation in Ref. ||309[] is obtained by neglecting this term, since VpM = when the 
interaction is momentum- independent, and setting C = 1.] Equations ( |5.5.4| ) and ( p.5.5| ) provide a 
coupled system for the quark's distribution function and self energies. 
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As remarked above, the derivation of a transport equation like Eq. ( |5.5.5| ) requires the existence of 
a quasi-particle mass shell. The instantaneous model ensures that, with p° = E{p,x,t). In this case 
and = 0, hence Eq. ( |5.5.5| ) can be written 



dl*. 



(5.5.7) 



using the Minkowski space metric conventions of Ref, 
described; e.g., in Ref. ||302|| 



, which is directly comparable with the equation 
However, the momentum dependence of Ai, here driven explicitly via 



Eqs. ( |5.5.4|) , precludes a simple spherical space- volume scaling solution |P06| , pi2|| 



The coupled system, Eqs. ( ^.5.4|) , (|5.5.5| ), illustrates some of the complexity to be anticipated in studying 
the re-emergence of DCSB and confinement in an expanding QGP. Even in this simple model, as in 



Ref. 11309 1, proceeding further requires a numerical solution. Such studies are in their infancy but 
the qualitative feature is plain: the nontrivial propagation characteristics of the dressed-partons will 
significantly affect /, and hence QGP evolution, as T ^ Tc and the distribution of partons begins to 
resemble a heat bath. 



6 Epilogue 

Continuum strong QCD is a broad field and we have only presented a snapshot. Nevertheless, our 
discourse should serve to introduce many of the topics currently occupying practitioners. Since we 
employed the Dyson-Schwinger equations as our primary medium, it should also provide an update of 
the progress that has been made with this tool in the last decade and identify the current challenges. The 
text provides the details. However, highlights of the progress include: the application of a one-parameter 
model for the dressed-gluon 2-point function in a description of the masses, decays and form factors 
of light pseudoscalar and vector mesons, and its simultaneous prediction of the critical temperature 
for quark gluon plasma formation and the properties of in-plasma correlations; significant progress in 
calculating the baryon spectrum, and leptonic and nonleptonic interactions involving baryons; a unified 
treatment of chemical potential and temperature, and their effect on the equation of state for a quark 
gluon plasma and hadron properties; and an incipient understanding of the evolution to equilibrium of 
a quark gluon plasma, and the dynamical influence of the re-appearance of confinement and dynamical 
chiral symmetry breaking on the plasma's expansion and cooling. On the other side, a primary challenge 
is to comprehend the origin of the infrared enhancement in the kernel of the QCD gap equation that 
is necessary to ensure dynamical chiral symmetry breaking. In summarising we cannot improve on 
Ref. |I[]: "This programme has a long way to go, but [we] hope you are convinced it has come far." 

While the Dyson-Schwinger equations have provided the backbone for our discussion, we have made 
connections throughout with the results of other methods. Where there is agreement there can be little 
doubt that the phenomena described are real. That is the rationale underlying the simultaneous pursuit 
of complementary methods. Disagreement, in fact or interpretation, provides a challenge, which should 
be met, but also an opportunity for dialogue that should not be missed. 
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